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ABSTRACT 


Analytical  studies  have  been  made  concerning  many  of  the  problems 
associated  with  the  expansion  of  extremely  high  temperature  and  pres¬ 
sure  plasmas  through  the  converging  section  of  a  hypersonic  wind  tunnel 
nozzle.  The  problems  in  this  study  relate  to  the  establishment  of  heat 
fluxes  to  the  nozzle  wall;  methods  of  removing  the  heat  load  or  protect¬ 
ing  the  nozzle  wall;  determining  the  stress  levels  in  the  nozzle  liner; 
analysing  various  materials  for  strength  and  thermal  properties. 

Ii.  all  cases  parametric  studies  have  been  made  resulting  in  design 
criteria  which  can  be  used  for  specific  conditions.  The  conditions  are 
established  from  the  performance  envelope  for  a  typical  low  density 
hypersonic  wind  tunnel.  Throughout  the  report  the  performance  envelope 
is  plotted  on  stagnation  pressure  and  temperature  coordinates  which  cor¬ 
respond  to  free  stream  velocities  and  altitudes  up  to  25  kilofeet/sec  and 
350  kilofeet  respectively.  Ultimate  use  of  the  information  will  depend 
upon  the  reliability  of  the  assumptions  made  which  in  many  cases  must 
be  determined  by  new  experiments. 
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I.  INTRODUCTION 


/Analytical  studies  have  been  made  concerning  many  of  the  problems 
associated  with  the  expansion  of  extremely  high  temperature  and  pressure 
plasmas  through  the  converging  section  of  a  hypersonic  wind  tunnel  nozzle. 

The  problems  in  this  study  relate  to  the  establishment  of  heat  fluxes  to 
the  nozzle  wall;  methods  of  removing  the  heat  load  or  protecting  the  nozzle 
wall;  determining  the  stress  levels  in  the  nozzle  liner;  analysing  various 
materials  for  strength  and  thermal  properties. 

In  all  cases  parametric  studies  have  been  made  resulting  in  design 
criteria  which  can  be  used  for  specific  conditions.  The  conditions  are 
established  from  the  performance  envelope  for  a  typical  low  density  hyper¬ 
sonic  wind  tunnel.  Throughout  the  report  the  performance  envelope  is  plotted 
on  stagnation  pressure  and  temperature  coordinates  which  correspond  to  free 
stream  velocities  and  altitudes  up  to  25  kilofeet/sec  and  350  kilofeet 
respectively.  Ultimate  use  of  the  Information  will  depend  upon  the  relia¬ 
bility  of  the  assumptions  made  which  in  many  cases  must  be  determined  by 
new  experiments.  For  example,  curves  of  heat  flux  for  a  laminar  boundary 
layer  and  for  a  turbulent  boundary  layer  are  shown  -  each  excluding  a 
portion  of  the  phase  envelope.  If  it  is  assumed,  and  further  verified 
experimentally,  that  the  boundary  layer  is  laminar,  the  heat  flux  from  a 
large  percentage  of  the  phase  envelope  can  be  handled.  However,  there  is 
no  known  published  experimental  information  on  the  transition  in  converging 
sections  to  assist  in  checking  the  reliability  of  an  assumed  boundary  layer. 
Thus,  the  heat  flux  from  both  types  of  boundary  layers  is  included  in  the 
report. 

The  problems  in  this  report  have  been  investigated  and  reported  by 
various  investigators  associated  wi th  the  contract.  Specifically,  Dr. 

Stevens  has  made  the  analysis  on  film  and  transpiration  cooling,  and  on 
heat  flux  for  laminar  and  turbulent  boundary  layers;  Hr.  Roland  has  made 
the  analysis  on  structural  design,  materials,  and  backside  coolant  flow 
conditions;  Professor  Mott  is  investigating  the  problems  of  transition  from 
laminar  to  turbulent  boundary  layer;  Dr.  Pasqua  has  studied  the  radiation 
heat  flux  exchange  between  the  plasma  and  the  nozzle  walls;  Hr.  Robinson 
has  made  the  analysis  on  a  new  configuration  -  a  sharp  edged  nozzle. 


II.  HEAT  FLUX  FOR  LAMINAR  AND  TURBULENT  BOUNDARY  LAYERS 


When  backside  cooling  of  the  nozzle  liner  is  considered,  the  amount 
of  energy  that  can  be  removed  by  the  coolant  is  limited  by  the  backside 
heat  transfer  coefficient.  The  limiting  energy  removed  by  backside  cool¬ 
ing  is  approximately  50  Btu/in^sec.  Therefore,  in  order  to  maintain  the 
Integrity  of  the  nozzle  liner,  the  heat  flux  into  the  liner  must  not  exceed 
some  critical  value.  The  energy  addition  to  the  liner  depends  upon  many 
factors  and  is  controlled  by  the  type  of  boundary  layer  between  the  plasma 
and  the  wall.  A  laminar  boundary  layer  is  preferred  over  the  turbulent 
boundary  layer  as  a  lower  heat  flux  exists  for  the  former  when  the  same 


1 


AEDC-TDR.  63-58 


free  stream  and  wall  conditions  exist. 

Several  analyses  have  been  made  for  heat  transfer  assuming  a  laminar 
boundary  layer  (  Jjj  )*.  However,  they  do  not  include  all  of  the 

effects  encountered  in  the  present  investigation.  The  following  analysts 
is  made  to  establish  the  heat  transfer  coefficient  for  a  laminar  boundary 
layer  on  a  flat  plate  considering  the  effects  of  compressibility,  high 
temperatures,  and  variable  fluid  properties  In  the  boundary  layer.  This 
flat  plate  model  will  provide  a  good  estimate  of  the  average  heat  flux  in 
the  converging  section  if  evaluated  using  flow  conditions  that  exist  at 
the  throat  of  the  nozzle.  The  flat  plate  model  also  serves  as  the  basis 
for  an  incremental  calculation  of  the  laminar  heat  transfer  coefficients 
along  the  converging  section  of  the  nozzle.  Briefly,  the  nozzle  is  divided 
into  a  series  of  flat  plates  each  subjected  to  the  local  fluid  properties 
of  the  nozzle.  The  flat  plate  equation  adjusted  for  variable  free-stream 
conditions  is  used  to  calculate  the  variation  of  the  heat  transfer  coefficient 
across  each  of  the  intervals. 


Flat  Plata  Modal 

Consider  the  following  sketch  and  definitions: 

P  =  density 
o  =  veloci ty 
T  -  temperature 
Ji  =  specific  enthalpy 
Cp  =  sped  f  ic  heat 
^  3  viscosi ty 
A  -  thermal  conductivity 
4)  *  hydrodynamic  boundary 
layer  thickness 
*  thermal  boundary  layer 
thickness  (Sr<^,) 

Subscripts: 

w  x  wal I 

ow  „  adiabatic  wall 
S  3  free  stream 

The  integral  momentum  equation  for  a  flat  plate  with  constant  free 
stream  velocity  is: 


^Underlined  numbers  in  parentheses  refer  to  similar  numbers  in  the 
bibliography. 
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which  can  ba  reduced  to 


a  *  &  [[  ('-a.)  a  i  +]  %L 


Define  a  new  variable  ^  such  that 

jf  dc 


Jy  .  £ 

d“ 


<  ■/* 


*.*<?•/,  o/  y  »  <T 


ds/  -  dst  d. (  ;  ~  J 

df  4r  Ps 


Therefore,  the  integral  momentum  equation  becomes 

*<£/ %""#l 


(  1  ) 


Assume  a  velocity  profile  of  the  form 


e?J 


■  o  *  d  .</  ct  •  c  •  o 


.t 

at  £  -  *  /-  6o f  *  e«f,' 


_ £.  c/o  «o=  6  t*-3*6l  o*  A  ^3e/ 

<4  <7*  ‘ 

+hos  d  ■  J_  f>  p  -  _/ 

-?-C  ' 
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The  velocity  profile  becomes 


and 


J, 

r  s  ,-jj  )jz  j>>  »  I!  «/, 

j  l  Oj'  c(f  (■  ZBo 


end 


^  )  «  Jl 


Substituting  equations  (  2  )  and  (  3  )  into  (  I  )  gives 
J>9 


or 


280  J  z  (  f>,  J^*'  5^ 

J.JJ*.  ill  £~  c/x 

/S  ui  A* 

<S,  =  4.64  (St  )2(<£~Jl  )* 

K  y  ps  <y,  y 


(  2  ) 

(  3  ) 


(  4  ) 


The  integral  energy  equation  for  the  same  model  is 


Introduce  the  variable  *1  and  express 


t*  -  Hi '  Hi  Hi 

Also  for  reasonably  low  wall  temperatures 


jl  ~  Cfi-r 


9ne/  <Jr  ) 
/w 


/  e/u  ) 


Then 


M/.h#  fig  *>■ 

o  J  * 
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where 


A  5 


"/* 


(*<*  ■(fr<t.) 


Rearranging 


k(&£l 


(  5  ) 


Assume  an  enthalpy  profile  as 


~'i*'  _  cl  *  6  O  +  co  ~+  ft 

A,-X~  ~  C  C 


*  .  .  •  *  *■  _ '  t 


at 


°*  t 


al  so 


A  m  &  ,  /  0*10/  0^4.  3  ^  ^  <»/  e£ 

which  leads  to 


where 


_  ±  £. _ i.  /_* ) 

"  2  d 

-T^ )  *•  —  f'*'* m.  -I  -  -l 

cr/j  /p»©  £  5  * 


(  6  ) 


f  5  ^ 

r* 

Thus,  using  equations  (  2  )  and  (  6  ) 

///-  ^4^  7J/  2T  _3.  / 

J  L  A,.^^  J  Us  c  *0  '  ' 

O 

The  integral  energy  equation  becomes 

(A ->V/5 «,*£(*/').  f£  (g) ^A*. 


5 


A  E  DC- T  DR. 63-58 


and  can  be  reduced  to 


Substituting  equation  (  4  )  for  <f,  into  equation  (  7  )  gives 

*  *•*'£]  - ' 

where 

Pr)  = 

/w 

then 

/*  7Stk^-['- 

Xq  «  unheated  starting  section 


Define  a  heat  transfer  coefficient  such  that 


or 


h  a  *  -2.  ) 

-2  J*;  (  ^  7 


(  9  ) 


Substitute  equations  (  4  )  and  (  8  )  into  equation  (  9  ),  assuming  x„  »  0, 
gives 
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where 


*•)*  - 


The  heet  transfer  coefficient  from  equation  (  10  )  is  used  to  determine 
the  heat  flux  for  low  free  stream  velocities  by 

(  II  ) 

however,  at  high  free  stream  velocities,  it  can  be  shown  that 


(  12  ) 


where 


13  ) 


Eckert  suggests  that  better  correlation  with  experimental  data  is  obtained 
if  the  fluid  properties  are  evaluated  at  a  reference  enthalpy 

,9 

_/•  *  *  O'  (  |1*  ) 


Equation  (  10  )  was  developed  assuming  the  following  conditions  or 
restrictions  for  the  flat  plate: 

1)  laminar  boundary  layer 

2)  constant  free  stream  conditions 

3)  variable  fluid  properties  in  the  boundary  layer 

4)  constant  wall  temperature  along  the  flat  plate. 

This  model  may  be  translated  to  a  nozzle  in  order  to  calculate  an 
average  heat  transfer  coefficient  for  the  convergent  section  by  assuming 
that  the  free  stream  conditions  of  the  flat  plate  model  are  equivalent  to 
the  nozzle  throat  conditions.  The  converging  section  of  the  nozzle  in  the 
following  sketch  and  the  flat  plate  model  at  x  »  xt  are  assumed  equiva* 
lent. 
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» 


It  will  be  shown  later  (by  comparison  with  an  incremental  calcu¬ 
lation)  that  this  flat  plate  model  yields  an  average  heat  transfer  coef¬ 
ficient  suitable  for  parametric  analyses. 


Incremental  Technique 

The  determination  of  the  local  heat  transfer  coefficient  in  a  nozzle 
of  the  type  considered  in  this  report  is  most  difficult  if  approached  by 
the  usual  techniques.  In  contrast,  the  Incremental  technique  presented 
here  is  basically  simple.  Many  difficult  problems  (theoretical  and 
calculational )  are  circumvented  by  a  direct  utilization  of  accepted  flat 
plate  results. 
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The  nozzle  Is  represented  as  a  series  of  m  incremental  flat  plates 
each  of  length  A*, ,  .  The  nozzle  free  stream  conditions  are  superimposed 
on  this  flat  plate  model. 


L 


AX. 


AX, 


1 


AX_  Thrtf* 

X _ __J _ ! _ L 


S/*/ton  q  t  JL  n  ft 


Let  us  then  assume  that  the  variation  in  heat  transfer  coefficient  across 
an  incremental  flat  plate  A3£„  can  be  represented  as 

Hm,  =r  /?„  + 

where  represents  the  total  change  in  the  flat  plate  heat  transfer 

coefficient  across  the  AX*  interval.  For  the  general  problem,  it  must  be 
assumed  that  variations  in  the  free  stream  flow  parameters  03  and  ft  over 
will  affect  aA^.  It  follows  that  A  h„  can  be  written  as 


Let  h  for  an  incompressible  laminar  boundary  layer  (small  tempera¬ 
ture  differences)  be  represented  by 

h  =  0.332  A  (&)*(#'/* 

Thus 
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The  fluid  properties  of  the  boundary  layer  could  be  assumed  as  con¬ 
stant  over  the  entire  heat  transfer  surface  (evaluated  at  some  reference 
temperature  T#)  or  as  constant  over  each  particular  Interval  AK*  and 
evaluated  at  some  local  reference  temperature  Tm*  . 

This  technique  requires  that  the  heat  transfer  coefficient  be  known 
at  some  position  Xo  as  an  initial  condition,  i.e.,  the  stagnation  point 
or  a  sharp  edge  consideration. 

A  nozzle  subject  to  the  following  conditions 

T®  «  4000  k,  Tw  *  800  k,  P°  *  34  atm. 

X  -  1.3,  isentropic  expansion 

was  selected  for  an  example  problem.  Calculations  of  the  local  heat  trans¬ 
fer  coefficient  were  made  using  the  incremental  technique.  Also,  estimates 
of  the  average  heat  transfer  coefficient  were  obtained  for  several  vari¬ 
ations  of  the  flat  plate  model.  Figure  I  shows  a  comparison  of  the  results 
of  these  calculations  with  the  predictions  of  Cohen  and  Reshotko  (  1  )  -  the 
solution  to  this  problem  is  given  in  considerable  detail  in  the  Cohen  and 
Reshotko  report. 

A  good  correlation  between  calculation  techniques  was  achieved  for 
the  sample  problem.  It  is  proposed  that  nozzle  heat  transfer  coefficients 
be  calculated  by  incremental  techniques  when  accurate  local  values  are 
required.  The  precise  nozzle  geometry  and  flow  conditions  are  necessary 
input  data  for  this  calculation.  Further,  the  extremely  simple-to-use  flat 
plate  model  should  provide  a  good  estimate  of  the  average  coefficient  over 
the  entire  converging  section.  This  calculation  is  independent  of  the 
throat  size  so  that  results  can  be  applied  (with  some  caution)  to  any  of 
the  considered  nozzles  for  a  given  flow  condition  if  the  result  is  modified 
by  the  simple  (xref/xt)^  ratio.  A  parametric  analysis  of  heat  flux  for 
a  xref  =  3"  based  on  the  flat  plate  model  is  shown  in  Figure  2  in  rela¬ 
tion  to  the  performance  envelop. 

The  magnitude  of  heat  transfer  from  the  plasma  main-stream  to  the 
nozzle  walls  is  strongly  influenced  by  the  nature  of  the  boundary  layer 
which  develops  along  the  nozzle  wall.  A  laminar  boundary  layer  will  in 
effect  insulate  the  nozzle  wall  from  the  hot  plasma.  In  contrast,  the 
transverse  fluid  motions  characteristic  of  a  turbulent  boundary  layer  pro¬ 
motes  extreme  rates  of  heat  transfer.  It  is  felt  that  an  accurate  estimate 
of  the  heat  transfer  under  turbulent  conditions  be  known  since  the  nature 
of  the  boundary  layer  cannot  be  predicated  with  confidence  for  all  nozzle 
configurations.  Following  usual  practice,  the  rates  of  turbulent  heat 
transfer  were  calculated  according  to  the  approximate  method  of  Bartz  (  2  ) . 

The  usual  approximate  method  of  solution  for  a  turbulent  boundary 
layer  involve  considerable  calculation  complexity  -  usually  requiring  an 
elaborate  computer  program.  The  Bartz  (  2  )  approximation  equation  for 
for  the  heat  transfer  coefficient  was  obtained  by  effecting  a  numerical 
fit  of  a  simple  dimensionless  equation  containing  the  characteristic  para¬ 
meters  of  the  system  with  the  results  of  the  usual  approximate  boundary 
layer  calculations 
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o  o 
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X 


Cohen  and  Reshotko,  NACA  1289,  1956 
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Incremental  Technique  (Sharp  Edge  Entrance) 
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o  w 
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Position,  x,  in. 


Fig.  1 


Comparison  of  Laminar  Hoot  Transfer  in  o  Convergent  Nozzle 
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where  CT  is  e  correction  factor. 


The  heat  transfer  coefficient  obtained  from  equation  (  15  )  is  based  on 
a  temperature  difference,  l.e. ,  qw  -  h  (Tw  -  Tw).  The  temperature 
difference  seriously  underestimates  the  energy  available  for  heat  trans¬ 
fer  from  the  disassociated  gas.  Experience  has  shown  (  3  )  that  a  heat 
transfer  coefficient  based  on  an  enthalpy  difference  (otherwise  evaluated 
in  the  usual  manner)  provides  a  reasonable  estimate  of  the  heat  transfer 
from  a  gas  with  disassoclation  effects  present.  Following  this  procedure, 
the  turbulent  heat  transfer  was  calculated  by 


n 

(  16  ) 

where  Hj  - 

^  ' 

lifv/tory  (/S) 

(  17  ) 

yL  »  local 

plasma  enthalpy 

The  calculation  of  the  turbulent  heat  transfer  was  included  as  a 
sub-routine  to  the  parametric  analysis  of  transpiration  and  film  cooling. 
Chapter  IV.  The  results  are  shown  in  graphical  form  for  the  three  nozzle 
sizes  under  consideration,  Figures  3,  *♦,  and  5.  It  should  be  noted  that 
the  practical  limit  of  50  B/in*sec  excludes  most  of  the  desired  performance 
envelope  if  the  boundary  layer  is  turbulent. 


Stagnation  Pressure,  ata 


Fig.  3  Convective  Heat  Load  (Throat);  Turbulent  Boundory  Layer; 
Twall  =  700nK  D* -  0.042  " 


AEDC-TDR-63-58 


III.  THERMAL  RADIATION 


The  problems  of  thermal  radiation  are  more  pronounced  in  the  entrance 
to  the  nozzle  than  in  the  throat  because  of  higher  temperatures  and  larger 
volumes  of  gases.  The  analysis  for  radiation  heat  flux  assumes  a  particu¬ 
lar  model  involving  an  emitting  and  absorbing  gas  enclosed  by  black,  non¬ 
radiating  walls.  The  assumption  of  black  walls  is  conservative  as  there 
could  be  a  reduction  in  absorbed  energy  at  the  walls  if  the  reflectivity 
was  high  and  gas  velocity  high.  However,  it  is  doubtful  whether  a  reflec¬ 
tive  surface  could  be  maintained  because  of  contamination.  The  walls  have 
been  assumed  non-radiating  because  of  their  relatively  low  temperatures. 

The  assumed  configuration  for  the  entrance  and  the  converging  walls 
of  the  nozzle  is  one  of  a  sphere  as  shown  in  Figure  6. 


Fig.  6  Modal  (or  Co*  Thermal  Radiation  Study 


The  emission  of  a  gray  gas  of  volume  V  to  V  +  dV,  radiating  uniformly  in 
all  directions  is 

c/g  _  4  ;U  J~  7"  oV  <j4io 

'  47T 

where  M  is  the  mean  absorption  coefficient.  The  volume  element  is  ex¬ 
pressed  as 


d \/  Z/T **  J \/r?  ft*  at'&c/r' 
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Assuming  exponential  absorption  of  emitted  energy,  the  amount  of  energy 
per  unit  solid  angle  absorbed  on  the  black  walls  becomes 

*  ir 


t-ff 


#  *•  *  ,  I  / 

C7*  T  /"  ^ 


c  o 


Noting  that 

Z%  =.  &  *■  r  -  ZrECo s <Z$ 


<dau  =  a//9  Ces<t> 

Z 


Co 5  *.  >€*-  reoz  4> 

Z 


z  jr  _^r  te'+r'-- xt-*****) 

2  ^  g-T( 3">  ^  g _ _  c/rc/4  (18) 


the  energy  absorbed  per  unit  wall  area  becomes 

^  JJ  (Gx>rx’-  zer  C030)  3/>* 

o  o 

In  order  to  integrate  equation  (  18  )  define 
A  =  yuu  f  E+r  -  zrK.  cos  <?) 

fi*  -i 


thus 


2  A  c//Z  =•  Zy*-  r  rS  a/4> 


/n  $  c/d  —  ^-eJX_ 
>fV€ 


then 


/  >i  )  *■  / 

zj2o~r V <*- ^r<cu */$  g'  o/A 

A.#'-/*) 
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and 


(e-  rcos  4)  = 


3  a  -JL C 

z  a*?.'  ze. 


resulting  In 


z  O+fl) 


ZX^  ZB  y 


or 


J-  • 

c rr« 


Integrating  the 


X 


variable  gives 


1 

<rr 


-s*-*Cn./s) 


(  19  ) 


or 


^Ul/Z  (TT  * 


(  20  ) 
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Equation  (  20  )  Is  evaluated  for  a  particular  value  of>T^  by  graphically 
Integrating  a  Y  versus  A  curve.  The  results  of  such  Integrations  give 
Figure  7. 

In  order  to  determine  the  mean  absorption  coefficient,  Information 
from  Kivel  and  Bailey  (  4  )  Is  used.  They  have  expressed  the  intensity  of 
radiation  in  wave  length  zJ  as 

3  —is 

JTV  -  2 yU*  l  o-  r  ^  J  o  &  o4z 


where 

Cy  H.  ^  •  of  do  h 

K  t  J  *z-r 


and  3--2J-  i  ntensi  ty  per  unit  wave  length. 

Thus 


o 
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integrating 

JZ*  =.  ZyU,  Z 


3  ! 

(&y 


JZ  ^  A3*  yZ  L  0~T 


Using 


4t- 

£  U~T 


the  mean  absorption  coefficient  becomes. 


^  * 


€ 

/.  e> *  l 


Values  of  total  emissivity  per  cm,  ~  ,  are  taken  from  the  report 
(  2  )  and  the  resulting  values  of  mean  absorption  coefficient  jui  are 
shown  in  Figure  8  as  related  to  temperature. 


Using  the  information  from  Figures  7  and, 8,  and  choosing  a  radius  of 
1.91  cm  (0.75  in),  enables  one  to  plot  f  YtAfi  as  a  function  of 
temperature  and  this  is  shown  in  Figure  9.  *  The  heat  flux  per  unit  area 


O^T*yL 


can  now  be  calculated  and  the  results  are  shown  in  Figure  10.  The  heat 
flux  information  is  also  shown  in  Figure  II. 


A  limiting  value  of  heat  flux  from  this  analysis  exists  when  the  gas 
does  not  absorb  any  of  the  emitted  radiation.  This  situation  occurs  when 
<®_>ta  in  equation  {  18  )  is  one,  i.e.^e^o  (note  that^Z  is  not 
considered  zero  in  the  emission  process).  Using  the  condition  that  eT** 
is  one,  equation  (  18  )  reduces  to 


(  21  ) 


which  is  the  heat  flux  per  unit  area  of  an  emitting  and  non-absorbing 
gas.  Results  from  equation  (  21  )  are  shown  in  Figure  10  which  shows  that 
the  heat  flux  is  Increased  for  the  non-absorbing  gas. 

The  analysis  for  thermal  radiation  heat  flux  assumes  that  the  con¬ 
vergent  walls  of  the  nozzle  see  a  spherical  volume  of  emitting  and  absorb¬ 
ing  gas;  the  gas  exhibits  and  exponential  absorption  characteristic;  the 
amount  of  energy  emitted  by  the  gas  is  proportional  to  the  mean  absorption 
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(19)  and  (20) 
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coefficient  of  the  gas;  and  that  tha  mean  absorption  coefficient  of  the 
gas  can  be  determined  by  integrating  tha  Intensity  par  unit  wave  length 
over  all  wave  lengths.  The  final  results  shown  in  Figure  12  indicate 
that  radiation  contributes  an  important  part  to  the  total  heat  flux  as 
tha  pressure  or  density  increases.  Tha  effects  of  thermal  radiation  can 
be  reduced  significantly  by  decreasing  the  radius  of  the  planum  or  stag¬ 
nation  region.  Since  thermal  radiation  would  be  the  major  contributor 
to  tha  heat  flux  In  the  plenum  region  it  becomes  a  necessity  to  decrease 
tha  radius  of  this  section  to  as  small  a  value  as  possible  without  affect¬ 
ing  other  parameters. 


IV.  TRANSPIRATION  AND  FILMCOOLINC  OP  HYPERSONIC  NOZZLES 


Laminar  boundary  layer  stability  considerations  may  restrict  the  use 
of  pure  'backside11  cooling  techniques  to  a  small  portion  of  tha  desired 
performance  envelope.  Other  cooling  techniques  -  "transpiration"  and 
"film"  •  may  allow  nozzle  operation  over  a  considerably  larger  portion  of 
the  envelope. 

Film  cooling  Is  accomplished  by  Introducing  the  coolant  at  soma  point 
upstream  of  tha  critical  heat  transfer  areas.  The  Joint  actions  of  the 
system  pressure  gradient  and  the  shear  stress  of  the  plasma  boundary  layer 
establishes  a  coolant  film  over  the  surface  to  be  protected.  A  schematic 
Illustrating  the  basic  features  of  film  cooling  Is  shown  In  Figure  1}.  The 
magnitude  of  tha  "apparent"  heet  flux  mey  be  decreesad  or  Increased  over 
that  with  backside  cooling  alone  because  of  the  Interactions  of  the  coolant 
with  the  plasma  boundary  layer.  The  Influence  of  this  effect  upon  the  heat 
load  Imposed  on  the  coolant  film  will  be  described  by  a  suitably  defined 
"blockage  factor."  The  blockage  factor  Is  defined  as  the  ratio  of  the  con¬ 
vective  heat  transfer  to  the  ft im  to  that  transferred  to  the  wal 1  when 
there  Is  no  film.  It  mey  be  necessary  to  augment  film  cool  Ingwl  th  back¬ 
side  cooling  to  ensure  stability  of  the  film  and  to  achieve  a  maximum 
cooling  effect. 

Transpiration  cooling  1s  similar  In  concept  to  film  cooling  except 
that  the  coolant  Is  Introduced  as  required  along  the  entire  cooled  surface, 
Figure  14,  The  energy  Is  first  transferred  from  the  plasma  to  the  porous 
non  Is  wall  end  then  to  the  coolant  as  It  flows  through  the  porous  wall. 

The  diffusion  out  from  the  wall  of  the  eoolant  vapor  acts  to  reduce  the 
magnitude  of  the  hast  flux,  Again  a  blockage  factor  will  be  used  to  quanti¬ 
tatively  describe  this  effect.  Thr  nature  of  tha  transpiration  process 
should  ensure  •  large  end  dependable  reduction  In  heet  flux,  backside 
cooling  will  not  In  general  be  compatible  with  transpiration  cooling. 

The  basic  mechanisms  of  film  and  transpl ration  cooling  can  be  repre¬ 
sented  by  a  single  generalised  model.  The  magnitudes  of  a  suitably  defined 
blockage  factor  and  a  backside  .cool  I ng  heat  flux  q*  can  be  adjusted  to 
provide  a  representation  of  alther  film  or  transpiration  eoollng  adequate 
for  parametric  analyses.  Tha  results  of  the  parametric  analyses  Indicate 
that  both  transpl ration  and  film  cooling  appear  to  offer  promise  of  extend- 


27 


Fig.  14  Sketch  of  Transpiration  Cooling 
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ing  the  envelope  of  nozzle  operation.  It  should  be  noted  that  an  actual 
design  Mould  require  a  more  precise  calculation  for  Its  basts. 


PARAMETRIC  ANALYSES 

The  flow  schematic  shown  In  Figure  15  represents  a  generalized  energy 
balance  for  a  differential  section  of  the  boundary  layer  region  -  either 
transpiration  or  film  cooling.  The  net  energy  transferred  to  the  coolant 
in  the  differential  element  is  carried  off  by  a  dm  amount  of  coolant  which 
experiences  both  an  enthalpy  increase  and  a  change  in  state.  The  blockage 
factor  V  is  so  defined  that  changes  In  heat  transfer  due  to  coolant  dif¬ 
fusion  into  the  plasma  and  increased  coolant  requirements  because  of  coolant 
loss  (directly  into  plasma  stream  before  experiencing  a  complete  enthalpy 
increase  as  a  film  coolant)  are  properly  accounted  for.  Consider  an  energy 
balance  on  the  element  of  coolant  film  shown  in  Figure  15. 


d 


rr\ 


where 


dA  = 
dm  - 
hui  = 


^9  = 

Aw  = 


Sc  - 

Sr  - 

- 

= 

<P  = 
9  - 


(  23  ) 


differential  nozzle-surface  area 

differential  coolant  requirement 

local  turbulent  heat  transfer  coefficient  based 

on  an  enthalpy  (plasma)  difference 

local  enthalpy  of  coolant  vapor  evaluated  at 

the  local  saturation  temperature 

initial  coolant  enthalpy 

local  stagnation  enthalpy  (plasma) 

enthalpy  (plasma)  corresponding  to  the  local 

saturation  temperature  of  the  coolant 

local  turbulent  convective  heat  flux  (without 

film) 

local  radiative  heat  flux 

local  heat  flux  removed  by  backside  cooling 

local  heat  flux  from  the  film  to  the  wall 

transmissivity  of  the  water  film 

"apparent"  local  convective  heat  flux  to  the  f iM «n 

1  oca  I  turbulent  convective  heat  flux  to  the  wall  (wi th- 

out  film) 
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The  total  coolant  requirement  Is  obtained  by  integrating  equation  (  22  ) 
over  the  entire  nozzle  area 


nn. 


h9~  Ve 


(  24  ) 


Equation  (  2k  )  was  evaluated  for  various  combinations  of  stagnation 
temperatures  (  T®  )  and  stagnation  pressures  (  P®  )  for  the  three  nozzle 
sizes  under  consideration  (  0  =  0.042,  0.272,  0.784  in.)  by  use  of  an 
IBM  1620  digital  computer.  The  program  includes  the  calculation  of  the 
plasma  flow  which  is  uniquely  determined  by  a  T° ,  P® ,  0  combination  for 
the  assumed  isentropic  (frozen)  expansion  through  the  nozzle.  The  cool¬ 
ant  requirements  were  calculated  as  a  coolant-flow  to  plasma-flow  ratio. 
These  results  are  presented  in  parametric  form  In  Figures  16  through  27. 

It  should  be  noted  that  results  are  not  shown  for  system  stagnation 
pressures  in  excess  of  ISO  atmospheres.  Thermodynamic  plasma  properties 
at  higher  pressures  could  not  be  found  in  the  literature.  The  calculations 
can  easily  be  extended  to  the  higher  pressures  when  the  data  becomes  availa¬ 
ble.  The  data  contained  in  these  figures  should  be  regarded  as  intermediate 
data  which  require  careful  interpretation  in  their  use. 


The  parametric  results  for  pure  transpiration  cooling  (  qr  =  0  )  can 
be  extended  to  include  radiation  heat  loads  and  to  predict  film  cooling 
requirements  by  addition  of  the  value  of  the  second  integral  of  equation 
(  24  )  to 


no  - 


rri 


Again  let  us  simplify  the  integration  (conservatively)  by  replacing  the 
local  values  of  qw  and  hR  by  their  nozzle  throat  values  and  qr  by 
qr(T#  P°  D° ) 

4r  (T#  P°  0®)  =  upstream  radiation  load  corresponding  to  the 

stagnation  temperature  T® ,  stagnation  pressure  P® , 
and  entrance  diameter  D® . 


Then 


nn  =  no.  + 

rntns. 


[  frCrW)  - 

hfc. 


(  25  ) 
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Fig.  17  Flow  Ratios  for  Transpiration  Cooling;  Blocking  Foctor  =  0.5; 
D*  -  0.784  " 


34 


Stagnation  Pressure,  ata 


Fig.  19  Blocking  Factors  for  Tronspirotion  Cooling;  M  (coo!ont)/M  (gas)  *  0.1; 
D*  -  0,784" 


I 


Stagnation  Pressure,  atH 
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Fig,  21  Flew  Ratios  for  Transpiration  Cooling;  Blocking  Foctor  »  0.5; 
D*  -  0.272" 
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Stagnation  Pressure,  ata 
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Fig.  23  Blocking  Foctori  for  Transpiration  Cooling;  M  (coolant)  M  (got)  »  0.1; 
0*  -  0.272 " 
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Fig.  25  Flow  Ratio*  for  Transpiration  Cooling;  Blocking  Factor  «  0.5 
0*  -  0.042" 
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Fig.  26  Flow  Ratios  for  Transpiration  Cooling;  Blocking  Foctor  =  0.1; 
D*  -  0.042” 
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The  magnitude  of  q  is  not  necessarily  determined  by  the  available  back 
side  cooling,  recall  equation  (  2 3  ) 

iL  —  Swc.  *  ^ 

JL  (conservative  approximation) 


(26) 

The  permissible  values  of  qwc  are  restricted  to  relatively  low  values 
because  of  the  poor  wall  heat  flux  convective  mechanisms  which  would  exist 
in  the  coolant  film.  Stronger  mechanisms  within  the  film,  i.e.,  turbulent 
exchanges,  would  most  likely  destroy  the  film.  The  removal  of  qw£  by 
backside  cooling  should  Improve  the  stability  of  the  coolant  film  -  the 
loss  of  which  would  certainly  result  in  burnout.  Also,  under  high  radiation 
loads,  q*  will  comprise  primarily  0q_  which  must  be  removed  by  backside 
cooling  (film  cooling  only)  because  of  the  limited  energy  convection  mecha¬ 
nisms  which  exist  within  the  film.  Without  backside  cooling,  large  radi¬ 
ative  loads  would  raise  the  temperature  of  the  nozzle  walls  until  the  film 
literally  explodes  away  from  the  wall  with  subsequent  burnout.  It  should 
be  emphasized  that  transpiration  cooling  can  handle  all  radiative  and  con¬ 
vective  loads  without  backside  cooling  within  the  practical  limits  of 
pumping  coolant  through  the  porous  walls  and  the  compati bl 1 1 ty  of  the  en¬ 
suing  coolant  vapor  with  the  plasma. 


SUMMARY 

The  studies  to  date  have  demonstrated  that  either  transpiration  cool¬ 
ing  or  film  cooling  will  extend  the  envelope  of  nozzle  operation.  However, 
It  should  be  noted  that  the  full  extent  of  the  new  operating  envelopes  and 
the  relative  merits  of  the  two  concepts  have  not  been  fully  established. 
Many  questions  are  still  unanswered  -  this  information  would  be  required 
for  any  nozzle  design  involving  either  concept.  Problems  requiring  further 
study  include: 

1)  determination  of  the  blockage  factors  ^ranj.  and  ^  film 

2)  backside  cooling  requirements  using  film  cooling 

3)  mechanisms  of  handling  the  radiative  load 

k)  interaction  of  the  coolant  vapor  and  plasma  -  special  design 
problems 

5)  effect  of  injected  coolant  on  nozzle  dimensions 

6)  determination  of  the  possible  "burnout"  conditions. 
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V.  NOZZLE  LINER  STRESS  ANALYSIS 


There  ere  three  cases  of  interest  in  the  stress  analysis  of  a  nozzle 
throat  liner,  namely,  (I)  an  unrestrained  liner  (2)  a  partially  restrained 
liner  (3)  a  fully  restrained  liner. 

By  an  unrestrained  liner  It  is  meant  that  the  liner  is  mounted  rigidly 
at  one  end  to  an  entrance  plenum  or  an  exit  expander  but  is  not  restrained 
by  a  sleeve  or  ribs.  It  Is  thus  free  to  expand  long! tudinal ly  and  radially 
except  for  self  restraint. 

By  a  partially  restrained  liner  it  is  meant  that  the  liner  is  mounted 
rigidly  at  one  or  both  ends  and, in  addition,  it  is  restrained  radially  and 
possibly  longitudinally  also,  by  a  surrounding  sleeve.  It  is  coupled  to 
this  sleeve  by  ribs  between  the  liner  and  sleeve.  The  sleeve  is  thin  enough 
that  it  expands  some  as  the  liner  expands  and  thus  the  liner  expansion  is 
only  partially  restrained. 

By  a  fully  restrained  liner  it  is  meant  that  the  radial  and  longi¬ 
tudinal  expansion  of  the  liner  are  completely  constrained  by  the  ribs, 
however,  bending  of  the  liner  between  ribs  is  possible.  In  the  present 
analysis  this  case  may  then  be  studied  by  increasing  the  constraint  in  case 
(2)  to  the  limit. 

There  are  two  fundamental  areas  of  stress  analysis,  namely,  elastic 
and  plastic.  Until  a  material  is  stressed  past  the  yield  point  it  is  con¬ 
sidered  to  follow  Hooke's  Law  (stress  proportional  to  strain),  and  the  analy¬ 
sis  for  simple  geometries  is  straight  forward.  When  the  yield  point  is 
passed  the  stresses  are  redistributed  by  yielding.  The  magnitude  of  the 
stress  actually  developed  then  depends  on  the  stress-strain  relationship  of 
the  material  past  the  yield  point.  The  problem,  so  far  as  the  nozzle  throat 
liner  is  concerned,  is  then  one  of  determining  if  the  yield  point  will  be 
exceeded  and  if  it  will  the  number  of  cycles  of  operation  before  failure 
must  be  predicted. 

To  determine  if  the  material  will  be  stressed  past  the  yield  point  for 
a  given  set  of  conditions  the  elastic  analysis  must  be  made.  The  stresses 
thus  calculated  are  compared  with  some  selected  yield  criterion  to  determine 
if  yielding  will  occur. 

For  the  elastic  analysis  the  geometry  was  simplified  somewhat.  Except 
for  the  thrust  load  the  liner  and  sleeve  were  considered  to  be  concentric 
cylinders  separated  by  a  coolant  annulus  and  ribs.  The  yield  criterion 
selected  was  the  critical  octahedral  shear  stress  (_§_). 

In  cases  where  the  yield  point  of  the  material  was  exceeded  the  number 
of  cycles  to  failure  was  predicted  using  the  equation  developed  by  Coffin 


/V 


A 


CortS-AlnJ 
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a/  ■  number  of  thermal  cycles  to  failure 
m  plastic  yield,  or  plastic  strain 
*  constant  which  varies  from  0.25  -  0.75. 


ELASTIC  ANALYSIS 
Cose  1  •  Unrestrained  Liner 


A  thick  walled  cylinder  is  assumed  in  calculating  hoop  and  thermal 
stresses  in  this  analysis.  Tr.e  following  general  formulae  for  thick  walled 
cylinders  are  available: 


CTr  »  &  Mi  -  t>x +  ( fo-t*) 

t>u-  a  r'-(b'-a') 

Ur  3  tangential  stress  due  to  pressure  (at  inside  or 
outside  radius  depending  on  choice  of  ),  psi 
CL  =  inside  radius,  inches 
b  «  outside  radius,  inches 
/«U'  a  inside  pressure,  psi 
P*  a  outside  pressure,  psi 
r  *  C  o.  Oy-  =  stress  at  inside  radius 

J  &  f*r  C~T  -  stress  at  outside  radius 

-  <=**  AT  //-  Zb'L  /n  h) 

/r*a-  zO-xs)/nb.  1  b'a1- 


(  27  ) 


(  28  ) 


e  at-  /  /_  z  a 1 


29  ) 


Orr  =  tangential  thermal  stress  from  self-restraint,  psi 
e*  *  thermal  expansion  coefficient,  in/ln  ®F 
/J  a  Poisson's  ratio 
b  a  outside  radius,  inches 
O  a  inside  radius,  inches 
£  a  modulus  of  elasticity,  psi 
at  =  temperature  difference  across  liner,  ®F. 

Calculations  of  the  airstream  pressure  variation  along  the  nozzle  axis 
show  that  the  pressure  remains  almost  constant  until  a  very  short  distance 
from  the  minimum  throat  diameter.  It  then  drops  rapidly  at  this  point  and 
continues  a  gradual  decrease  in  the  expanding  section.  An  approximate 
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thrust  load  on  the  nozzle  may  thus  be  calculated  using  the  stagnation  air 
pressure  and  the  projected  liner  area  in  the  direction  of  the  axis. 

Three  end  conditions  for  this  case  may  be  considered.  In  the  first, 
the  nozzle  is  fastened  at  the  exit  end  and  free  to  expand  thermally  at 
the  entrance.  In  the  second,  the  nozzle  Is  fastened  at  the  entrance  and 
free  to  expand  thermally  at  the  exit,  in  the  third  condition,  the  nozzle 
is  fastened  at  both  ends  and  thus  does  not  have  free  axial  expansion. 

Since  there  must  be  smooth  aerodynamic  matching  of  the  nozzle  with  the 
following  expansion  section  in  a  wind  tunnel,  the  second  condition  may  be 
eliminated.  Considering  the  third  condition,  it  is  noted  that  In  the  case 
of  the  smallest  throat  diameter,  since  the  expansion  forces  must  balance, 
there  would  be  very  high  compressive  stresses  at  the  throat.  This  would 
be  reduced  some  by  bending,  due  to  the  geometry,  but  this  would  change  the 
throat  diameter  and  probably  produce  buckling  at  the  small  throat.  For 
these  reasons  only  condition  one  is  analyzed. 


Fig.  28  Sketch  of  Linor,  Rib,  and  Sloovo 
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(a)  Compressive  thrust  stress  at  the  throat  (approximate): 


<J1 

{fit- n>A\  -  r"  ' 

Pr 

CfiC-fi,') 

P* 

a 

ai rside  pressure 

X 

coolant  pressure 

fi 

X 

plenum  inside  radius 

fit 

X 

throat  inside  radius 

fix 

X 

throat  outside  radius 

(b)  The  hoop  stresses  are  given  by: 


^  c, \  -  kpl.  .  ffi*-  ^  fix 

Vjj,»  m  - - V-  - 

r  -  r  r  *-/>  x  rx) 

Vz  fit  rtt  (.fix  -fit  / 


(  30  ) 


(  31  ) 


(Note  that  subscript  11  refers  to  oosition  at  the  inside  radius  and  12  at 
the  outside  radius;  t  refers  to  the  airside  and  b  refers  to  the  coolant 
side.  See  Figure  28). 


’bUZ 


fii  -  fix  Pb  ,  ( fit  - Px  )  fit  fix. 


f  _  /<■ 
fit  'u 


r,xC^- 


letting 


t  x 

-A-  ^ 

/-  '  /-  *- 

'Z  '»! 


and 


* 


fix 


r  _  r 

/a.  '// 


= 


/X  *t 


then  C 1/^ /f  —  fit  —  %  ^x  fib 


(  32  ) 

(  33  ) 

(  34  ) 

(  35  ) 

(  36  ) 
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•nd  C£/Z  js  Z  -  2,  p* 


(  37  ) 


(c)  Thermal  stresses  due  to  self-restraint  are  given  by: 


®thu 


^TMU  - 


oi  £  AT 


2(/-jU)  tn  '2k 

n, 


oi  £  AT 


Z(i-M)ln  Hh. 


/"J£ ) 


(  38  ) 

(  39  ) 


(These  apply  both  longitudinally  and  tangentially). 

(d)  Radial  pressure  stresses  on  the  airside  and  backside  are: 


Oiy,  =  -  /=V  (  *0  ) 

~  ~  hi,  (  **’  > 


The  summation  of  the  appropriate  stresses  from  above  gives  the  principle 
stresses  at  the  inside  and  outside  of  the  liner. 

(e)  The  principle  tangential  stresses  are: 

Inside  radius  - 


0Z..  =  2  ft  -  Z  E,  to  +  «-^T: - (/-Z  /n  )  (  42  ) 


Outside  radius  - 


O^.  B  ZZ^-Z.b  t-  ?-£  A T~-  e ?3  ^ -g*-)  <  ) 

"  zo-^)  In  (  n,  / 


50 


r 

i 

i 

1  AE  DC- TDR. 63-58 

f 

(f)  The  principle  longitudlnel  stresses  ere: 

Inside  redius  - 

azLII  .  +  JLE2L —  (/-  2^1  i»  z±)  (  i*  ) 

(  nS  -  O  In  2> 

'll 

Outside  redius  - 


.  -Cn-KXQ- ntf  +  *££Z  . (-  »,  '»  2±)  <  «  > 


(g)  The  principle  redial  stresses  ere: 
Inside  radius  - 


o?„,  .  -  f>,  <  ■*  > 


Outside  radius 

»  -  A*  (  «  > 


The  octahedral  shear  stresses  are  then  calculated  using  the  general 
formula: 

_____  ^  _  'll. 

°°CT  =  i  i  49 ) 

where  CTj^  CT^  CT^  are  the  principle  stresses  at  the  point  considered. 

Co*o  2  >  Partially  Restrained  Linar  •  Longitudinal  Rib* 

As  in  Case  I,  it  is  assumed  that  the  liner  and  the  surrounding  sleeve 
are  thick  walled  cylinders  in  this  case  for  restraint  by  longitudinal  ribs. 
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For  restraint  by  cl rcumferential  rtbs  it  is  assumed  the  liner  is  a  thin 
welled  cylinder,  as  the  solution  for  this  case  was  readily  available  and 
the  approximation  is  good  except  at  very  small  liner  diameters. 

The  same  notation  is  used  In  this  case  as  In  Case  I.  When  new  terms 
are  Introduced  they  are  defined.  As  in  Case  I,  the  formulae  for  hoop  and 
tharmal  stressas  for  thick  walled  cylinders  (equations  27,  28,  and  29)  are 
used. 


(a)  For  calculation  of  liner  and  sleeve  tangential  stresses  assume 
a  pseudo-pressure  P|  between  the  ribs  and  the  liner  and  ?2  between  the 
ribs  and  the  sleeve 


(  50  ) 


p2  and  P)  are  thus  equivalent  uniformly  distributed  loads.  (Subscript  22 
refers  to  the  inside  of  the  sleeve). 


77- 'Tz  «  2<Tlt7rrrt 

S>z  -  ^zx. 


<f,e.  unit  strain  of  outside  of  liner,  in/in 
unit  strain  of  inside  of  sleeve,  in/ln 


Since 


CTUii 


(  51  ) 

(  52  ) 


(  53  ) 

(  5**  ) 


(Subscript  I  refers  to  the  liner,  2  to  the  sleeve). 


now 

in 


%  -(?**•  ?>)  r>z 


K<\ 


L  rH 


%  -  +  G)J 

p  *■_  r  x- 

riz  rn 


(  55  ) 
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*  Z,  pt  -  2  «,.(  Pbi-P.)  (  57  ) 

in  the  same  manner 

«.  -  e,  (  Pb*  R)  (58) 


Since  the  outside  pressure  (gauge)  on  the  sleeve  is  zero 


<>«  -  <'***'  ) 


(  59  ) 


where  a  outside  sleeve  radius 


Letting  £.  *  /*'t  *  *** — and  substituting  from  equation  (50) 

Sa*  ~ 

CW  =  *  Pt  ^  )  (  60  ) 

Substituting  from  equation  (  5*+  ) 

0*1  *+  Pb  i-  *<-  P  (  61  ) 

E,  <\*  — 

Equating  07t ’ s  in  equations  (  58  )  and  (  59  ) 

Z  ?3  -  *,(Pb+P,)~  —  ^  ^  —  P, 

Solving  for  Pj 

R.  -  pb(g. ) 

*?.  — 


(  62  ) 


(  63  ) 
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substituting  from  equation  (  50  ) 


Px  * 


*± S  p-- 

2, 


*■ 


€, 


(  64  ) 


Substitutions  into  equations  (  57  ) .  (  58  )  and  (  60  )  give  G7t  ,  07j_ 
O’i'K,  respectively.  ‘ 

(b)  The  hoop  and  longitudinal  thermal  stresses,  due  to  self-restraint, 

are: 


e(8  6T 

zc '-»*■)  /» 

'a.'*. 


^TH>1.  s 


at  6,  ar 
Z  (/-***.)  U 

r^x. 


('-**>  u!t) 


(  65  ) 

(  66  ) 


(c)  The  hoop  stresses  in  (b)  above  are  modified  by  the  force  between 
the  liner  and  the  ribs  caused  by  the  thermal  expansion  of  the  liner.  Assum¬ 
ing  pseudo-pressures  P^  between  the  ribs  and  liner  and  Pj  between  the 
ribs  and  sleeve 


^nth  -  ^thu  —  ^  (  67  ) 

^NTIT  m  ®TWIX  _  Pe.  (  68  ) 

*♦  PS  {  69  ) 


s.  net  thermal  hoop  stress 
Since 


TV-  H  (  70  ) 

rtx 

<W  -  *4  (  71  ) 

rl.l_ 
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but 


cn 


urn 


v'**- 


(  72  ) 


Substituting  from  equations  (  65  ) ,  (  66  )  and  (  68  ) 


(  73  ) 


Solving  for 


P*-- 


(  74  ) 


Substituting  from  equation  (  74  )  into  equations  (  67  ) ,  (  68  ) 
and  (  71  ) 


rr  _  <t~  ??  / 

^Hr  II*  UT#»-  <,rt,[ - - - 

W,  ♦  f-4. 

( . - i - 'j 

^  V. 

r  & _ ' 

\  2.  a  JL_  / 

•  4  €. 


^UTIZ  “  ^THIL 


(  75  ) 


(  76  ) 


(  77  ) 


(d)  Beam  stresses  are  calculated  from  the  general  formulae: 
(1)  Over  a  rib 


*■*  ¥(±) 


(  78  ) 
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»  bending  tensile  stress,  psi 
AP  »  load  per  unit  area,  ps  1 
L  *  distance  between  supports,  inches 
t  *  beam  depth,  inches 

(2)  Midway  between  ribs 


+1 

ll 

b* 

(  79  ) 

4  V  t  / 

P.+P*. 

(for  liner) 

(  80  ) 

APS--  P*i-P<r 

(for  sleeve) 

(  81  ) 

then  inside  the  liner  over  a  rib 


(  32  ) 


Dr  s  distance  between  ribs,  inches 


outside  the  liner  over  a  rib 


‘ -Ir^-hUr) 


(  83  ) 


inside  the  liner  midway  between  ribs 


(  84  ) 


outside  liner  between  ribs 


(  85  ) 
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(e)  The  radial  pressure  stresses  at  the  inside  end  outside  of  the 
I iner  are: 

(t)  Over  a  rib,  inside  and  outside  respectively 


-  "Be  (  86  ) 

*■  “  )  (  87  ) 

■  ri b  ci rcumferentlal  thickness,  inches 
(2)  Between  ribs,  inside  and  outside  respectively 


^Rl!2.  * 

-"Re 

(  88  ) 

- 

-"Pe, 

(  89  ) 

(f)  If  the  liner  is  constrained  longitudinally  by  the  sleeve  a 
uniform  compressive  stress  will  be  introduced  in  the  liner.  In  this  case: 


4  *  ^ 


now 


&Tm 


S,  =.  y.  ^tlL 

rise  in  liner  mean  temperature 

Sx  -- 

E  , 


from  a  force  balance 


7T(^ 

-  ** 7r'(^~  n')<rTu 


(  90  ) 

(  91  ) 


(  92  ) 


(  93  ) 
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from  equations  (  90  ) ,  (91  )  and  (  92  ) 


^  ^  e,  +  <r;Li  ±_ 


(  94  ) 


substituting  for  from  equation  (  93  ) 


^u- 


-°i,  e,  -  07u 


ii.  /  r'*~-  ' 

v  ^sV-  V, 


(  95  ) 


let 


2jl  -<r. 

r  r"  ^ 

*1  rll 


(  96  ) 


solving  for  ^TU  t  in  equation  (  95  ) 


-  - 


«*,  £T, 

/*■?<•  4-' 


(  97  ) 


substituting  into  equation  (  93  ) 


C^T-LT.  ~ 


<*■  A-r^  g,  z,, 

/-A- 


(  98  ) 


(g)  The  pri  nciple  tangential  stresses  are: 
(I)  At  the  inside  radius  over  a  rib 


^PCill  »  >  ^urii  ^ 


(  99  ) 


(2)  At  the  outside  radius  over  a  rib 


(  100  ) 
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(3)  At  the  intide  rediut  between  ribs 


<^PC.  III  * 


^MTII 


»»r 


(  101  ) 


(4)  At  the  outside  radius  between  ribs 

^PetiT.  -  +■  <^isv*.7_  (  102  ) 

(h)  The  principle  longitudinal  stresses  are: 

(1)  At  the  inside  radius 

^u.  --  ->■  l  '03  ) 

(2)  At  the  outside  radius 

-  ^TU  ^TH  |  z.  (  104  ) 


(i)  The  principle  radial  stresses  are  given  in  (e). 


Cot*  2  -  Partially  Rattrainad  Liner  -  Circumferential  Ribs 


The  stresses  for  circumferential  ribs  were  first  analyzed  using  the 
pseudo-pressure  method  used  for  longitudinal  ribs.  A  more  exact  analysis 
is  that  of  Timoshenko  (  £  )  for  an  infinite  cylinder  with  ciroumferential 
restraining  bands.  This  latter  analysis  was  used  for  the  final  calculations 
and  is  presented  below.  It  assumes  a  thin  walled  cylinder. 


(a)  Bending  stresses  at  ribs  and  midway  between  ribs  must  be  calcu¬ 
lated.  When  a  distributed  radial  force  is  applied  around  a  line  on  a 
cylinder  the  resistance  to  bending  is  analogous  to  that  of  a  beam  on  an 
elastic  foundation.  The  differential  equation  for  the  deflection  of  the 
strip  is 

D  djL 

y  »■  J 


(  105  ) 
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where 


D*  All. 


the  flexure)  rigidity  of  the  wall,  Ib/ln 


(  106  ) 


■  wall  thickness,  inches 
r  ■  cylinder  mean  radius,  Inches 
^  *  deflection  from  Initial  position,  inches 
X  *  distance  along  a  line  in  the  surface  parallel 
to  the  axis  of  the  cylinder,  inches 

Defining 


B  = 


(&S-  ( 


3-  ) 


>U 


(  107  ) 


the  general  solution  of  equation  (  105  )  is 


«*  *■  C  8xJ  G  { e Cos  3k  +■  £  Sin  *X 


.ax 


(  108  ) 


A,  C,  E,  and  F  are  constants  to  be  determined  from  known  deflections, 
slopes,  etc.  at  various  values  of  x. 

For  the  case  of  a  single  concentrated  load  P  lb.  acting  at  x  =  0  on 
a  beam  of  unit  width  and  thickness  t  inches  (part  of  the  cylinder  of 
radius  r) , 


U  SB.  C.  ( Cos  &TC  +  &X  ) 

A7  =■  -  —  e  ( 3/n  fly  —  Cos  f3x  3 

4-C> 


(  109  ) 


(  HO  ) 


where  M 


bending  moment  (lb  in  ) 


thus 
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crs=  1  1  _P_  s.  (  S  m  B  x  -  ccs  "&*.  ) 
Z  K  ' 


(  H2  ) 


(When  x  «  0  the  quantity  in  brackets  is  negative,  so  the  positive  sign 
gives  the  stress  at  the  outside  of  the  cylinder). 

Since  the  general  differential  equation  is  linear  the  bending  moment 
for  multiple  loads  may  be  calculated  by  the  principle  of  superposition, 
and  thus,  the  stresses  may  be  calculated  by  superposition,  it  is  necessary 
to  take  only  a  few  loads  to  get  the  approximate  result  for  an  infinite 
beam,  as  the  terms  become  essentially  negligible  after  Bx  cs  7.0. 

The  general  solution  for  the  bending  stresses  may  now  be  written  as 
an  infinite  series  using  equation  (  112  ) 


z  at 


CO 


R.C  (Si* 


(3h£*|_  Cos  alt,, 


>1 


(  113  ) 


The  terms  for  negative  n  are  the  loads  and  their  distances  on  one  side 
of  the  point  chosen  as  x  *  0  for  the  stress  calculation.  The  terms  for 
positive  n  are  the  loads  and  distances  on  the  other  side.  For  n  >  0, 

x  -  0,  P0  is  the  load  at  x  «  0.  The  absolute  value  signs  apply  to  x 

since  it  can  be  shown  from  the  general  equation,  equation  (  108  ),  that 

the  contribution  from  a  load  P|  at  a  distance  xj  on  one  side  of  x  *  0 

is  the  same  in  sign  and  magnitude  as  that  of  the  same  load  at  the  same 
distance  on  the  other  side  of  x  »  0. 

The  deflection  may  also  be  written  as  an  infinite  series  using  the 
same  notation 


<r-  - 

J  zst 


to 

m 


Pn  C  ( <3/«  O/xJ  +  CCS  %/*«!)  I 


(  114  ) 


Equation  (  114  )  is  used  to  calculate  the  force  P  at  each  rib. 
First  assume  all  the  ribs  are  the  same,  so  Pn  -  Const.  »  P.  Then 

®  fl/v  / 

Hz  -CO 


L 


7 

8 /TCj  +  Ccs  <3/jr„/ ) 


(  115  ) 
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The  radial  expansion  of  tha  linar  dua  to  pressure  and  increase  in  mean 
temperature  (unrastralnad)  is 


Ar*  (3- +  r«  trm 

£t 


(  1 16  ) 


A  rib  with  cross  sactional  araa  Ar  and  Modulus  of  alasticlty  Er  is 
incraasad  in  radius  by  a  fcrca  P  (Ib/lnch  of  cl rcuMfaranca)  an  amount 


Ar 


Pr v 


Ar  £r 

Than  tha  dafiactlon  of  tha  cylindar  is 

U*  .  rocATrr,  -  _3rl 

Using  aquations  (  115  )  and  (  118  ),  and  letting  Ef  •  Ej, 

c*llU(s,in 4. cov *\xr,\)  »  A~3)r.  ErL 
Z€,t  Z_  E,t  *r<s, 

V\«-«0 


(  117  ) 


(  118  ) 


(  119  ) 


solving  for  P 


('t-PL)r'’ 

e,  * 


*■  r«AT-„ 


j3 

Z€, 


— —  ^[e~9IXJ(a/n  8/xJ  t  CCS 


(  120  ) 


J1 

Are, 


nm-co 


Multiplying  numerator  and  denominator  in  equation  (  120  )  by 


At 


P 


($-/%)+ 


~2  ^  s'"  *  Ccs  ^"O]  -J- 


(  121  ) 


n%- 
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For  the  case  of  P  -  Const. 


-Ofx«l 


_  -J  <  V 

=  t  *p  )  e  (sm  air* l  -  cox  a(t*| ) 
zm  *■  l— 


(  122  ) 


n«-« 


Since  P  is  calculated  at  a  rib,  with  constant  distances  between 
ribs  (dj )  equetion  (  121  )  may  be  written 


P*. 


_ cp*-> %)+  _ 

+  ^  »  Utijt  *■  c*>3  *  ~ 

^  <*«/  ^ 

Equation  (  1 22  )  may  be  written,  for  the  bending  stresses  at  a  rib 


(  123  ) 


=  - - -  2  y  «  ('omn/sJJ  -  Cos  r>  i  (124 

* at  L  »t;  J 


°S.*«  •  ~  *T,., 


(  125  ) 


For  bending  stresses  midway  between  ribs 


0»i» 


126  ) 


*  ~  <^isin. 


(  >27  ) 


(b)  Hoop  stresses  under  and  between  the  ribs  are  calculated  using 
the  deflection  equation  (  H5  ).  The  general  hoop  stress  formula  for  a 
thin  wailed  cylinder  is: 


A  « 

r  e. 


(  128  ) 
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SO 


cr 

tfeoP 


(  <29  ) 


Now  A  is  the  deflection  of  the  cylindrice!  surface.  When  the  pressure 
difference  (  Pt  -  )  and  the  rise  in  mean  temperature  A  Tm  occur  the 
deflection  is 


A 


f*r 


(  130  ) 


and  the  hoop  stress  is 

cr  _ 

^T’  - * - 

The  ribs  cause  a  deflection  inward  of 


(  131  ) 


A 


ft 


(Z-Vr*- 

e.t 


+  r-t  Arm 


Pr * 

Are, 


(  132  ) 


and  a  stress  of 

<T.  =-  _  r^crmet  + 

*  * 


Pr 

Ar 


The  net  hoop  stress  under  the  ribs  is  thus 


(  133  ) 


*•  ^ - *Tm  *'  (  >34  ) 

*r 

The  deflection  caused  by  the  ribs  between  ribs  is 

x  00 

£e7~L C  +  Co*(l!~>)OJx']  (  135  ) 

ffce 
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The  net  hoop  stress  between  ribs  Is  thus 


(*-5k 


Z*  )04 


(  136  ) 


(c)  The  thermal  stresses  due  to  self-restraint  are,  for  a  thin  wall, 
both  circumferentially  and  longitudinally 


at  e, 

Z  ( /  -  ju  ) 


(  137  ) 


(d)  The  radial  pressure  stresses  are: 

(1)  Inside  radius  under  a  rib 


Bill 


*  _  Pu 


(2)  Outside  radius  under  a  rib 


tr  =  rib  thickness,  inches 


(3)  Inside  radius  between  ribs 


-  -  5 


(4)  Outside  radius  between  ribs 


<£/«  -  ~  P* 


(  138  ) 


(  139  ) 


(  140  ) 


(  141  ) 


(e)  The  calculation  of  longitudinal  stresses  from  the  thrust  load 
with  circumferential  ribs  is  not  so  straight  forward.  In  the  case  of  no 
restraint  the  entire  thrust  is  taken  by  the  liner.  In  the  case  of  longi¬ 
tudinal  ribs  essentially  all  the  thrust  Is  taken  by  the  sleeve  through 
the  ribs.  In  the  case  of  circumferential  ribs  the  thrust  component  of  the 
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pressure  which  will  be  taken  by  the  ribs  depends  on  the  number,  thickness, 
redial  width,  end  constraint,  and  angle  of  the  ribs.  If  the  ribs  are 
pieced  at  all  points  perpendicular  to  the  liner  they  will  take  the  thrust 
load  largely  in  compression,  but  if  they  are  placed  perpendicular  to  the 
flow  axis  they  will  take  the  thrust  load  mainly  in  bending. 

It  was  at  first  proposed  that  tha  ribs  be  placed  perpendicular  to 
the  liner  but  difficulties  In  fabrication  end  assembly  led  to  the  con¬ 
clusion  the  ribs  should  be  placed  perpendicular  to  the  flow  axis  If 
circumferential  ribs  were  used.  This  causes  a  problem  in  determining  the 
distribution  of  the  thrust  between  the  tensile  (or  compressive)  stresses 
in  the  liner  and  the  bending  stresses  in  the  ribs. 

For  a  model  to  use  to  approximate  these  stresses,  assume  a  contour 
as  shown  below  in  cross-section 


Section  n 


i 


Fig.  29  Sketch  of  Noxile  with  Rib*  Perpendicular  to  Flow  of  Cat 
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A  force  balance  may  be  made  on  the  i ti:  taction.  Forces  to  the  right 
are  taken  as  positive. 


O  «■  ft  4  5'»  4  -  *  -[<77  W{-  < 


■<£ 


a‘(SLiHi 


(  142  ) 


^  «  pressure  difference,  air  minus  coolant,  psi 
-4/  .  area,  Inches 
ft  •  rib  restraining  force,  pound 
w<-  m  liner  circumference  at  position  of  07  ,  inches 
£  -  liner  thickness  at  position  of  <7?  ,  Inches 


(Note  that  the  pressure  force  on  the  expansion  end  is  negative).  Since 
the  rib  is  a  circular  plate  fixed  in  the  liner  and  sleeve,  in  the  exact 
analysis  of  the  relationship  between  force  and  displacement,  one  obtains 
the  equation 


Set* 

k,  a* 


(  143  ) 


where 

S  m  rib  displacement 
£  m  rib  modulus  of  elasticity 
tr  *  rib  thickness 
a  =  outside  radius  of  rib 
b  t  Inside  radius  of  rib 

Kt  ■  some  function  of  (&  )  given  in  Table  I. 

This  same  result  may  be  closely  approximated  by  an  analytical  ex¬ 
pression  based  on  a  simple  beam  analysis.  Treating  the  rib  as  a  beam  of 
length  (a  -  b)  and  width  ZTT F 


r_  F(a-b) 

^  Zrrr  +*  £ 

(  144  ) 

e _  zrrr<f  £  4* 

(+-*) 3 

(  145  ) 

Setting  equations  (  143  )  and  (  143  )  equal  gives 


K 


(a-h)* 

Zrrr  a* 


(  146  ) 
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The  expressions  for  K|  for  severel  definitions  of  r  ere  given  below: 


a  +t> 

z 


GL-b> 


In 


t 


[af-  T^y^af 


Some  velues  of  K,  ,  as  given  by  Timoshenko  (  t>  )  and  as  calculated  from 
the  above  expressions  are  tabulated  below: 

TABLE  I 

THE  FUNCTION  K, 


r 

a 

b 

1.10 

1.25 

a 

.000120 

.00127 

b 

K, - 

.000132 

.00159 

a  +  b 

2 

K| - 

.000132 

.00142 

a  -  b 
in  a/b 
Table 

K, - 

K, - 

.000125 

.00141 

.00129 

1.50 

2.0 

3.0 

4.0 

5.0 

.0059 

.0198 

.047 

.067 

.082 

.0088 

.0396 

.  141 

.268 

00 

0 

.0071 

.0265 

.071 

.107 

.136 

.0072 

.0275 

CO 

rw. 

0 

; 

.122 

.164 

.0064 

.0237 

.062 

.092 

.114 

These  values  show  that  the  best  choice  for  r  is  the  mean  radius.  They 
also  show  that  this  choice  should  give  forces  which  are  within  />-IO%  of 

those  obtained  from  equation  (  I A3  )  for  -2 _  <  1.25.  Equation  (  145  ) 

was  used  with  7  *  g_t>  for  this  analysis. 
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Equation  (  142  )  may  now  be  wrlttan 

Pi  Qi  2/n  &t  a  f  +  r  w  +.  <r0i/_*  +*'+'  ) 

(Ai-iTr  {  *  7  (  147 ) 


-  cr  w-.;  $  Cos  (-*■■+ 

A*  $c 

(  148  ) 

/C  7T  £”  tj.. 

UT<  *  - * - 1 - “* - 

(  149  ) 

(  150  ) 

The  general  equation  for  the  force  balance  is  now 

c,  -  s.  j;  *  4  <?  -  o-, 

(  >51  ) 

In  order  to  solve  these  equations  <f  must  be  expressed  in  terms  of  (7 
or  vice  versa.  In  general  Hooke's  Law  is 


E  *  modulus  of  elasticity,  psi 
0“  -  stress,  psi 
4  «  strain,  inches 
X  ■  length,  inches 

If  d|  is  the  amount  the  i*h  segment  has  elongated  due  to  tensile 
stresses 


(  152  ) 


j.  3  (jj 

1  ZE 


(  153  ) 


The  axial  movement  is  then 


x.  ^ 

1  zc  cos  ec 


(  154  ) 
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Stresses  in  the  liner  and  rib  displacements  may  be  related  using  equation 
(  151*  ).  The  average  stress  in  a  section  is  taken  as  the  mean  of  the 
stresses  at  the  ends.  Ribs  removed  from  a  fixed  end  must  yield  the  total 
elongation  in  the  axial  direction  of  ell  sections  between  that  rib  and  the 
fixed  end,  plus  one-half  the  elongation  of  their  own  section,  thus 


/  /PTr  07.,  )  lj 
Z  (  *  /  £  Co*  0C 


Jen-/ 


If**) 


-J±±- 


e  co  s  e. 


i+t 


(  155  ) 


This  equation  holds  for  all  except  the  section  numbered  n  where 


ecese,. 


(  156  ) 


When  these  equations  are  substituted  into  equation  (  151  )  a  series  of 
(  n  -  1  )  equations  in  n  unknowns  is  obtained.  The  "extra"  unknown  is 
the  stress  at  the  exit  end, 


The  "extra"  stress,  07  ,  must  be  picked  to  satisfy  the  condition 

that  the  integrated  displacement  over  the  length  of  the  liner  is  zero, 
that  is,  the  ends  are  fixed.  This  is  done  as  follows: 


(I)  The  end  stress  07  is  arbitrarily  specified,  together  with 
the  other  stresses,  and  equation  (  1 5 1  )  solved  by  relaxation  of 
residuals.  This  gives  a  value  of  OI  ,  and  <f  is  calculated  from 
equation  (  155  ). 


(2)  The  fixed  exit  end  condition  requires  that 


157  ) 


Letting  the  calculated  value  of  the  left  side  of  this  equation  be 


Jt, 

€  ems  es>, 


(  158  ) 


(3)  If  V  >0  the  liner  extends  too  far  in  the  direction  of 
flow,  and  to  shorten  it  CT;  must  be  more  compressive  so  O';  must  be 
reduced.  If  ¥<0  ,  must  be  increased.  X  may  be  made  as  close 

to  O  as  desired  by  changing  Oj  in  the  proper  direction  and  repeating 
the  calculation  in  (I). 
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Calculations  were  made  for  a  representative  liner  configuration  with 
a  Mo-TI  liner  and  ^  200  atm.  entrance  air  pressure.  The  liner  had 
20  ribs  .025  inches  thick  and  spaced  .25  inches  apart  along  the  axis. 
Results  are  given  in  Table  II.  The  computer  program  used  for  the  calcu¬ 
lations  was  written  in  Fortran  II  language  and  is  Included  in  Appendix  A. 
The  calculations  were  made  on  an  IBM  1620  Digital  Computer. 

(f)  If  the  liner  in  this  case  Is  fixed  at  the  ends,  compressive 
longitudinal  stresses  will  be  developed.  For  restraint  by  the  sleeve  only 
these  will  be  the  same  as  In  Case  I,  Part  (f) ,  thus 


O'tu  -  - 


E, 


and 


6*4  f:) 

tTm  £,  gf 

(g)  The  principle  tangential  stresses  for  this  case  are: 
(I)  At  the  inside  radius  over  a  rib 


rn. 


(  159  ) 

(  160  ) 


^cu. 


CT" 

'“'NHl 


(  161  ) 


(2)  At  the  outside  radius  over  a  rib 


^Pciii  ”  ^siti  ~  ^rw  (  162  ) 

(3)  At  the  inside  radius  between  ribs 

^peili  =  (  163  ) 

(4)  At  the  outside  radius  between  ribs 

^Pcizt  m  +•  (  164  ) 
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TABLE  II 

RESULTS  FOR  A  REPRESENTATIVE  LINER  CONFIGURATION 


Section 

Number 

Assumed  A I r 
Pressure  (psi) 

Dlemeter 

(inches) 

Liner  Thickness 
(Inches) 

Tensile  Stress 
(psi) 

1 

2800 

1*32 

.034 

200 

2 

2800 

1.24 

.033 

510 

3 

2800 

1.16 

.032 

695 

4 

2800 

1.08 

.031 

895 

5 

2300 

1.00 

.030 

1,220 

6 

1400 

*92 

.029 

990 

7 

1300 

.84 

.028 

625 

8 

1200 

76 

.028 

45 

9 

1100 

.68 

.027 

O 

O 

CO 

t 

10 

1000 

.60 

.026 

-  440 

II 

900 

*52 

.025 

1 

Xr 

O 

O 

12 

800 

.44 

.024 

-  2,540 

13 

700 

•  36 

.023 

-  4,700 

14 

600 

.28 

.022 

-  8,220 

15 

500 

.  24 

.025 

-  7,940 

16 

400 

.60 

.028 

-  4,100 

17 

300 

1.04 

*031 

-  2,100 

13 

200 

1.41 

.033 

10,475 

19 

100 

1*50 

.034 

8,185 

20 

5 

1*50 

035 

7,930 
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(h)  The  principle  longitudinal  stresses  for  this  case  ere: 

(I)  At  the  inside  radius  over  a  rib 


cp  +  a—  _  cr  *  err 

TtteOiT  TU  TU  r  IS 


•I  l 


(  165  ) 


(2)  At  the  outside  radius  over  a  rib 


°TU  +  ^isiil 


(  166  ) 


(3)  At  the  inside  radius  between  ribs 


Gputl  *  ^T»*U»T  +  ^Tt|  “  °th  ♦  ^|,v 


(  167  ) 


(k)  At  the  outside  radius  between  ribs 


^TKWJT  *  ^n.1  ®TM  *■  *^SIXX. 


(  168  ) 


Note:  The  ribs  cause  a  small  additional  bending  stress  which  is  neglected 
here. 

(i)  The  radial  principle  stresses  are  given  in  (d). 


Cos*  3  -  Fully  Retrained  Liner 

This  is  the  case  of  real  interest  in  the  present  analysis  as  sample 
calculations  of  Case  2  show  that  unless  the  sleeve  is  made  thinner  than 
the  liner  little  is  gained  stresswise. 

This  case  is  calculated  from  the  equations  developed  in  Case  2  (Longi¬ 
tudinal  Ribs)  by  letting  the  outside  radius  of  the  sleeve  be  much  larger 
than  the  inside  radius,  thus,  using  a  very  thick  sleeve.  In  Case  2  - 
(Circumferential  Ribs)  this  case  is  calculated  by  assuming  a  very  large 
rib  cross-sectional  area. 

Since  a  large  number  of  variables  were  involved  these  equations  were 
programmed  in  Fortran  il  language  for  the  IBM  1620  Digital  Computer.  This 
program  is  attached  and  described  in  Appendix  A. 
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CONCLUSIONS 

Parametric  stress  studies  and  coolant  flow  considerations  show  that 
longitudinal  ribs  offer  the  best  possibilities  for  backside  cooling.  If 
circumferential  ribs  are  placed  perpendicular  to  the  flow  axis  high  bend¬ 
ing  stresses  are  introduced  in  the  ribs  and  liner  due  to  the  thrust  load. 

The  problem  of  getting  proper  coolant  flow  at  each  section  is  also  compli¬ 
cated  and  causes  design  and  fabrication  difficulties. 

The  maximum  design  pressure  is  500  atmospheres  in  the  entrance  tube. 
This  would  produce  a  pressure  of  approximately  250  atmospheres  at  the 
throat.  Calculations  were  made  using  these  pressures  and  assuming  a  heat 
load  of  60  B/in‘sec  for  a  Ho-TI-Zr  liner.  A  backside  cooling  coefficient 
of  150,000  B/hr  ft*  °f  was  assumed.  The  bulk  coolant  temperature  was 
assumed  to  be  85*F.  A  heavy  sleeve  was  used  and  the  liner  thickness  varied. 
The  results  of  these  calculations  are  presented  graphically  on  pages  75 
through  78>  Some  properties  of  Ho-Ti-Zr  are  shown  graphically  or.  page  79. 

Several  observations  may  be  made  from  these  graphs:  the  critical 
octahedral  shear  stress  for  yielding,  feer-  ,  is  exceeded  for  all  liner 
thicknesses  and  at  all  four  locations  for  the  1.5  inch  diameter  entrance 
sections;  at  two  positions  ?Ur  Is  exceeded  for  the  0.78  inch  diameter 
throat;  the  octahedral  shear  stresses  for  the  0.27  inch  diameter  throat 
may  be  made  •£  ^critical  bT  making  the  liner  thickness  6.  10  mils;  the 
octahedral  shear  stresses  for  the  0.042  inch  throat  arei.  T'crjt|caj  for 
a  liner  thickness  4  14  mils.  It  should  be  noted  that  although  the  curves 
are  not  plotted  below  10  mils  liner  thi ckness,  al I  the  octahedral  shears 
approach  infinity  as  the  liner  thickness  approaches  0.  An  optimum  thick¬ 
ness  for  the  liner  may  be  picked  from  the  curves  to  minimise  the  plastic 
strain,  and  for  this  thickness  the  mean  longitudinal  compressive  stress, 
which  is  the  maximum  mean  stress,  may  be  found.  These  values  are: 


Liner  Diameter 
(Inches) 

1.50 

.784 

.272 

.042 


Liner  Thickness 
(Inches) 

.034 

.018 

.012 

.006 


Rib  Spacing 
(Inches) 

.226 

.120 

.059 

.014 


Nean  Longitudinal 
Stress  (psi) 

-65,000 

-52,000 

-47,000 

-40,000 


The  mean  longitudinal  stress  is  compressive,  and  does  not  exceed  the  ulti¬ 
mate  strength  of  the  liner  material  for  any  of  the  conditions. 

A  brief  calculation  may  be  made  to  predict  the  operating  cycles  before 
failure  for  these  nozzles.  The  maximum  principle  stress  occurs  midway 
between  ribs  at  the  inside  of  the  liner  and  at  the  1.5  inch  diameter  section. 
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Longitudinal  Ribs  2 

Q  «  60  Btu/in,  sec 
Material  -  Mo-Ti~Zr 
Diameter  -  1.5  in. 

Pressure  -  500  atm  2 

Backside  Coefficient  -  150,000  Btu/hr  ft  F 

S1M1  -  Octahedral  Shear  Inside  Liner  over  a  Rib 
S1M2  -  Octahedral  Shear  Outside  Liner  over  a  Rib 
S2M1  -  Octahedral  Shear  Inside  Liner  between  Ribs 
S2M2  -  Octahedral  Shear  Outside  Liner  between  Ribs 


Liner  Thickness,  mils 


Rib  Spacing  for  Indicated  Liner  Thickness,  mils 
fi9.  30  Strait  vanvi  Linar  Thicknan;  Diomatar  »  1.5" 
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Longitudinal  Ribs 

Q  -  60  Btu/in.^  sec 
Material  -  Mo-Ti-Zr 
Diameter  -  0.78  in. 

Pressure  -  250  atm  9 

Backside  Coefficient  -  150,000  Btu/hr  ft  °F 


89  128  159  185 

Rib  Spacing  for  Indicated  Liner  Thickness,  mils 

Fig.  31  Stress  Versus  Liner  Thickness;  Diameter  «  0.78" 
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9 


Longitudinal  Ribs 

Q  *  60  Btu/in.‘ 
Material  -  Mo-Ti-Zr 
Diameter  -  0.27  in. 
Pressure  -  250  atm 


sec 


2ot 


120 


100 


Backside  Coefficient  -  150,000  Btu/hr  ft  F 


n 

i 

o 


m 

a 


m 

x 

« 

u 

GO 


80 


60 


40 


20 


J. 


10  20  30  40  50 

Liner  Thickness,  mils 


60 


54  78  99  118 

Rib  Spacing  for  Indicated  Liner  Thickness,  mils 


Fig.  32  Stress  versus  Liner  Thickness;  Diameter  -  0.27" 
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Longitudinal  Ribs 

Q  -  60  Btu/in.2  sec 
Material  -  Mo-Ti-Zr 
Diameter  «  0.042  in. 
Pressure  «=  250  atm 


24  40  55  69 

Rib  Spacing  for  Indicated  Liner  Thickness,  rails 


Fi*.  33  St,...  v.r.u.  Lin.,  Thiefcn.,,;  Diom.t.r  -  0.042  ■' 
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Fig.  34  Strength  versus  Temperature  for  Mo-Ti-Zr 
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The  computed  value  of  this  stress  is  -232,000  psi ,  which  exceeds  the  yield 
strength  of  the  materiel  by  172,000  psi  when  the  inside  temperature  is 
I470*F.  The  modulus  of  elasticity  at  this  temperature  is  35,000,000  psi. 
The  plastic  strain  is  thus 


4  ~  172,000  .  0.00492  In/In 

p  35,000,000 

Using  an  average  value  of  0.36  for  the  constant  In  Coffin's  (  6  )  equation 

"  ‘[-Mr]  2 

and  N  —  5000  cycles  to  failure. 

Though  not  presented  graphically  in  this  report,  calculations  for 
Ho-Ti-Zr  liners  at  lower  pressures  and  heat  loads  indicate  that  these  same 
liner  thickness  and  rib  spacing  values  are  acceptable.  Calculations  may  be 
rapidly  made  for  any  assumed  conditions  using  the  Nozzle  Stress  Program. 

In  the  actual  nozzle  design  it  ts  proposed  that  the  ribs  will  terminate 
just  before  reaching  the  throat  as  the  pressure  drops  sharply  at  this  point. 
This  simplifies  the  nozzle  construction  and  reduces  the  coolant  pressure  loss 
in  the  longer  expansion  section,  it  also  alleviates  the  problem  of  close 
rib  spacing  at  the  throat  of  the  ,042  inch  nozzle.  The  ribs  may  be  termi¬ 
nated  when  the  liner  diameter  is  .3  inch  and  a  rib  spacing  of  ~  . 1  inch. 

A  sketch  of  a  typical  nozzle  design  is  shown  on  page  81.  In  this  sketch 
the  ribs  are  shown  attached  to  the  sleeve  rather  than  the  liner.  This  offers 
some  advantage  in  so  far  as  construction  and  longitudinal  compressive  stresses 
are  concerned.  However,  there  may  be  some  distortion  at  the  throat  and  some 
heat  transfer  problems  on  the  backside.  It  is,  therefore,  recommended  that 
the  ribs  be  machined  integral  with  the  liner  if  costs  are  not  prohibitive. 


VI.  MATERIALS  STUDY 

The  selection  of  the  best  material  to  use  for  a  nozzle  liner  for  high 
pressure-high  temperature  air  must  be  based  on  the  following  considerations: 

(a)  Ability  to  stand  the  imposed  stresses.  A  stress  parameter  will 
include  as  factors: 

(1)  Thermal  conductivity 

(2)  Strength  at  working  temperatures 

(3)  Modulus  of  elasticity 

(4)  Thermal  expansion  coefficient. 
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(b)  Resistance  to  oxidation  at  elevated  temperatures.  This  factor 
includes: 

(1)  Oxidation  resistance  of  the  material  chosen. 

(2)  Possibility  of  protection  by  coatings. 

(c)  Fabricabi 1 1 ty,  which  Includes: 

(1)  Mach i nab i 1 1 ty. 

(2)  Joining  methods. 


STRESS  PARAMETER 

it  is  noted  that  the  thermal  stress  in  a  structural  member  across 
which  there  is  a  temperature  difference  is  directly  proportional  to  the 
temperature  difference,  the  thermal  expansion  coefficient,  and  the  modulus 
of  elasticity.  For  a  given  heat  load  the  temperature  difference  is  inverse 
ly  proportional  to  the  thermal  conductivity.  A  thermal  stress  parameter 
should  also  include  the  ultimate  strength.  Thus: 


(  169  ) 


%  *  stress  parameter,  B/hr  ft 
K  =  thermal  conductivity,  B 

hr  ft  °F 

£  3  modulus  of  elasticity,  lb/in2 
°<  =  thermal  expansion  coefficient,  in/in  ®F 
Su  *  ultimate  strength,  lb/inz 

Figure  36,  page  33,  shows  this  parameter  plotted  vs  temperature  (*F) 
for  the  following  materials:  Cu-Be  alloy,  Cu-Zr  alloy,  Mo-0.5  Ti-0.07  Zr, 
AM  350  steel,  modified  H - 1 1  steel.  It  is  evident  from  the  curves  shown 
that  above  approximately  900'F  the  Mo-TI-Zr  alloy  is  a  much  better  material 
to  use  so  far  as  the  thermal  stresses  are  concerned  and,  at  the  high  heat 
loads  expected  in  the  nozzle,  these  stresses  are  dominant. 


OXIDATION  RESISTANCE 

The  oxidation  rate  of  pure  copper  is  given  by: 
w‘>  kt 

where  w  *  weight  of  oxide  formed  in  grams/cm^ 
■i  =  time  in  seconds 
yk  =  function  of  temperature  in  °F 
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values  of  k  are  {_]_): 

k  x  1010 

Temp.  #F 

Pure  O2 

Air 

750 

.044 

- 

950 

.440 

- 

1 100 

3.24 

- 

1300 

16.0 

8.03 

Assuming  that  the  rate  In  air  at  the  lower  temperatures  is  half  that  In 
pure  O2  the  penetration  in  cm/hr  Is: 


5 


(it) 


4 


%  m  cm/hr  of  oxidation 
Jl  ■  function  of  temperature 
t  -  time  in  seconds 
e/  =  density  in  grams/cm’ 

p\  _  'Q~/0X  3  *  *  ><>') J 

J,fof  ~  S.<56 

P.)  ~  /O  Cm/hr. 

'ifO*  ' 


Even  if  the  oxide  is  continuously  removed  the  rate  of  loss  is  very  small 
at  IIOO'F,  and  the  inside  temperature  of  the  copper  alloys  would  not  exceed 
this  temperature  except  at  extremely  high  heat  loads  (  >  200  B/in^  sec). 
Oxidation  of  the  copper  alloys  does  not  appear  to  be  a  large  problem  if 
they  are  used. 

The  Mo-Ti  and  Mo-Ti-Zr  alloys  are  not  as  resistant  to  oxidation  as 
the  Cu  alloys.  The  oxide  of  Mo  is  also  volatile,  so  there  is  no  self  pro-' 
tection.  There  are  metallic  coatings  available,  however,  such  as  electro¬ 
plated  or  hot  dipped  Ni,  which  will  protect  Mo  in  the  I300-1500°F  range. 
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For  protection  at  much  higher  temperatures,  say  2000-2500°F,  a  slllclde 
coating  known  as  Chromalloy  W-3  Is  available.  It  Is  thus  possible  to 
provide  oxidation  protection  for  Ho  and  Mo  alloys  for  several  hours  at 
the  temperatures  of  interest. 


FABRICABILITY 

Fabricability  is  an  Important  consideration,  as  the  optimum  thickness 
of  the  liner  at  some  positions  is  10-15  mils.  The  thrust  load  requirements 
must  be  met  by  using  ribs  on  either  the  liner  or  the  sleeve.  If  the  ribs 
are  a  part  of  the  liner  they  must  be  machined  on  this  surface  and  thus,  the 
machlnabi 1 i ty  of  the  material  is  very  important. 

The  copper  alloys,  Including  Cu-Be  and  Cu-Zr,  have  a  machlnabi I i ty 
index  of  20  (based  on  free-cutting  brass  =  100).  Cu-Be  and  Cu-Zr  may  both 
be  joined  by  soldering  and  brazing  (  8  ) . 

Molybdenum  and  its  alloys  Mo-Ti  and  Mo-Ti-Zr  machine  much  like  cast 
iron  (  8  ) .  Tungsten  carbide  tools  should  be  used  (  8  ) .  The  machinabili- 
ty  index  is  approximately  the  same  as  that  of  the  copper  alloys. 


CONCLUSION 

The  following  observations  are  made: 

(a)  If  the  backside  cooling  technique  (forced  convection  or  boiling), 
the  heat  load,  and  the  liner  thickness  are  such  that  the  temperatures  at 
the  inside  of  a  Cu-Be  liner  remain  below  300°F,  Cu-Be  is  the  preferred 
liner  material.  Under  these  conditions  it  has  the  largest  stress  parameter 
and  does  not  require  oxidation  coating. 

(b)  If  the  backside  cooling  technique  (forced  convection  or  boiling), 
the  heat  load,  and  the  liner  thickness  arc  such  that  temperatures  at  the 
inside  of  a  Cu-Zr  liner  are  between  800°F  and  1000°F,  and  if  It  is  desired 
to  avoid  the  added  expense,  complication,  and  uncertainty  of  coating  the 
liner,  Cu-Zr  is  attractive. 

(d)  If  the  heat  loads  and  backside  cooling  technique  (forced  convection 
or  boiling)  are  such  that  the  conditions  in  (I)  and  (2)  are  not  met,  Mo-Ti 
or  Mo-Ti-Zr  are  plainly  superior. 

As  the  heat  load  is  increased  for  any  given  backside  cooling  tech¬ 
nique,  the  inside  liner  temperature  rises.  Studies  of  the  heat  load  in 
the  upper  region  of  the  envelope  show  that  with  radiation  included  this 
load  may  easily  be  several  hundred  B/inZsec.  This  would  produce  inside 
surface  temperatures  of  >  300°F  using  Cu-Be.  Since  Mo-Ti  and  Mo-Ti-Zr 
have  higher  stress  parameters  than  Cu-Zr  and  this  advantage  increases  as 
the  heat  load  is  increased,  it  is  recommended  that  Mo-Ti-Zr  be  used  in 
order  to  push  as  far  as  possible  in  the  envelope.  At  lower  heat  loads 
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and  Inside  temperatures,  metallic  coatings  should  be  used.  At  higher 
heat  loads  and  temperatures  (  ~  2000° F),  Chroma  1 1 oy  W-3  coating  should 
be  used. 


VII.  SHARP  EDGED  NOZZLE 

Since  the  heat  transfer  rate  across  a  laminar  boundary  layer  is 
small  compared  to  the  heat  transfer  rate  across  a  turbulent  boundary 
layer,  it  is  desirable  to  design  a  nozzle  in  which  a  laminar  boundary 
layer  exists  as  far  downstream  as  possible.  For  this  reason,  the  sharp 
edged  nozzle  shown  in  Figure  37,  is  proposed  for  an  analsis.  The  primary 
advantage  of  this  type  of  nozzle  is  the  short  distance  of  the  upstream 
air-metal  contact  relative  to  a  conventional  nozzle  approach  section. 

From  the  above  arguments,  it  is  concluded  that  it  will  be  desirable 
to  have  the  radius  of  curvature  of  the  nozzle  as  small  as  possible;  there¬ 
fore,  it  becomes  necessary  to  attain  the  shape  of  the  sonic  line  for  any 
proposed  nozzle  of  this  type  since  this  information  is  required  for  the 
design  of  the  downstream  section  of  the  nozzle.  Other  questions  which 
need  to  be  answered  about  the  proposed  nozzle  are:  (I)  how  much  water  is 
carried  downstream  by  the  air,  (2)  does  the  water  jet  transport  air  out 
of  the  main  water  stream,  and,  if  so,  how  much,  and  (3)  what  influence  will 
the  eddies  have  on  the  design  of  the  nozzle?  The  analysis  which  has  been 
accomplished  to  date  is  briefly  outlined  in  the  following  sections. 

The  first  region  of  analysis  was  confined  to  the  water  air  Interface, 
see  Figure  37.  Since  the  water  emerges  from  the  jet  at  a  high  velocity 
(  ''-'125  ft  per  sec),  it  was  assumed  that  the  water  velocity  was  large 
relative  to  the  adjacent  air  velocity.  This  assumption  implies  that  the 
emergent  water  can  be  treated  as  a  wall  jet  flowing  into  a  semi - i nf ini te 
medium.  The  Reynolds  number  of  the  jet  based  on  the  diameter  of  the  inlet 
annulus  is  about  6  x  10  ;  therefore,  the  problem  is  that  of  a  turbulent 
wall  jet  (with  the  implicit  assumption  ot  equivalent  static  pressures  in 
the  air  and  water  streams).  Some  analytical  work  and  experimental  results 
were  found  for  wall  jets  for  the  case  of  a  fluid  emergent  into  a  large 
medium  which  is  filled  with  the  same  fluid,  but  no  information  was  found 
for  the  case  of  dissimilar  fluids.  Since  the  static  pressure  in  the  two 
streams  is  assumed  to  be  equivalent,  it  was  concluded  that  the  only  driv¬ 
ing  force  to  cause  mixing  would  be  the  momentum  flux,  and  hence,  the  only 
physical  properties  of  significance  would  be  the  density  and  viscosity. 

On  this  basis,  it  was  decided  to  treat  the  jet  as  water  emergent  into  water 
and  calculate  a  Reynolds  number  of  the  fluid  crossing  into  the  jet.  This 
Reynolds  number  is  based  on  the  component  of  velocity  normal  to  the  jet 
axis.  Assume  that  the  air  "sucked1  into  the  water  jet  would  be  of  such 
magnitude  that  it  would  have  an  equivalent  Reynolds  number  to  that  cal¬ 
culated  above. 

M.  B.  Glauert  (_SL)  developed  an  analytical  solution  for  the  wall  jet 
assuming  a  similarity  solution.  His  results  show  that  the  characteristic 
jet  width  is  essentially  proportional  to  the  distance  along  the  axis 
parallel  to  the  jet  wall  raised  to  the  first  power.  An  experimental  in- 
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vestlgation  of  a  wall  jet  was  carried  out  by  P.  Bakke  (  10  )  and  Inde¬ 
pendently  by  W.  H,  Schwarz  and  W.  P.  Cosart  (II). 

From  the  above  references  It  was  concluded  that  the  jet  would  spread 
at  about  3-5°;  therefore,  the  dimension  "a"  (see  Figure  37)  is  known. 

Since  ''a"  is  chosen  to  Include  the  water  "spread'1,  there  will  be  no  water 
carried  into  the  air  stream  due  to  mixing;  however,  the  model  chosen  re¬ 
quires  entrainment  of  the  surrounding  fluid  into  the  jet.  It  should  be 
observed  that  the  model  under  analysis  neglects  diffusion  of  the  water  into 
the  air  stream;  however,  it  is  felt  that  the  water  transported  into  the  air 
stream  by  this  mechanism  will  be  small  because  of  the  short  distance  of  water 
travel  in  the  nozzle. 

Referring  to  Figure  37,  it  Is  seen  that  the  length  of  the  water  jet 
is  less  than  three  equivalent  diameters;  therefore,  the  velocity  profiles 
will  not  be  susceptible  to  a  similarity  type  solution  since  the  similarity 
solution  is  valid  only  at  far  downstream  conditions.  Since  the  entrance 
stream  is  not  subject  to  similarity  analysis,  the  conclusions  drawn  might 
be  questioned;  however,  direct  observation  of  an  air  jet  under  analysis  at 
The  University  of  Tennessee  shows  that  the  air  wall  jet  pulls  in  air  (smoke) 
from  the  surrounding  medium  in  the  entrance  region  as  well  as  far  down 
stream.  On  this  basis,  it  is  felt  that  the  conclusions  drawn  from  a  similar! 
ty  analysis  are  quantitatively  correct.  No  attempt  will  be  made  to  use  the 
similarity  velocity  profile  to  attain  the  quantity  of  air  pulled  into  the 
water  jet. 

If  eddies  are  developed  In  the  region  near  the  throat  it  is  possible 
that  they  will  change  the  effective  throat  geometry;  therefore,  it  would 
be  difficult,  if  not  impossible,  to  calculate  the  shape  of  the  sonic  line 
(hence  design  of  the  downstream  section  of  the  nozzle).  Since  air  is  being 
pulled  into  the  water  jet  the  eddies  will  be  depressed,  and  conceivably  can 
be  controlled  to  the  point  where  there  are  essentially  no  eddies  formed  in 
the  section  near  the  throat.  On  this  basis,  it  was  concluded  that  the  for¬ 
mation  of  eddies  would  be  of  no  consequence  in  the  design  of  the  nozzle. 

Before  considering  the  preceding  problems  in  a  more  rigorous  fashion, 
it  was  decided  to  examine  the  shape  of  the  sonic  line  at  the  throat  for 
different  throat  diameters  and  wall  geometries,  and  then  determine  what 
flow  regime  exists  in  the  boundary  layer. 

The  determination  of  the  sonic  line  is  not  complete,  but  the  equations 
to  be  solved  have  been  developed.  The  development  of  the  equations  to  be 
used  for  determining  the  shape  of  the  sonic  line  is  as  follows. 

The  solution  to  the  potential  flow  problem  in  a  region  near  the  throat 
outlined  here  follows  the  procedure  described  by  Oswatitsch  and  Rothstein 

(_LL>- 


Use  the  nomenclature  as  presented  in  Figure  38. 
t inui ty  is  wri tten  as 


The  equation  of  con- 


(  170  ) 
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Fig.  37  Sketch  of  Sharp  Edged  Nozilo  (Not  to  Scale) 
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where  u  end  v  ere  the  velocity  components  In  x  end  y  directions, 
end 


cO 


a 


(  171  ) 


Solve  equetion  (  170  )  for  v  by  making  en  approximation  forfiO.  Take 
the  value  of  fio  to  be  the  stream  filament  value  designated  by  pti/s 
The  equation  of  energy  can  be  written  as 


J.  * 

2 


y 

r-/ 


JL  _£* 


(  172  ) 


where  V  Is  the  ratio  of  .pecific  heats  and  the  superscript  0  refers  to 


stagnation  conditions.  In  equation  ( 
at  this  condition,  then 

c\  jLSL  jl* 

V-J 

Let 

a  J&L 

c. 


72  )  let  T-^O  ,  and  define  u)*c 

(  173  ) 

(  174  ) 


and 


(  175  ) 


Since  the  main  flow  Is  assumed  to  be  non-vtscous,  the  flow  Is  isentropic, 
I  .a. , 


(  176  ) 


Substitute  equations  (  173  )  through  (  176  )  Into  equation  (  172  ),  then 
equation  (  172  )  becomes 

X  Y'l 

^  (f)  (  177  > 
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Since  equation  (  177  )  is  in  a  simple  form  the  quantities 
be  used  to  make  all  velocities  and  density  non-dimensional 
consequently,  the  primes  will  be  dropped. 

The  equation  of  continuity  for  the  stream  filament  is 


a/*  / 


where  f  is  as  indicated  in  Figure  38. 

Differentiation  of  equation  (  l 78  )  yields 

Equation  (  1 70  )  becomes  (using  equation  (  179  )  ) 


4- (rv 'i)  -  4 

au  r 


dy 

Since  there  is  symmetry  about  the  centerline,  v  =  o  at 
Integration  of  equation  (  180  )  (with  o *  /  )  yields 

Assume  u  to  be  given  by, 

V( \y)  *  +  -7  y*-* . 

Potential  flow  ( i rrotational  flow)  requires  that 

o)c/  -  c)  V 
dX 

Using  equation  (  1 8 1  ),  the  first  two  terms  of  equation  ( 
carrying  out  the  operations  indicated  in  equation  (  1 83  )» 


U  -  Om 


i *  i  (v,  ff  -  )  4 


c  and  p*  will 
henceforth ; 

(  178  ) 

(  179  ) 

(  180  ) 

y  ■  o. 

(  181  ) 

(  182  ) 

(  183  ) 

182  ),  and 
the  result  is 

(  184  ) 
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The  parameter  In  equation  (  184  )  Is  attained  from  the  equation  of 

continul ty 


M*  21 Tj 


rf  Up  if  Jy 

a 

From  equation  (  171  )  and  equation  (  177  ) 

U-  > 

,  V  \  i'"/  L  11  F~! 

Aa.  / /-  UJ  J  m.  (  /-  U  -  V  J 


(  185  ) 


(  186  ) 


Def I ne 


/ 

*-/ 


G>(uyK)  s  u(/-u%-v*‘) 

then  equation  (  185  )  reads  as 

f 

A?  *  2/rf  0(a /*)  y  e/tj 


Solution  of  equation  (  188  )  yields 

-  ±(“s  ff"+  tf') 

ig  equi 
)  ylel 


(  187  ) 


(  188  ) 


(  189  ) 


Solving  equation  (  171  )  for  uQ  and  substituting  the  result  into  equation 
(  184  )  yields 


^  '// ') 

*  i  (“,  f{ "*  ffk -  *4/') £ 

For  the  nozzle  presented  in  Figure  38,  f  is  found  to  be 

'■i 


(  190  ) 


(  191  ) 
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Equations  (  190  )  and  (  1 9 1  )  have  not  been  solved;  however,  the  solution 
of  these  equations  should  yield  the  flow  field  In  the  region  near  the 
throat,  and  in  particular,  the  shape  of  the  sonic  line  with  the  radius 
of  curvature  of  the  nozzle  as  a  parameter.  D*  and  us  will  be  specified 
by  the  design  criteria. 


VIII.  BACKSIDE  COOLING  STUDY 

The  heat  loads  on  the  liner  of  any  proposed  nozzle  are  calculated  to 
be  very  high.  The  heat  load  postulated  for  the  stress  analysis  is 
60  B/ln^sec  or  30  x  10b  B/hr  ft*.  |f  only  backside  cooling  without  film 
or  transpiration  assistance  is  used  the  loads  In  the  upper  portion  of  the 
envelope  may  exceed  this  value  several  times.  This  not  only  causes  problems 
in  liner  materials,  temperatures,  and  stresses,  but  also  in  problems  with 
the  removal  of  the  heat  at  the  backside  of  the  liner.  Unless  a  high  heat 
transfer  coefficient  on  the  coolant  side  can  be  attained,  an  excessive 
rise  in  the  liner  temperature  will  occur. 

It  is  proposed  that  water  be  used  as  the  backside  coolant  since  it  is 
readily  available  in  sufficient  quantities  at  high  pressures.  The  thermo¬ 
dynamic  and  heat  transfer  properties  are  known  and  several  equations  are 
available  which  give  a  good  prediction  of  the  heat  transfer  coefficient  at 
various  flow  conditions.  Water  boils  at  fairly  low  temperatures  (at  low 
pressures)  offering  the  possibility  of  a  high  heat  transfer  coefficient 
with  a  low  flow  rate  and  pressure  drop. 

Two  heat  transfer  modes  are  of  interest  using  water  as  the  coolant, 
namely,  forced  convection  and  boiling.  Forced  convection  exists  at  all 
positions  along  the  coolant  channel  when  the  static  pressure  is  high  enough 
to  cause  the  saturation  temperature  to  equal  or  exceed  the  liner  coolant 
surface  temperature.  Boiling  heat  transfer  exists  when  the  liner  surface 
temperature  is  allowed  to  exceed  the  saturation  temperature  of  the  coolant 
by  an  amount  sufficient  to  cause  vaporization  of  the  coolant  at  spots  on 
the  surface.  This  is  called  nucleat'  boiling,  and  in  conjunction  with  the 
forced  fluid  flow,  will  be  referred  to  as  surface  boiling. 


FORCED  CONVECTION 

In  order  to  design  the  coolant  annulus  for  forced  convection  cooling, 
it  is  necessary  to  calculate  the  static  pressure,  bulk  temperature,  surface 
temperature,  and  heat  transfer  coefficient  at  each  position  along  the  coolant 
channel.  The  inlet  pressure,  flow  rate,  and  coolant  channel  geometry  must 
then  be  adjusted  to  assure  that  the  heat  will  be  transferred  without  caus¬ 
ing  the  surface  temperature  to  exceed  the  saturation  temperature  of  the 
coolant  at  any  location.  Since  the  properties  of  the  coolant  such  as  vis¬ 
cosity,  Prandtl's  number,  and  density  are  functions  of  the  bulk  temperature 
and  these  vary  with  the  flow  rate  and  the  heat  load  distribution  along  the 
channel,  a  stepwise  fluid  flow  and  heat  transfer  calculation  along  the 
channel  is  required. 
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The  static  pressure  distribution  Is  calculated  from  the  general 
equation 


T>  -  -  APj 


(  192  ) 


whe  re  "Rl  ■ 

*R  - 

a?  - 

» 


static  pressure  at  exit  of  a  section,  psi 
static  pressure  at  entrance  to  a  section,  psi 
loss  of  pressure  due  to  frictional  drag,  psi 
loss  of  pressure  due  to  acceleration,  psi. 


The  frictional  pressure  loss  Is  calculated  from  the  Blasius  formula 
for  fully  developed  turbulent  flow  with  fluid  properties  evaluated  at  the 
bulk  temperature, 


(  193  ) 


where  Ax  • 
Dm  * 

<°  - 
V  - 

9b  * 

R,  . 

M.  - 


length  of  channel  section,  inches 

hydraulic  diameter  of  channel,  Inches 

density  of  coolant,  Ibm/ft^ 

coolant  mean  velocity,  ft/sec 

32.2  Ibm  ft/lbf  sec2 

Reynolds  number 

coolant  viscosity,  Ibm/ft  sec. 


The  dynamic  pressure  loss,  APj  ,  15  calculated  from 

APj  -  ftVm.  ~  Q  V. 

2  9* 


(  19^  ) 


The  pressure  loss  due  to  the  sharp  bend  at  the  coolant  channel  entrance 
is  calculated  from  the  Blasius  formula,  equation  (  193  ),  using  an  equiva¬ 
lent  section  length  in  place  of  Ax  .  The  equivalent  section  length  used 
in  the  present  calculations  Is—  4S  DH  ,  thus, 


(  195  ) 


The 
equation 
the  film 


surface  heat  transfer  coefficient  is  calculated  using  the  Colburn 
modified  by  evaluating  Reynolds  number  and  Prandtl's  number  at 
temperatures 

h.  *  O.OZ3  C-  Cr  R. 

C  P  \ 


0.2 

.  v 


(  »96  ) 
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where  - 


C*  - 
G  ■ 


surfece  heat  transfer  coefficient,  B/hr  ft**F 
specific  heat  of  coolant,  B/lbm  *F 
mass  flow  rate,  Ibm/hr  ft* 

Reynolds  number  (film) 

Prandtl's  number  (film). 


The  bulk  temperature  along  the  coolant  channel  is  calculated  from 
the  heat  balance  equation  on  each  section, 


Tb 

2 


Tb 


QAa 

G  Ac  Cf, 


(  197  ) 


whe  re  "7^  * 


bulk  temperature  at  section  exit,  *F 
bulk  temperature  at  section  entrance,  *F 
liner  heat  load,  B/hr  ft* 
liner  section  surface  area,  ft? 
coolant  mass  flow  rate,  Ibm/hr  ft* 
coolant  channel  cross-sectional  area,  ft* 
coolant  specific  heat,  B/lbm  °F. 


The  calculational  procedure  is  to  choose  a  nozzle  geometry,  coolant 
flow  rate,  and  entrance  pressure  for  a  given  heat  load  distribution.  The 
bulk  temperature  and  static  pressure  distributions  are  first  calculated 
using  equations  (  192  ),  (  193  ),  (  >9^  ),  (  195  ),  and  (  197  ).  The  heat 
transfer  coefficients  are  calculated  using  equation  (  196  )  by  successive 
approximation  from  the  surface  temperature  calculation  where 


7;  *  3 


(  193  ) 


wi  th 


"Tj  =  surface  temperature,  °F 
7J,  =  bulk  temperature,  °F 
G)  =  liner  heat  load,  B/hr  ft* 

\  =  heat  transfer  coefficient  based  on  film 
temperature,  B/hr  ft*  °F. 


If  the  surface  temperature  exceeds  the  saturation  temperature  by 
more  than  10  °F  at  ant  section,  adjustments  are  made  in  entrance  pressure, 
channel  geometry,  and  flow  rate  and  the  calculation  repeated  until  the 
desired  conditions  are  achieved. 

A  general  computer  program  was  written  for  the  calculations.  It 
includes  the  boiling  heat  transfer  calculations  described  below.  The  com¬ 
puter  program,  entitled  Backside  Cooling  Program,  is  discussed  in  Appendix 
B. 
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BOILING  HEAT  TRANSFER 

The  problem  with  the  backside  cooling  Is  in  obtaining  high  heat  trans¬ 
fer  coefficients  in  order  to  hold  down  the  liner  coolant  side  surface 
temperature.  With  forced  convection  cooling  extremely  high  velocities  are 
required  to  achieve  heat  transfer  coefficients  from  75,000  to  100,000 
B/hr  ft*  *F.  This  causes  high  entrance  pressures,  large  frictional  pressure 
loss  in  the  coolant  channel,  large  pumping  power,  etc.  With  surface  boil¬ 
ing,  high  heat  transfer  coefficients  may  be  obtained  at  much  lower  veloci¬ 
ties  and  consequently  with  less  pressure  loss,  pumping  power  etc.  As  a 
safety  feature  it  should  be  noted  that  a  small  change  in  flow  rate  is  not 
so  critical  as  In  forced  convection  (provided  one  Is  not  too  close  to 
"burnout",  which  will  be  discussed  below). 

The  general  form  of  the  boiling  heat  transfer  curve  for  pool  boil¬ 
ing  can  be  found  in  any  standard  text  on  heat  transfer.  In  the  nucleate 
boiling  region  it  is  noted  that  the  heat  load  rises  very  rapidly  with  a 
small  change  in  ,  thus,  giving  a  large  heat  transfer  coefficient. 

A  point  is  reached  called  the  "burnout  1  point,  or  the  point  of  burnout 
heat  load  Qg0  .  For  a  condition  of  constant  Q  when  this  heat  load  is 

exceeded,  there  is  a  very  large  increase  in  AT^  .  The  hot  metal  surface 

melts  or  vaporizes  at  this  temperature,  and  thus,  the  term  burnout  point 
arises. 

Two  calculations  are  required  in  the  design  of  a  system  for  surface 
boiling  heat  transfer.  Calculations  must  first  be  made  to  determine  the 
surface  temperature  required  to  transfer  the  given  heat  load  at  a  given 

coolant  flow  rate  in  a  given  channel.  The  burnout  heat  flux  Qto  ,  must 

then  be  checked  to  see  that  for  these  conditions  it  exceeds  the  given  heat 
load. 


Since  the  heat  load  is  essentially  fixed  by  the  gas  temperature  on 
the  air  side  of  the  liner  the  first  calculation  reduces  to  a  determination 
of  the  liner  surface  temperature  for  a  given  bulk  coolant  temper _ure, 
channel  geometry  and  flow  rate.  The  method  of  Rohsenow  as  described  by 
Kreith  (  16  )  is  used  for  this  calculation.  The  heat  flux  is  separated 
into  two  parts,  the  convective  heat  flux,  and  the  boiling  heat  flux. 


Then 

onvtct/v*  ~t~  mourn 9 

(  199  ) 

Now 

^convwAv  r  *  ^ c  ( ) 

(  200  ) 

thus  bc(Ts-Th)  (  20,  ) 
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In  these  equations,  hc  is  calculated  from  equation  (  196  ).  Using  the 
data  of  HcAdams,  et  al  (  17)  the  boiling  heat  flux  may  be  calculated 
using  the  equation 


i.efo 


0.2 


(  202  ) 


Then 


QToW  =  0.2 4-  Kc(T,-Tte) 


(  203  ) 


Using  the  results  of  the  convective  heat  transfer  calculations,  T^,  hc, 
and  Tsat  are  known.  Equation  (  203  )  is  then  solved  by  trial  and  error 
for  T$,  the  surface  temperature. 

Equations  for  predicting  QgQ  for  various  flow  conditions  do  not 
show  very  good  agreement.  A  summary  by  Kreith  f  16  1  of  recent  correla¬ 
tions  gives  the  following  equations: 


4-*(0  4  4Aoo  ^ 


?! 


(  206  ) 


Qi 


a© 


7000  U 


(  205  ) 


?bo  -  CCtx.o’^at^) 


(  206  ) 


where  21"^= 


U  = 
Gr  = 


^sat  "  ^b 

velocity,  ft/sec 

mass  flow  rate,  lbm/ft2  hr 


c  and  m  are  tabulated  functions  of  the  static  pressure. 

Another  correlating  equation  by  McGill  and  Sibbitt  (  18  )for  a  tube 
with  L/D  =21  is 


a/,„xs*c. 


(  207  ) 


To  show  the  lack  of  agreement  in  the  above  equations  for  0go  assume 
for  equations  (  204  )  and  (  205  ) 

=  100  °F 

CJ  -  20  ft/sec. 
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These  values  are  In  the  range  supposedly  correlated  by  equations  (  20k  ) 
and  (  205  ). 

From  equation  (  204  ): 

6  2 

qBo  =  2.38  x  10  B/hr  ft4 

From  equation  (  205  ): 

A  2 

0Bo  =  3. 11*  x  10°  B/hr  ft 

From  equation  (  206  ) ,  using  the  same  values  for  u  and  A  Tsub  and  a 
pressure  of  500  psla 

QBo  =  0.806  x  I 06  B/hr  ft2 

From  equation  (  207  ) 

qBo  *  4.0  x  10^  B/hr  ft2 

The  total  variation  Is  thus  from  800,000  to  4,000,000  B/hr  ft2,  a  factor 
of  5. 


As  another  approach,  the  effects  of  sub-cooling,  velocity,  and  static 
pressure  may  be  considered  to  be  Independent  and  another  equation  proposed. 

Lowering  the  pool  temperature  below  the  saturation  temperature  increases 
°Bo  for  pool  boiling.  For  degassed  distilled  water,  lowering  the  pool 
temperature  to  A  Tsub  *  50*F  doubles  QBo  and  lowering  I t  to  A  T  b  *  100“F 

quadruples  .  An  equation  which  correlates  the  date  of  M.  E.  Ellion  (as 
given  by  Krelth,  not  taken  from  original  report)  Is, 

O#o  *  4*'0* +  6000  Al«ob  +  60  ATsub  (  208  ) 


or  -/  »  v 

w  4xio  ^1  +  O-OllS1  AT,oW-eO.OOOZ  )  (  209  ) 

Equation  (  209  )  is  plotted  in  Figure  39*  Ellion  also  gives  data  for  water 
with  some  detergent  added  to  lower  the  surface  tension,  the  same  type  effect 
as  would  be  obtained  from  gas  in  the  water.  An  equation  which  correlates 
this  data  well  is , 

$5o*  O.OV1  ATJ(>b)  (  210  ) 
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Equation  (  210  )  is  plotted  in  Figure  39. 

Experiments  by  M.  T.  Ciechelli  and  C.  F.  Bonilla  (  19  )  using  various 
liquids  showed  that,  for  pool  boiling,  increasing  the  static  pressure 
increased  QBq  up  i0  a  pressure  of  one-third  the  critical  pressure.  An 
equation  which  correlates  their  data  well  is, 

*'lot  ['*■*“(%) -*'(%)  ]  (2") 

where  Ps/p  =  ratio  of  static  to  critical  pressure.  Equation  (  211  ) 
is  plotted  in  Figure  40. 

Combining  equations  (  210  )  and  (  21 1  )  and  using  as  a  conserva¬ 

tive  choice  for  the  velocity  dependence  gives 

qSe*4x^i+o.o*»7ATlt%)[.f  ]ci  (  212  ) 

This  equation  is  plotted  for  three  values  of  ps  and  four  values  of  u  in 
Figures  41  ,  42,  and  43.  For  the  values  u  =  20,  ps  =  500,  and  &  TSub  =  100 
as  used  in  previous  calculations,  this  equation  gives 

QBo  =  x  I05  (i  +  2.17)  (1  +  4.63  -  2.64)  (2.71) 
or  Q0O  =  10.4  x  10  B/hr  ft 


Though  this  value  is  considerably  higher  than  the  highest  value  previously 
calculated  from  these  conditions  there  is  some  indication  it  might  be 
achieved.  Equation  (  205  )  was  proposed  by  Gunther  (  20  ) .  and  for  these 
flow  conditions  gave  0~  =  3. 14  x  10®  B/hr  ft  .  In  his  article  Gunther 
noted  that  unpublished  data  by  McKenney  for  a  tube  with  .012"  walls  gave 
burnout  limits  "well  above"  those  predicted  by  equation  (  205  )•  This  was 
attributed  to  the  increased  heat  capacity  of  a  heavier  tube  which  handled 
transients  better.  Since  the  nozzle  liner  will  have  a  minimum  thickness  of 
approximately  .012"  it  is  possible  that  equation  (  212  )  gives  a  better 
prediction  of  0,Bo  than  those  cited. 

Since  equation  (  212  )  is  unproved  and  there  is  a  large  scatter  in  QBo 
predictions,  using  suggested  equations,  several  values  of  QBo  are  calcu¬ 
lated  in  the  computer  study  so  a  comparison  may  be  made. 
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Fig.  42  Plot  of  Equation  (212);  Pressure  =  100  psio 
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Fig.  43  Piet  of  Equation  (212);  Pressure  =  300  ptia 
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CONCLUSIONS 

Studies  of  the  backside  cooling  were  made  using  the  computer  program. 
The  results  of  three  cases  of  particular  Interest  are  shown  on  pages  1 1*7 
through  155.  The  nozzle  used  for  the  studies  had  an  entrance  diameter  of 
1.5  inches,  a  46  degree  converging  section,  a  .25  inch  diameter  throat  and 
a  7.5  degree  diverging  section.  The  heat  load  distribution  assumed  is 
given  under  Q  in  the  Input  data. 

Since  subcooling  and  velocity  of  the  coolant  both  affect  the  boiling 
burnout  heat  load,  a  flow  rate  of  20  Ibs/sec  gave  good  results  for  the 
cas  s  with  boiling.  A  flow  rate  of  25  Ibs/sec  was  used  for  the  third  case 
whe.~  boiling  was  not  allowed. 

The  first  case  presented  shows  the  conditions  necessary  to  achieve 
approximately  the  required  boiling  burnout  heat  load  with  equation  (  212) 
used  for  0go.  These  are  summarized  as  follows: 

Coolant  Flow  Rate  »  20  Ibm/sec2 

Coolant  Annulus  Area  =  0.35  in 

Coolant  Velocity  »  133  ft/sec 

Maximum  Surface  Temperature  =  489  #F 

Surface  Coefficient  *  76,000  B/hr  ft2  ®F. 

The  second  case  presented  shows  the  conditions  necessary  to  achieve 
the  required  boiling  burnout  heat  load  with  equation  (  205  }  used  for  0go. 
These  are  summarized  as  follows: 

Coolant  Flow  Rate  =*  20  lbm/sec 

Coolant  Annulus  Area  =  0.200  in2 

Coolant  Velocity  *  232  ft/sec 

Maximum  Surface  Temperature  *  498  °F  _ 

Surface  Coefficient  =  75,000  B/hr  ft  ®F. 

The  third  case  presented  shows  the  conditions  necessary  to  achieve 
forced  convection  heat  transfer  at  all  sections.  These  are  summarized  as 
fol lows : 


Coolant  Flow  Rate  *  25  lbm/sec 

Coolant  Annulus  Area  3  0.160  in2 

Coolant  Velocity  =  363  ft/sec 

Maximum  Surface  Temperature  «  478  ®F 

Surface  Coefficient  =  78,000  B/hr  ft2  *F. 

Several  conclusions  may  be  drawn  from  these  calculations  based  on  the 
heat  load  and  coolant  used: 

(1)  The  maximum  heat  transfer  coefficient  obtainable  while  satis¬ 
fying  the  burnout  requirements  Is  80,000  B/hr  ft2oF. 

(2)  No  appreciable  reduction  in  surface  temperature  may  be  obtained 
by  boiling  as  compared  with  high  velocity  forced  convection. 
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(3)  The  required  Inlet  pressure  to  prevent  boiling  Is  ~  1620  psia 
et  e  flow  rete  of  25  Ibm/sec  end  e  cooling  ennui  us  of  0.16  In  . 

(4)  The  required  inlet  pressure  with  boiling  is  ~  325  psie,  and 
little  is  gained  by  using  equation  (  212  )  rather  than  equation 
(  205  )  for  the  burnout  calculation. 

(5)  There  is  room  for  doubt  that  this  heat  load  may  be  transferred 

by  boiling  without  burnout,  due  to  the  feet  that  equation  (  204  ) 
predicts  burnout  In  all  cases  calculated  where  boiling  is  allowed. 

It  Is,  therefore,  recommended  that  provisions  be  made  to  supply  the 
coolant  water  at  sufficient  pressure  to  prevent  surface  boiling  and  the 
coolant  annulus  be  designed  for  high  velocity,  high  pressure,  forced  con¬ 
vection  heat  transfer. 
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NOZZLE  STRESS  PROGRAM 


This  computer  program  is  written  In  tha  Fortran  II  languaga.  It 
accapts  data  which  dascriba  tha  Itnar,  rib,  and  slaava  as  to  dimensions, 
configuration,  materials,  haat  load,  ate.  It  than  calculates  tha  princi¬ 
ple  stresses,  octahedral  shear  stresses,  etc.  as  described  in  the  Output 
section  below.  The  equations  used  for  the  calculations  are  given  in  the 
Stress  Analysis  discussion  under  Case  I  and  Case  2. 


INPUT  DATA 


Flrit  Card 

The  first  data  card  Is  a  title  card  with  up  to  forty-nine  alphabetic 
or  numeric  characters  to  identify  the  case  being  studied. 


Second  Cord 

The  second  data  card  gives,  according  to  FORMAT  (E 1 2. 1 ,  EI2.I,  EI2.I, 
El 2. 1 .  EI2.I,  El 2. 1 )  the  following  data  in  order: 

PT  -  The  airslde  pressure  at  the  throat,  psi. 

Q  -  The  liner  heat  load,  Btu/hr  ft2. 

Oil-  The  liner  inside  diameter,  inches. 

DS  -  Liner  Inside  diameter  at  nozzle  inlet,  Inches. 

DC  -  Coolant  annulus  width,  radially,  inches. 

TCOOL  -  Coolant  bulk  temperature,  °F. 


Third  Card 

The  third  date  card  gives,  according  to  FORMAT  (E 12.1,  El 2.1,  EI2.I, 
EI2.1,  El 2. 1 ,  EI2.I)  the  following  data  in  order: 

Cl  -  Liner  thermal  conductivity  at  temperature  Til,  B/hr  ft  *F. 

C2  -  Liner  thermal  conductivity  at  temperature  TI2,  B/hr  ft  *F. 

C3  -  Liner  thermal  conductivity  at  temperature  T 1 3 ,  B/hr  ft  *F. 

Al  -  Liner  thermal  expansion  coefficient  at  temperature  Til,  in/ln  *F 

A2  -  Liner  thermal  expansion  coefficient  at  temperature  T12,  in/in  *F 

A3  -  Liner  thermal  expansion  coefficient  at  temperature  TI3,  In/in  *F 


Fourth  Cord 


The  fourth  data  card  gives,  according  to  FORMAT  (EI2.I,  EI2.I,  EI2.I, 
EI2.1,  EI2.1,  EI2.I)  the  following  data  in  order: 


El  -  Liner  modulus  of  elasticity  at  temperature  Til,  psi. 

E2  -  Liner  modulus  of  elasticity  at  temperature  T12,  psi. 

E3  -  Liner  modulus  of  elasticity  at  temperature  Tl 3 ,  psi. 

U  -  Poisson's  ratio  for  the  liner. 

ES  -  Sleeve  modulus  of  elasticity,  psi. 

ER  -  Rib  modulus  of  elasticity,  psi. 
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Filth  Cord 

The  fifth  data  card  gives,  according  to  FORMAT  (E12.I,  EI2.1,  E12.I, 

El 2.1,  EI2.I,  El 2.1)  the  following  data  in  order: 

HT  -  The  limiting  liner  high  temperature,  *F. 

HW  -  The  coolant-liner  heat  transfer  coefficient,  B/hr  ft^  *F. 

Til-  First  material  property  reference  temperature,  *F. 

TI2-  Second  material  property  reference  temperature,  eF. 

Tl 3-  Third  material  property  reference  temperature,  *F. 

TYPE-A  number  from  1-6  inclusive  which  specifies  the  configuration 
for  which  calculations  are  wanted.  These  are  as  follows: 

1.0  :•  circumferential  ribs  only 
2.0  *  longitudinal  ribs  only 

3.0  *  unrestrained  liner  and  circumferential  ribs 
4.0  a  unrestrained  liner  and  longitudinal  ribs 
3.0  x  unrestrained  only 

6.0  *  unrestrained,  longitudinal  ribs  and  circumferential  ribs. 


Sixth  Card 

The  sixth  date  card  gives,  according  to  FORMAT  (E12.I,  EI2.I,  E12.I, 
EI2.1,  EI2.1)  the  following  data  in  order: 

TR  -  Rib  thickness,  inches. 

0L  -  Distance  between  circumferential  ribs,  inches. 

AR  -  Cross-sectional  area  of  circumferential  ribs,  inches. 

THK-  Stepwise  increase  in  thickness  of  liner  for  calculations,  inches. 
TlO-  Starting  thickness  of  liner  for  calculations,  inches. 

Seventh  Card 

The  seventh  card  gives,  according  to  FORMAT  (E12.I,  E12.1)  the  follow¬ 
ing  data  in  order: 

OR  -  Distance  between  longitudinal  ribs,  inches. 

ST  -  Liner  thrust  stress  for  circumferential  ribs,  psi. 


PROGRAM  OPERATION 

The  program  operates  automatically  using  the  above  input  data  for  each 

case. 


Starting  with  a  liner  thickness  given  by  TlO  this  thickness  is  increased 
in  seven  steps  by  the  amount  THK.  The  sleeve  thickness  is  constant  at  1.0 
inches,  as  the  configurations  with  ribs  will  be  fully  restrained. 

For  each  liner  thickness  the  calculations  for  stresses  are  made  for 
zero  coolant  pressure  and  then  for  coolant  pressure  equal  to  air  pressure. 
Other  pressure  values  may  be  used  by  altering  statement  number  49.  If  the 
equal  pressure  case  is  not  desired  the  program  is  run  with  SENSE  SWITCH  2  ON. 
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The  program  is  so  Mritten  that  if  the  limiting  high  liner  temperature 
HT  is  exceeded  at  the  Inside  surface  the  maximum  possible  liner  thickness 
which  does  not.  cause  an  over  temperature  is  calculated  and  used  for  the 
last  thickness. 

To  speed  up  the  running  of  the  program  for  cases  which  do  not  involve 
circumferential  ribs  it  is  run  with  SENSE  SWITCH  3  ON.  This  by-passes  the 
summing  operations  necessary  for  circumferential  ribs. 


OUTPUT 


Unroxtrainod  NoxzU: 

PB  -  Coolant  pressure,  psi. 

Tl  -  Liner  thickness,  inches. 

Tl  -  Liner  inside  surface  temperature,  *F. 

TB  -  Liner  outside  surface  temperature,  *F. 

SPCI-Prlnciple  circumferential  stress,  inside  surface,  psi. 

SPC2-Princtple  circumferential  stress,  outside  surface,  psi. 

SPLI -Pri nclple  longitudinal  stress,  inside  surface,  psi. 

SPL2-Prlnclple  longitudinal  stress,  outside  surface,  psi. 

SIH  -Octahedral  shear  stress,  inside  surface,  psi. 

S2H  -Octahedral  shear  stress,  outside  surface,  psi. 

Longitudinal  Rib  NoxzU: 

PB  -  Coolant  pressure,  psi. 

T2  -  Sleeve  thickness. 

Tl  -  Liner  thickness,  inches. 

Tl  -  Liner  inside  surface  temperature,  *F. 

TB  -  Liner  outside  surface  temperature,  °F. 

DR  -  Distance  between  ribs,  inches. 

SIPCI-PrincIple  circumferential  stress,  inside  surface  over  rib,  psi. 
S!PC2-Principle  circumferential  stress,  outside  surface  over  rib,  psi. 
S2PCI-Principle  circumferential  stress,  inside  surface  between  ribs,  psi. 
S2PC2-Principle  circumferential  stress,  outside  surface  between  ribs,  pse. 
SPLI-  Principle  longitudinal  stress,  inside  surface,  psi. 

SPL2-  Principle  longitudinal  stress,  outside  surface,  psi. 

SR|2-  Pressure  between  liner  and  rib,  psi. 

SSC22-S!eeve  bending  stress  (negligible),  psi. 

SLS  -  Sleeve  longitudinal  stress  (negligible),  psi. 

SIHI-  Octahedral  shear,  inside  surface  over  rib,  psi. 

S|M2-  Octahedral  shear,  outside  surface  over  rib,  psi. 

$2Ml-  Octahedral  shear,  inside  surface  between  ribs,  psi. 

S2M2-  Octahedral  shear,  outside  surface  between  ribs,  psi. 

Circumforontiol  Rib  NoxzU: 

PB  -  Coolant  pressure,  psi. 

T2  -  Sleeve  thickness,  inches. 

Tl  -  Liner  thickness,  inches. 
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Tl  -  Liner  inside  surface  temperature,  *F. 

TB  -  Liner  outside  surface  temperature,  *F. 

OL  *  Distance  between  ribs,  inches. 

SILCi-Principie  circumferential  stress,  inside  surface  over  rib,  psi. 
SILC2-Principle  circumferential  stress,  outside  surface  over  rib,  psi. 
S2LCI-Principle  circumferential  stress,  inside  surface  between  ribs,  psi. 
S2LC2-Princlple  circumferential  stress,  outside  surface  between  ribs,  psi. 
SILLI-Princlple  longitudinal  stress,  inside  surface  over  rib,  psi. 
SILL2-Prlnciple  longitudinal  stress,  outside  surface  over  rib,  psi. 
S2LLl-Principle  longitudinal  stress,  inside  surface  between  ribs,  psi. 
S2LL2-Prlncip1e  longitudinal  stress,  outside  surface  between  ribs,  psi. 
SIHI  -Octahedral  shear,  inside  surface  over  rib,  psi. 

S I M2  -Octahedral  shear,  outside  surface  over  rib,  psi. 

S2MI  -Octahedral  shear,  inside  surface  between  ribs,  psi. 

S2H2  -Octahedral  shear,  outside  surface  between  ribs,  psi. 

B  -  Parameter  used  In  stress  calculation  (See  derivation). 

P  -  Force  per  inch  of  length  on  ribs,  lb/in. 

PI  I  -  Bending  stress  on  liner  at  rib,  psi. 

SR2I-  Pressure  between  liner  and  rib,  radial,  psi. 

The  program  and  a  sample  output  for  longitudinal  rib  calculation  follow 
on  page  1 14. 
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c  c 

c 

c 

c 

c 

c 

c 

c 

c 

90 

91 

92 

93 
100 
101 
102 
103 
1  04 

1  1  o 
1 1 1 
1  12 
1  13 
1  14 
1  15 
120 
121 
126 
126 

129 

130 


2 

0 


3 

7 


H • C •  R0LAN0  *  ♦  NOZZLE  STRESS  PROGRAM 

FIRST  INPUT  IS  TITLE  WITH  UP  TO  49  ALPHAMERIC  CHARACTERS 
TYPE  6*0  ■  uNRESTRA I NED* LONG  I TUDINAL •  AND  CIRCUMFERENTIAL  RIBS 
5*0  «  UNRESTRAINED  ONLY 

4*0  ■  UNRESTRAINED  ANU  LONGITUDINAL  R I  OS • SENSE  SWITCH  3  ON 
3.0  ■  UNRESTRAINED  AND  CIRCUMFERENTIAL  RIBS 
2*0  ■  LONGITUDINAL  RIBS  ONLY*  SENSE  SWITCH  3  ON 
1*0  ■  C IRCUMFERENT I Al  R I  PS  ONLY 
SENSE  SWITCH  2  ON  BYPASSES  AUTOMATIC  PRESSURE  LUuALlZATJQN 
FORMAT  (49m  ) 

FORMAT (E 12. 1 12. 1 *E12.1*E12.1*E12.1 *E12.1 ) 

FORMAT  ( E  1  2  •  1  .t  12. 1  *C.12. 1  .El  2.  I  «c  12*1  ) 

FORMAT ( I OHI NPUT  DATA) 

FORMAT (6X.2HPB. 12X.2HT1 .9X*2hT1 *VX*2HTB) 

FORMAT ( 3X.4HSPC I « 7X • 4HSPC2 • 9X . 4HSPL 1 *8X • 4HSPL2 • 9X • 3HS 1 M • 9X *3HS2M) 
FORMAT ( 19HUNRE STRAINED  NOZZLE) 

FORMAT (F12*3.F12*3.F12«3.F12*3.F12*3*F12.3) 

FORMAT (FI  2.3. F 12.3.F 12.3 *F 12.3) 

FORMAT (23HLONG1 TUD I NAL  RIB  NOZZLc ) 

FORMAT (6X.2HPB. 12X.2HT2, 1 0X.2MT 1 .8X.2HT 1  *  1 0X.2HTb. 13X*2hDR ) 

FORMAT ( 3X.5HSI PCI  • 7X . bMS 1  PC 2 • 7X • 3HS2PC 1  .7X.6HS2PL2) 

FORMAT (3X .4HbPL 1 .8X . 4HSPL 2 . BX • 4H3R 1 2 . 7x • 5HSSC22 • 9X  ♦ 3HSLS ) 

FORMAT (3X.4HS1M1  . 8X • 4HS 1  M2 . 8X . 4Hb2 M 1  .bX.4HS2M2 ) 

FORMAT  (F  |  £.  U.J.h  i 

FORMAT <26mC IRCUMFERENT I AL  R.O  NOZZLE) 

FORMAT (6X.2HPB. 12X.2HT2. 10X. 2HT I  »8X.2MTI  *10X.2HTB»  I3X.2HDL  ) 

<6X. 1HB' 1  I  X.  1 HP.9X. JHP1  1  .bX.4HSR21  ) 

(5X.5HS1LL1 . 7X.DHSll.L2.7X.bH5)2LLl  *7X.5HS2LL2> 

(5X.5HS1LCI  .  7X.5HSILC2. 7X.5HS21.C  1 .7X.Sms2LC2> 

( E  1  2  •  1  •  E  1  2  •  1  ) 

1*1 .100 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
DO  51 
READ  90 
PUNCH  90 

READ  91  ..PT.Q.Ul  1  *OS*UC.TCOOL 
READ  91 *C1 *C2.C3.A1 , A2.A3 
READ  91 .E 1 *E2*E3.u*Eo.tR 
READ  91 .HT .Hw. T I  1  . T I2.T 1 3. TYPE 
READ  92. TR.DL. AR.THK , T1 0 
READ  130.  DR*  ST 
PUNCH  93 

PUNCH  91 .PT.0.D1 1 *DS*OC *TCOOL 
PUNCH  91  .Cl  .C2.C3.A1  •  A2  »  A  3 
PUNCH  91  .E 1  * E2.E3.U*ES*ER 
PUNCH  91  .HT.HW. T I  1  . T I2.T! 3*TYPE 
PUNCH  92  * TR.DL • AR. THK  .T 1 0 
PUNCH  130*  DR*  ST 
COUNT  »  0. 0 
T  1  «  T  1  0 
DO  52  J« 1 *0 
IF  < J-l )  8.8*2 


T  1  «T  1  4-THK 
D1 2*DI 1 +2 • ♦ T 1 
DRCAL*SQRTF (D1 2*T1  ) 
IF  (DR-DRCAL)  3*7*7 

dr«orcal 

RL»LOGF (D12/D!  1  ) 


1 

2 

3 

4 

5 

6 

7 

8 
8A 

9 
10 
1  1 
12 
1  3 

14 

15 

16 
1  7 
18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
4t j 

49 

50 

51 

52 

53 

54 

55 
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0L«AT»0L/TR 

56 

tlrat*tr/dl 

57 

TRL-DC/6. 

58 

IF  ( TRL-OW/4 •  >  4*4*5 

59 

* 

GO  TO  6 

60 

5 

TRU-DR/4. 

61 

6 

0T1 «0*D1 1*RL/ (24.*C2 » 

62 

TB*TC00L+0/HW 

63 

T !  *TB+DT 1 

64 

TM*DT1 /LOGF ( T I/T8 ) 

65 

IF  (TM-T12)  10*30*20 

66 

10 

DO  15  K*1 *3 

67 

TC*<TM-TI  1  >/(TI2-TI  1  ) 

68 

CM»Cl +TC*(C2-Cl ) 

69 

OTt ■0*D1 l*RL^(24.*CM) 

70 

TI«TB+OT! 

71 

TM»0T1/L0GF(TI/TB) 

72 

1 F  (TM-TI2)  15*30*20 

73 

15 

CONTINUE 

74 

GO  TO  30 

75 

20 

DO  25  L* 1 *3 

76 

TC»<TM-TI2)/(TI3-TI2) 

77 

CM»C2+TC*(C3-C2 ) 

78 

DTI «0*D1 1*RL/(24.*CM) 

79 

T  I  *TB+DT 1 

BO 

TM»DT1/L0UF(TI/Te J 

81 

IF  (TM-TI2)  10.30*25 

82 

25 

CONTINUE 

83 

GO  TO  35 

84 

30 

TC*  ( TM—T I 1  )  / ( T I 2  — T I  1  ) 

85 

AL»AI4TC*<A2-A1 ) 

86 

EL*EM>TC*(E2-E 1 ) 

87 

IF  (TI-HT)  40*40,32 

88 

32 

C0UNT«C0UNT*1  • 

89 

33 

T1 *T1-.001 

90 

012*01  1  +2«*T  1 

91 

RL*L0GF <012/01 I > 

92 

DT 1 *0*01 l*RL/<24.*CM> 

93 

T  1  • TB+DT 1 

94 

IF  (TI-HT)  34.34*33 

95 

34 

IF  (COUNT-1.5)  40*40*38 

96 

38 

GO  TO  51 

97 

35 

TC  • ( TM—T 1 2 )/ <  T  I 3— T I 2 ) 

98 

AL*A24TC*<A3-A2 ) 

99 

El_*E2*TC*<E3-E2 ) 

100 

IF  (TI-HT)  40*40*32 

101 

40 

T2*  1 .0 

102 

OTM«TM-TCOOL 

103 

022*01  2-f  2.*DC 

104 

023*D22*2.*T2 

105 

RMN* ( D 1 1 ♦ T I )/2. 

106 

IF  (SENSE  SWITCH  3)  81,41 

107 

41 

B-SQRTF (SURTF (2.73/ ( ( RMN**2 ) * ( T 1**2 >  >  > I 

108 

0*0. 

109 

F«|  • 

1  10 

SUM! «-.5 

1  1  1 
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SUM2»0<  112 

00  80  LM=1 <8  113 

ARG1=B*D*DL  114 

ARG2=B*F*DL/2.  115 

SUMI»SUM|+  (2.718**(-ARGI  )  )*(C0SF(ARG1  )  +S  INF  (  ARG1  )  )  116 

SUM2«SUM2+(2.7I8**(-ARG2)  )  *  (  COSF  (  ARG2  )  +s  INff  (  ARG2  )  )  117 

0*0+1.  118 

F*F+2.  119 

80  CONTINUE  120 

A*B»OL  121 

CH8*.3*( (2.718**A)+<2.718**(-A) II  122 

CB*COSF (B*DL )  123 

SH8«.5*< (2.7!8*»A)-<2.718**(-A) >  >  124 

SB*SINF (B*OL )  125 

SI  P8*0 .  126 

012S0* (012**2 )  127 

Dl 1  SO* (Oil **2 )  128 

022S0* (022**2)  129 

023S0*( 023**2)  130 

DSS0*(0S**2>  131 

21 ■ (D12SQ+01 1  SO  >/ (D12SQ-D1 1  SO)  132 

22*012S0/ ( D1 2SQ-01 ISO)  133 

Z3«01 1 SO/ (01 2S0-D1 1  SO )  134 

Z4 • ( D22S0+0 1 2  SO ) / ( 022S0-D 1 2SQ )  135 

ZS=022S0/ (D22SO-D1 2S0 )  136 

Z6=Dl2Ea/(D22SQ-0l2S0 )  137 

27= ( D23SQ+D22SQ ) / (023S0-022SQ )  138 

26= (01 2 $0-01  ISO)/) Q23SQ-D22SQ )  139 

Z9XDSS0-DI  1  SO ) / (0 l 26Q-D1 ISO)  140 

Z10»  <0SSQ-D12Sa>/(012SQ-01 ISO)  141 

FI =AL*0Tl*EL/<2.*< 1 .-U)*RL )  142 

F2«| «-2.*Z3*RL  143 

r4=l .-2.*Z2*RL  144 

DO  50  N« 1 <2  145 

IF  (N-l )  67.67,68  146 

66  IF  (SENSE  SWITCH  2)  69.67  147 

69  GO  TO  50  | 48 

67  IF  (TVPE-2.5)  70.60.60  149 

70  I F ( TYPE- 1 <5 )  46.42.42  150 

60  PS»PT/.526  151 

26  SPC1 =Z1*PT-2.*Z2*P8+F1*F4  152 

SPC2=2<*Z3»PT-Z1*PB+F1*F2  153 

SPL1 2* (PS/ (OS-Dl 1 ) )*( ( <328*DS-D1 1 >*Z9)  | 54 

SPL 1 3= ( PS/ ( OS-D 1 l ) )*.3148*( (0S**3)-(D1 1**3) )/(D12SQ-01 ISO)  155 

SPL 1 =F 1 *F4— ( SPL 1 2+SPL 1 3 ) +PB*Z 10  156 

SPL2*F 1 *F2- ( SPL 1 2  +  SPL I  3 ) +PB*Z 1 0  157 

S 1 M«  SQRTF < ( (SPC1-SPL1  )**2 )♦ ( (SPC l +PT ) **2 ) + (  ( SPL t  +PT )**2 ) )/3.  158 

S2M«S0RTF( ( (SPC2-SPL2 1**2 )♦ ( ( SPC2+P6 )**2 )♦ ( ( SPL2+PB ) »*2 ) )/3.  159 

PUNCH  102  160 

PUNCH  100  161 

PUNCH  104.PB.TI .TI ,TU  162 

PUNCH  101  163 

PUNCH  103.SPC1 .SPC2.SPLI .SPL2.S1M.S2N  164 

29  IF  (TYPE-5.5)  61.42.42  165 

61  IF  (TYPE-4.5)  62.49.49  166 

62  IF  (TYPE-3.5)  46,42 .42  167 
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42  PlN*2#*Z3*PT-2i*PB-2?*EL*D22*PB/<tS*D12)  166 

PI  D*Z1+Z7*EI_/ES  169 

P 1 =P 1 N/P | D  170 

P2  *P 1 #D 1 2/D22  171 

SI  1  “Z 1 *PT-2 • *Z2*  »PB+P 1  )  172 

S12*2.+23*PT-Z1 * (PB+Pl )  173 

S22«27*(Pb+P2  )  1  7 '♦ 

PlPaFl *F2/ (Z1 +2 7*EL/ES>  1  /b 

SI  1 P«F1 *F4-2.*Z2*P1P  176 

S12P*F1 #F2-21 *P1P  177 

S22P*Z7#P1P*D12/D22  1 7ti 

SB1 1 **S*P1 * < 1DR/T1 ) **2 )  179 

SB  1 22  -SB 1  1  160 

Sbl  1P=.S*P1P*<  (DR/T1  >**2)  161 

SB  1 2P  =  -SB  1 1 P  162 

SRI  2  =  -<Pl  -fPlP  )*  (DP/TP  )  163 

SLT 1  1  *F 1 *F4  164 

SLT  1  2*F  1  #F  2  lt>5 

S 1 L  T  s-Al*DTM*EL/ ( 1 •♦Z6*EL/tS )  166 

S22PP  =  2  7*  <Pb+ ( J12/022 )* (PI +PIP  >  )  167 

SB22s-*6*(022/Q12)*(Pl+P!P)*( UW/T2>**2)  166 

Sd222*-«5*SB22  169 

S2LT =AL*DTM#tL/ (EL/tS+l ./Z6)  1 9Q 

S1PC1 =S1 1 ♦SI 1P+S61 1 ♦Sbl IP  <  191 

SlPC2*S12  +  S12P-fSB12  +  Sbl2P  192 

S2PC  1  =%  1  1  t-bl  IP-f  •b»5612+«b#SH12P  193 

S2PC2  =  C  12  +  S12P*  •  b*SL1l  1  ♦  •£5*Srt  1  IP  194 

SPL 1 =SLT 1 1 ♦SILT  195 

SPL2*SLT t2*MLT  196 

SSC22*S22P>S6222  197 

SlS=  S2LT  190 

51  Ml  *SQPTF  <  (  (  SI  PC  1  -SPL1  )**2  >♦ (  (  Si  PC  l  «-PT  >**2  >  +  <  (  SPL  1  +PT  >**2  )  >/->•  199 

S  1M23  =  SCIRTF  (  (  (  S  1  PC2  -  SPL2  >  *  *c  )  M  <  S 1  PL2-SR 1  2  )  **2  >♦<  (  SPL  2“S«  1  2  )  **2  )  )  200 

S1M2*S1M2 J/3#  201 

52  v  1  s  SOWTF  (  <  (  -SPLl  l  **2  >  ♦(  <  S2PC1  +^T  )  **2 )+ (  l.'PLl+PT  )**2  )  )/ 3*  2Cc 

S2M2  =  S9RTF  (  (  (  S2PC  2- 6»  'l.  2  )  »  »  2  )  «•  (  (  bHl.?  +  Po  )**2  )♦  (  (  S2PL2  +  Pr> )  *  #2  )  )/ J*  ^0  r 

PUNCH  110  204 

PUNCH  1  l  1  c:Oi 

PUNCH  1  f  •  3 ♦ PB  ♦  T  2  •  T  1 •  T  I  • T  6 • DR  20c 

PUNCH  1 12  207 

PUNCH  104,MPC  1  *:.1®C2.S2PC  1  ♦L-.2PC2  cOe 

PUNCH  1 1 3  209 

PUNCH  l  15*SPL1  t6PL2»SHi2f  L->C22t.->L.>  21  C 

PUNCH  114  211 

PUNCH  1  04  «  S  1  Ml  •_  1  M2  f  ot'M  1  ,  -x.tfc  2  \  c 

|F  <TYPc.-b«b)  64  «  46 • 46  21. 

64  CO  TO  49  2W 

46  PM  =9T-PH«-AL*UTM*tL*  T  1  /*,N.N  2H 

P-PNl/ld*iUMl+Tl/ARl  2lC 

PMC  =  P*  ( SHO—SB ) /  (4«*6»  KHtJ-Co)  )  21 

PI  1  =6«*PM0/ < T 1 **? )  2 1 k 

SlLLl  =E  1  *F  4-fPl  1  -fST  2  1 1 

SlLL2=r  1  #F2-Pi  1  +ST  22f 

s2ll i -r i *\ 4-,h»pi i +st  22: 

S2LL2  =  F  1  ♦F2+«b*Pl  1-fST  22< 

SW2 1 =-P/Tk 
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S 1LC 1 «P*RMN/Aft-AL*OTM*EL+F 1 *F A  224 

S1LC2«P*RMNXAR-AL*0TM*EL+Fl *F2  225 

SLC2  I  •  <PT-PB>»RMN/Tl+Fl*F2-P*b*RMN#i»uM2/Tl  225 

SLC23« (PT-PB  >*RMN/T l +F 1 »F4-P*b*RMN*SUM2/T I  227 

S2LC 1 “SLC23  22B 

S2LC2-SLC2 1  229 

SIM1»S0RTF( < (S1LL1-S1LC1 1**2 1+1 ( S 1 LL I +PT ) »«2 >* ( ( S 1 LC 1 +PT ) **2 * >✓3.  230 

SIM*.SQRTF< < ( S l LL2-S1 LC2  1 **2 1  +  { (SILL2-SR21 )«*2  >  +  ( (S1LC2-SR21  )**2)  )  231 

SlM2*SlM4/3.  232 

S2M1  «SQRTK  (  <  (S2LL1-S2LC1  >*«2l+<  I  52LL  1  +PT  ) **2  )  ♦  I  (  S2LC 1  4-PT  ) **2  )  >/3.  233 

S2M2«SQRTF( <  <S2LL2-S2LC2>»*2>+( ( S2LL2»PB >*«2 ) ♦ (  ( 52LC2+PB ) **2 > )/3.  234 

PUNCH  120  235 

PUNCH  121  235 

PUNCH  103.PB.T2.T1 .T1 .TB.OL  237 

PUNCH  1 29  238 

PUNCH  104.SILCI  .SILC2.52LC1  .S21.C2  239 

PUNCH  128  240 

PUNCH  1 C4 • S 1 LL 1 • S 1 LL2 . S2LL 1  «  S2LL2  241 

PUNCH  114  2*2 

PUNCH  104.S1M1 .S1M2.S2M1 .S2H2  243 

87  PUNCH  126  24* 

PUNCH  1C4.8.P.P1  1  «  SR2 1  245 

49  PB-PT  245 

50  CONTINUE  247 

52  CONTINUE  24B 

51  CONTINUE  249 

END  250 
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H.C. ROLAND  *  NOZZLE-STRESS  MO-TI  100ATM  .7BDIA 
INPUT  DATA 


t  5000.0E-01 

20000.0E+03 

78400. 0E-05 

15000. 0E-04 

12500. CE-05 

85000. OE-03 

78000. CE-03 

65000. CE-03 

52000. 0E-03 

30000. 0E-10 

30000. CE-1 0 

30000. OE-IO 

7000.0E+03 

35000. OE+03 

22000. 0E+03 

31 00O.0E-O5 

30000. 0E+03 

47000. 0E+03 

85000. OE-Ol 

15000. CE+01 

80000. 0E-03 

15000. 0E-01 

35000. OE-01 

20000. OE-04 

30000. 0E-06 

30000.0E-05 

10000. 0t-C6 

bwOOO .OE-07 

1000C.0E-06 

LONGITUDINAL 

RIB  NOZZLE 

PB 

T  2 

T  1 

T  1 

TB 

DR 

0.000 

.150 

.010 

435.298 

2 18.333 

.089 

S1PC1 

S1PC2 

S2PC1 

S2PC2 

37273.467 

-58848.215 

-66838.347 

45263.599 

SPL  1 

SPL2 

SR  12 

SSC22 

SLS 

—5048 1 .115 

-8013.607 

-5160.472 

1 630.772 

1278.018 

SI  Ml 

SI  M2 

S2M 1 

S2M2 

35906.323 

24663.672 

27760.480 

23455.556 

LONGITUDINAL 

Rib  NOZZLE 

PB 

T  2 

T  1 

T  1 

TB 

DR 

0.000 

.150 

.015 

543.416 

21B.333 

•  1  10 

S1PC1 

S1PC2 

S2PC1 

S2PC2 

6922.805 

-44673.096 

-79588.817 

41838.526 

SPL.  1 

SPL  2 

SR  12 

SSC22 

SLS 

-65091 .374 

-1964.038 

-7829. 1 87 

3367.455 

21B0.655 

S  1  M  1 

SI  M2 

S2M  1 

S2M2 

32146.945 

1 8903.071 

33914.763 

20201 .713 

LONGITUDINAL 

RIB  NOZZLE 

PB 

T2 

T  1 

T  1 

TB 

DR 

0.000 

.150 

.020 

651 .199 

218.333 

.128 

S1PC1 

SI  PC2 

S2PC1 

S2PC2 

-13598.21 l 

-38471 .732 

-95691 .718 

43621 .775 

SPL  1 

SPL2 

SR  12 

SSC22 

SLS 

-79037.995 

4389.749 

-1 1 368.618 

8602.869 

3208.340 

SI  Ml 

SI  M2 

S2M1 

S2M2 

34060.160 

1 7701 .253 

41  044.070 

19610.878 

LONGITUDINAL 

RIB  NOZZLE 

PB 

T  2 

T  1 

T  1 

TB 

DR 

0.000 

•  150 

.025 

758.614 

2 1 8.333 

.  1  44 

S1PC1 

S1PC2 

S2PC  1 

S2PC2 

-30158.264 

-34726.437 

-1  12463.090 

4737b. 395 

SPL  1 

SPL2 

SR  12 

SSC22 

SLS 

-92445.664 

1 0936.377 

-  1 5833.24 1 

82  7b .086 

4336.720 

S 1  M  1 

SI  M2 

S2M1 

S2M2 

37965.063 

18733.969 

4B286.856 

20346.836 

LONGITUDINAL 

RIB  NOZZLE 

PB 

T  2 

T  1 

T  1 

TB 

DR 

0.000 

•  1  50 

.030 

665.640 

218.333 

.  159 

SI  PCI 

S1PC2 

S2PC1 

S2PC2 

-44708.344 

-31 756.228 

-129121 .060 

52656.492 

SPL  1 

SPL2 

SR  12 

SSC22 

SLS 

-105400.650 

1  7601 . 1 54 

-2 1219.544 

1 1 335.277 

5547.330 

SI  Ml 

SI  M2 

S2M1 

S2M2 

42616.963 

2 1 224 .229 

55422. 706 

21887.034 

LONG  I TUDI NAL 

RIB  NOZZLE 

pe 

T2 

T  1 

T  1 

TB 

DR 

0.000 

.  1  50 

.0  35 

972.263 

2 18.333 

•  1  72 

S1PC1 

S1PC2 

S2PC1 

S2PC2 
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-58061 .400 

-28939. 46C 

-145401 .560 

58380.745 

SPL1 

SPL2 

SR  12 

SSC22 

SLS 

-1 17966.190 

24331 .083 

-27498.272 

14744.236 

6625.731 

S  1  M  1 

S  1  M2 

S2M 1 

S2M2 

47553.453 

24779.245 

62382.976 

23 943.656 

longitudinal 

PB 

RIB  NOZZLE 

T  2 

T  1 

T  1 

TB 

•C  .CCC 

.  1  5C 

.040 

1076. 475 

216.333 

S1PC1 

S1PC2 

S2PC1 

S2PC2 

-70692.1 00 

-26029.401 

-161219.540 

64496.063 

SPL1 

BPL2 

SR  12 

SSC22 

SLS 

-130190.720 

JICbV. 1 30 

-34026.264 

1 8469.955 

6160.367 

S1M1 

SI  M2 

S2M 1 

S2M2 

52587.39 3 

29163.606 

691 49. 1 26 

26336.91  1 

LONGITUDINAL 

PB 

Rio  NOZZLE 

T  2 

T  1 

T  1 

IB 

O.OOG 

.  1  50 

.045 

1 1 64 .274 

216.333 

S1PC1 

S1PC2 

S2PC  1 

S2PC2 

-82769. 36C 

-22932.020 

—  1 76359 .600 

70656.231 

SPL1 

SPL2 

SR  12 

SSC22 

SLS 

-1421 12.050 

J7b40 .264 

-42563.592 

22460.049 

»54 1 .645 

SIM! 

M  M2 

S2M1 

32M2 

57636.750 

34  2  £  7  .  C  66 

7:j721  .n2  3 

28950. 004 

DR 
•  1  85 


DR 

.196 
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RIB  STRESS  PROGRAM 


This  program  is  written  In  the  Fortran  II  language  and  may  be  run 
on  any  digital  computer  designed  to  use  this  compilation  language. 

A  sketch  of  the  nozzle  liner  and  sleeve  with  circumferential  ribs 
is  shown  below.  A  typical  section  is  labelled  according  to  the  program 
naming  of  variables. 


4 


INPUT  DATA 

The  input  data,  in  order  of  read-in  are: 


First  Card 

Title  card,  which  may  contain  up  to  50  letters  and/or  numbers  to 
identify  the  case. 


Sacrnd  Cord 

This  card  is  divided  into  three  10-space  columns  and  one  3~*pace  column 
according  to  the  Fortran  FORMAT  form  (E10.I,  E10.1,  E10.1,  13).  It  contains 
in  order: 


TR  -  The  rib  axial  thickness  in  inches,  same  for  all  ribs. 

EL  -  The  liner  modulus  of  elasticity,  1 b/ in2. 

TEST  I  -  The  convergence  criterion  for  each  iteration  for  stress 
distribution.  A  reasonable  value  is  from  I  to  100.  Since 
this  number  is  compared  with  the  residual  in  the  force  balance 
equation  it  is  of  the  order  of  the  forces  involved,  which,  for 
each  section  is  several  hundred  pounds. 
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N  -  The  number  of  sections.  Including  Number  I  es  e  dummy  section 
et  the  exit  end,  Number  2  the  first  section  et  the  exit  end, 
end  Number  N  the  lest  section  et  the  entrance  end.  Sections 
Number  2  and  N  do  not  have  ribs,  end  ere  only  half  es  tong  as 
their  neighboring  sections.  Up  to  end  Including  99  sections 
may  be  used. 

Cardt  3  through  (N  +  2) 

These  cards  each  contain  In  order  the  following  values  for  sections 
I  through  N  according  to  FORMAT  (EIO.I,  EIO.I,  EIO.I,  EIO.I,  EIO.I): 

PA  -  The  airside  pressure  on  the  section  In  lb/ln2. 

PC  -  The  coolant  side  pressure  on  the  section  in  Ib/ln  . 

DL  -  The  liner  diameter  at  the  entrance  end  of  the  section  in  inches. 

CA  -  The  radial  width  of  the  rib  for  the  section.  Inches. 

OX  -  The  length  along  the  liner  of  the  section,  inches. 

Cards  (N  +  3)  through  (2N  +  2) 

These  cards  each  contain  in  order  the  following  values  for  the  sections 
I  through  N  according  to  FORMAT  (EIO.I,  EIO.I,  EIO.I,  EIO.I): 

TH  -  The  angle  between  the  liner  and  the  nozzle  axis,  taken  at  the 

middle  of  the  section  in  radians.  Note  that  for  sections  on  the 
exit  side  of  the  throat  this  angle  is  given  a  negative  value  so 
the  force  of  the  air  pressure  on  these  sections  will  be  negative. 

THCK-The  thickness  of  the  liner  at  the  entrance  end  of  the  section, 
inches. 

ER  -  The  modulus  of  elasticity  of  the  rib  in  the  section.  This  value 
is  made  0.0  for  a  section  with  no  rib.  It  will  usually  be  the 
same  for  other  sections.  Units  are  lb/ln2. 

S  -  The  initial  stress  value  inserted  for  the  section.  Any  value 

may  be  used,  but  good  original  estimating  will  shorten  convergence 
time  si ightly. 


PROGRAM  OPERATION 

After  the  data  are  entered  the  computation  will  proceed  automatically 
until  the  convergence  criterion  for  the  residuals  is  satisfied  for  all 
sections,  at  which  time  the  exit  end  liner  displacement  will  be  printed 
out.  The  computer  will  then  pause.  The  operator  may  examine  this  displace¬ 
ment.  If  it  is  positive  he  prepares  a  single  card  for  a  new  value  of  S(, 
the  exit  end  stress,  which  is  less  than  the  original  value.  The  FORMAT  is 
(EIO.I).  He  places  this  card  in  the  machine  and  starts  the  computer.  It 
reads  the  new  $j  value  and  again  calculates  the  stresses.  If  the  displace¬ 
ment  as  calculated  is  negative  the  operator  increases  S|  for  the  next 
iteration.  The  end  displacement  is  thus  made  to  approach  zero  as  closely 
as  desi red. 
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i 

*  If  a  readout  of  rib  number,  residual,  end  stress  is  desired  during 

the  relaxation  of  residuals  SENSE  SWITCH  3  Is  turned  on.  When  the  end 
.  displacement  Is  close  enough  to  zero  the  operator  turns  on  SENSE  SWITCH 

2.  This  causes  a  print  out  of  rib  number,  residual,  stress,  section  dis¬ 
placement,  and  force  (in  pounds)  on  the  rib. 

The  FORTRAN  program,  entitled  H.  C.  Roland,  RIB  STRESS  PROGRAM,  follows 
on  pages  124  through  126.  Sample  data  are  given  on  page  127,  and  the  final 
print  out  on  page  128. 
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C 

c 

c 

c 

c 

c 


C  H»  C.  ROLAND  **  RIB  STRESS  PROGRAM 

SENSE  SWITCH  3  ON  GlVuS  READ  OUT  OF  RIB  NO. tftES 1  DUAL .  AND  STRtSS 
INPUTS  *  TR-RIb  THCK  4  TL-LINR  THCK,  EL-LlNER  MOD  '..F  FLAJ-.T. 

TEST  1-1  ST  CCNV  CR1T»  TEE  T2-2ND  CON'/  CR1T,  N-NO  OF  UTRE.SSES 
PA(I)-AIR  PRESS..  PC(I)-COOL.  PR=-SS..  LL<I)-DIA.  lin-r. 

CA  (  1  J-COOL  .  ANNULUS.  DX(1)-SECT.  LC*TH  ON  L  I  NtR  .  c(I)-cTRESS 
TH  (  I  )  —SECT  .  ANG..  TUCK  (  1  > -L  I  NER  THCK.  ER<!)-RI->  EL  Ai-T  MOO. 

GIVE  TH  NuG.  ON  OUTLET  SIDE  OF  THROAT 
GIVE  ER  <11=0.  FOR  SECTS.  U.2.  N 
ICO  FORMAT  (50H 

101  FORMAT  (1  OH  INPUT  DATA) 

102  FORMAT (E10.1.E1U.1,:  10.1.1J) 

10  3  F  ORM AT (ElC.l.ElC.l.plO.l.ElC.l.tlO.l) 

1  C  4  F  OR  Mi  AT  ( E  l  0 . 1  ,  t:  1  0  •  1  , ..  1  C  .  1  «  E  1  0 . 1  ) 

1  C  5  FORMAT  < 1 3H I TERAT I  ON  NO.) 

1C  6  FORMAT  (8HL00M  NO.) 

107  FORMAT  (16) 

108  FORMAT  (7HR1H  NO  •  •  4X  .  8HRtS  1  DUAL  .  OX  .  6HSTRES  ?  .  7>  .  vi-L'l  LT  A  , 

1  7X.5HFCRCE  ) 

1~9  FORMAT  (  14  .  •>X.-  1  2. 1  »T  12  .  1  .F  1  2.8  ,F  1  £.  1  ) 

1  1  0  FORMAT  (  1  7hpo.jv.LL  ■'  COVRLE Tt. _  ) 

111  FORMAT  <  2v . 1 h 1 . 1 Dx . L 1 o • I  ) 

112  FORMAT  U7HS.NJ  ;-! -PwAC'_MENT-.Fl2.e) 

01  MENS  1  CN  PA  (20  )  ,HC  (2C  )  »L>L  (  :.C  )  .CA  (?c  )  ,;;x  (2C.  )  .Th(£0  )  .Thun  (20  )  « 

1  ER  (20  )  .  (2*.  )  ♦U'M(2'.  >♦«<-.  )  .ait/  )  .  Cu:.R  (  2C  )  .  /  ( 2  0  )  .  m  (  2C  )  .  U  (  2  0  )  . 

2Y(20  )  .  SUV  (20  )  .  •iur-'-l  <21  )  .i.'t-L  (20  )  .f  <  '  r  > 


LF  =  C 


RE AC  ICC 
PUNCH  ICO 
PUNCH  1  •'  1 

READ  1  C  2 .  TR,  --L.TEST1,  N 
PUNCH  IT.  TR.  EL.  TEST  1  .  N 
SO  E  I  =  1  .(» 

READ  1  C  3  .  PA(i).  DC  (  |  )  ,  DL  (  1  )  «  0**(I).  UX  (  1  ) 
s  punch  i  c  3 »  faid.  *  .in.  dl(I>.  a  <  i  i ,  .  x  <  i  > 
DO  1C  I  X  -  1  , N 

RE  A  j  1  C  «  .  TH(lx),  TtiC*(!x).  -  R  (  I  »  )  .  L  t  I  X  ) 
ir  punch  1L4.  Th(1x),  thlmixi.  •  f  t  1  *  )  •  ;>  (  I  x  ) 


NO  =  N— 1 
DO  IE  J- 1 .N 

IF  < J-l  )9,9,  1  1 

9  DM < J ) =3. 142* (DL ( J > +CA ( J ) ) 

A ( J )=3. 142*DX ( J ) *OL ( J ) 

GO  TO  3 

11  DM  (  J  )  =  1  • 5  7* ( UL ( J ) +  UL (J-l  )  +2 • *  C A  (J)  ) 
A ( J )  =  1  .571 *DX ( J ) * (DL ( J ) +DL (J-l  )  ) 

3  IF(J-N>6.7,7 

6  COS ( J > — COSF (  (Tn  1  J l+TH ( J+l  )  )/2. ) 

GC  TO  8 

7  COS ( J ) =COSF ( TH ( J ) ) 

8  COSR ( J ) =EL*COSF ( TH ( J  )  ) 

W ( J ) =3 . 1 42*DL ( J ) 

C ( J ) = ( PA ( J ) -PC  <JI)*A(J)»SINF(TH(J)) 
D  (  J  )  =Vv  (  J  )  «  THCK  (  J  )  *COS  (  J  ) 

Y ( J ) =  DM ( J  >  *  ( TR*  *  3  >  *ER ( J  > / ( C  A ( J ) *  *3  > 
15  CONTINUE 
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L?=I.P+  1 
I  T  -  I'! 

CONTINUE 
Go  1C  r;Astr,N 
—  ^  u  -  o  • 

IT  <NA-\>  45.47416 

X  2."  u=NA.NC 

ADU  =  .  L»  (  5  (  L  )  +5  C_+  1  )  )  *CX (L+ 1  ) XCOSR  CL+ 1  ) 

SUN  (  \ A  )  =  AiJD 
CCNT I \Jr 
00  55  NF=2.N 

SUM  1  <  6K  )  s[iX  (MF  )  *.25“  (  S  <  Mr  )+?.  <  AF—  1  )  )/COiP  (  MF  ) 

CC\T  I  NUr. 

DO  2  Ms?.'; 

5  !  M  >  sy  ( |/.  )*  (  SUM  1  (  W  )  +5U1'1  CM)  )+u  C  -i  )*5  (  M  >-D  (M-l  )  *£  ( >*»—  1  )-C  C  *>'•  > 
1 T- i T+l 

IF  (SENSE  SWITCH  3)  26.27 
PUNCH  1 C  6 
punch  k?,  lp 

PUNCH  1 05 
PUNCH  107, IT 

p'jnch  i  rn 
punch  in,  SCI) 

r\n  Tc.  ..*■->  =  7  ^  fvj 

PUNCH  1  O,  MP  .  U  (  MP  )  ,  4  (  vf  ) 

pji  *  f  -z  c  ,  n 


1  r  (  fihi.h  c 

(  MC  )  )  -T-:  ST  1  )  26 . 2rt  .  2  ■> 

26 

CCVT  1  NU'; 

CO  TO  -f 

20 

MX  •=  if 

DO  3'' 

*  ♦  ‘4 

1 F  ( Ah  i  i 

••  (  -  1  ( tv  t  M  y  1  )  )  . 

.v'  ♦  3  £ 

32 

*» 

X 

11 

N 

•> 

3  " 

con r  i  vui 

s  CM>  |sU  ■. 

•  )  -6  (  .1 V  )/(,'•»  (  Y  (  ■'  Y 

...  y  t r.;  v 

GO  TC  1  < 

CO  NT  |  NU- 

DC  N4  MY-O.N 

DHL  (MY  )  =i,UM  1  .  My  )  *MP  ({.•'>  ) 

r  (  MY  )  =  Y  (MY  )  *  L.'Nl.  (MY  ) 

CCNT I NUl 

0 T  =  05 L  (  2  )  *  •?.'■>*  (  S  (2  )  +  S  (  I  )  )  *ux  (  2  )  X  -O^U  (  2  ) 

PUNCH  112.  I  IT 

PAUfF 

IP  (  r,  [.  N  0,  y  0  .v  1  TCh  2  )  o  4  .  e>  b 

SENSE  SWITCH  2  ON  I F  PPCoLL  M  C  ‘L  E  T  -  .  Pk'fcSS  START 

SWITCH  2  OFF --CHANG!.  SCI)  DPI-'.  SIGN  TO  END  UlS°.b  PRESS  START 

GO  TO  66 

READ  1  <"  4  ,  SCI) 

GO  TO  16 
PUNCH  1 1 0 
PUNCH  1 Cfc 
PUNCH  1 C  7 .  LP 
PUNCH  1 06 
PUNCH  107.!“ 

PUN oh  1  OH 
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PUNCH  1114  S  <  1  ) 

00  6C  MD*2.N 

6C  PUNCH  lC9tMO.R(f>T).Si(«D)iDiL(''li:iltF(MO) 
END 


H  C  ROLAND  Rib  STRESS  ** 
INPUT  DATA 

250.0E-04  450.0E+05  400 
500.0E-02  600.0E-01  103 

100.0E-00  600.0E-01  134 

200 • OE— OC  7CC.0E-01  116 

300.0E-OC  700  »0E  — C 1  1CB 
400 •OE— 00  80O.CE-01  100 

500 • OE— 00  80C.OF-01  920 
6C0.0E-00  BOC.OE-Ol  B40 
700 • OE-CO  9CC.OE-01  760 
800 • 0E-00  9CG.ut-01  bbC 
900.0E-G0  lOG.ut-00  toOO 
1  00 .OE+O 1  1CO.GE-OC  520 
UO.OE+Cl  IIO.CE-CC  440 
1  20  •  OE  +  O 1  110* OE— 00  360 

150.0E+C1  12C .OE-CO  280 

140. OE+O 1  120.0E-00  240 

280 ■ OE+O 1  140.0E-00  600 

280 .OE+O 1  1 40 .CE— 00  104 

280 .OE  +  O 1  14  0. OE  —  00  141 

280  .OE  +  O  1  14  0j.CE  — 00  150 

280 . OE+C 1  14(.0t-00  150 

-140.0E-C3  34C.GE-G4  COO 
-140.0E-03  3JO.OL-04  0 OC 

—  140. OE  — C 3  320.JE-04  460 

—  140. OE  —  03  310.CL-04  450 
-140.0E-C3  3CC.0L-04  460 

—  140.0E  —  0  3  290.Ut.-04  46  c 

—  14  0.  OE  — 03  280.01.-04  45T, 
-140.0E-03  28n .OE  —  04  464 

—  140. OE  — 03  270. CE  —  04  450, 
-140.0E-03  260.0E-04  450, 

—  140.0E  — 03  250 • CE  —  04  450, 

—  140.0E  —  03  24G.0E— 04  450, 

—  140.0E  —  03  230 .Ot  —  04  450, 
-140.0E-03  220.0E-04  450, 

500 . OE— 03  250.0E— 04  460, 
845.0E— 03  280.0E-04  450, 
845 . OE  — 0  3  310.0E-04  450, 
670  «  OE  — C  3  330.0E-C4  450, 
000 • OE— 99  340.0E-04  450, 
COO • OE— 99  35C.0F-04  OGO, 


MO— T I  200  ATM  .240 1 A 
,0E-0 1  20 

1 CE-02  100.0E-O3  130.0E-04 
i0E-02  100.0E-O3  126.0E-03 
i0E-02  1 00 . OE-03  252.0E-03 
iOE-02  10U. 0E-03  252. 0E-03 
>0E-02  100. 0E-03  252. 0E-03 
,0E-03  100. 0E-03  252. 0E-03 
.OE-03  100.0E-O3  252.0E-O3 
, 0E-03  100.0L-03  252.0E-03 
. OE-03  100.0t-u3  252. OE-03 
iOt-03  lUO.OE-oj  252.0E-03 
, OE-03  1C0.CL-U3  252. OE-03 
, 0 E  —  C 3  100.CL-03  252. OE-03 
.0E-03  100. OE— 03  252 .OE-03 
,OE-C  3  ICC. OE-03  252.0E-C3 
, OE-03  100. OE-03  252. OE-03 
,CE-C3  ICO. OE-03  350.0E-C4 
,  CE  —  02  ICO. OE-03  350.0E-04 
, CL  — 02  loC.OE-i-3  290. OE-03 
■0E-02  100. OE-03  250. OE-03 
,CE— Oc  100. OE-03  125.0L-03 
uu -99-200 .Oc  +  O 1 
Ot -99- I vU.Ci+U 1 
ot+06-120. 06+01 
i  Ot  +06— 650 . Ot  —  00 
,Ot  +  C6-50'J.nt-00 
fit  +  06— 490  »  Ot  — 00 
Cl +C6-5CO.CE-00 
cl +05-550. OE-CO 
OE+05-3CG.OE-CO 
Ct+05— 440.CE— oO 
Ot+05— 400.C2— 00 
06+05-300. Ct+01 
0L  +  05-4  70  »0t.  +  o  1 
ot_  +  05— 750.UL  +  O1 
0t+05— 89C.0t+01 
Ot  +  06— 4 10. Ot  +0  1 
Ot  +  C.6-21  G.Ot  +  0  1 
Ot  +  05  1  13. Ot 
ot+o-o  bJo.cioi^i 

Ot— 99  700.0L+0 1 
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PROBLEM  COMPLETED 
LOOP  NO* 

9  • 

ITERATION  NO* 

1 


RIB  NO* 

residual 

STRESS 

DELTA 

FORCE 

I 

200.0 

2 

39*8 

509.8 

.00000298 

0.0 

3 

35*0 

696.9 

.00000078 

2.2 

4 

32*6 

894  .6 

-.00000317 

-8.5 

5 

36*4 

1221.7 

—  »  0000oB4 1 

-21.1 

6 

-21.1 

99  1  .8 

-.00001453 

-34.0 

7 

-29*0 

626.7 

-.00001995 

-43.1 

B 

-35.0 

44.5 

-.00002319 

-46.0 

9 

-6.6 

-300.0 

-.00002378 

-43.0 

10 

15.7 

-440.0 

-.00002237 

-36.5 

1  1 

31.5 

-400.0 

-.00002014 

-29.3 

12 

-33.7 

-2541 .1 

-.00001479 

-18.9 

1  3 

10.1 

-4700.0 

-.00000039 

-.4 

14 

32.3 

-8218.9 

.00002810 

26.0 

15 

-29.6 

-7942.0 

.00007215 

57.3 

16 

-1  .9 

-4100.0 

.000 1 u 1 4b 

1  16.4 

1  7 

26.6 

-2100.0 

.000 1 9680 

216.9 

18 

24.1 

10472.8 

.00009141 

267.4 

19 

31  .0 

8182.9 

.00004828 

165.7 

20 

-3.1 

7930.0 

.00001  1  18 

0.0 
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APPENDIX  B 

FORTRAN  PROGRAMS 
FOR 

BACKSIDE  COOLING  CALCULATIONS 
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LINER  BACKSIDE  COOLINC  COMPUTER  PROCRAM 


The  program  described  below  is  written  In  Fortran  II  language.  It  is 
designed  to  calculate  the  bulk  temperature,  static  pressure,  convection 
heat  transfer  coefficient,  surface  temperature,  boiling  heat  transfer  coef¬ 
ficient  and  boiling  burnout  heat  load  along  a  nozzle  liner  coolant  annulus 
with  longitudinal  ribs.  The  program  propar  is  preceded  by  a  brief  dis¬ 
cussion  of  the  equations  used  for  calculation,  in  this  discussion  the 
notation  used  is  the  same  as  that  used  in  the  computer  program. 


COOLANT  ANNULUS  DESCRIPTION 

The  coolant  flows  axially  in  an  annulus  around  the  nozzle  liner.  The 
coolant  side  surface  of  the  liner  is  a  surface  of  revolution  about  the 
nozzle  centerline  which  Is  designated  the  x-axis.  For  calculation  of  the 
required  variebles  in  a  stepwise  fashion,  the  liner  is  divided  axially  into 
small  finite  sections.  A  description  of  each  section  as  to  length,  radius, 
angle  of  the  liner  surface  with  the  x-axis  and  distance  across  the  coolant 
annulus  perpendicular  to  the  liner  surface,  provides  the  necessary  liner 
and  annulus  description.  Using  these  variables,  the  liner  heated  area, 
coolant  annulus  flow  area  and  flow  path  length  along  the  liner,  are  calcu¬ 
lated. 

A  cross-sectional  view  of  a  typical  liner  and  annulus  section  is  shown 
below  with  the  necessary  variables  identified. 
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If  the  section  length  Dxn  is  made  small  compared  to  the  radius  of 
curvature  of  the  liner,  the  heated  area  of  the  section  is 


A  HO  Tn  m  2  w  Dx« 

Cms  (THm) 


(  I  ) 


The  coolant  annulus  flow  area  is 


ACcci^  m  Tr(  *  *2„)  P2n 


(  2  ) 


Now 


*  <€7„  *■  D2„  Co*(Ttf„) 


(  3  ) 


So 


*  zrrz!„  D2n+  rrc»i(T«„)  dz„ 


(  4  ) 


The  above  equation  for  the  coolant  flow  area  applies  only  in  sections 
where  there  are  no  ribs.  If  there  are  RN  ribs  of  thickness  RT,  the 
cool  ant  flow  are  i s 


fiCOOLn  «  MB/*  §2„  *  ^rce>5(T^)DZf>  -  AEV  KT  02 r 


(  5  ) 


If  a  constant  flow  area  is  desired  let  AC00Ln  »  OCOOL  »  constant. 


Then 


*  . 

7TCas(TMn)  DZ„  +(a7re/n.  gA/  ftr)  D2n  -  DCcot  -  o  (  6  ) 


and  D2, 


,  -(znttt*.  zn  KT)+lz'r*f*-**KT)++rrco*(TH»)occoc^’  ^  y 


2TTCCi(-T»„) 
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BULK  TEMPERATURE,  VELOCITY  AND  PRESSURE  CALCULATIONS 

Given  a  coolant  water  flow  rate  and  a  heat  load  distribution  along  the 
liner,  the  bulk  temperature  at  the  nth  section  Is  calculated  from 


2'CP-  PAK  Z-Cp.  PLtZ 


where 


a  bulk  temperature  at  section  n,  *F 
dp*  ■  heat  load  on  liner  at  section  n,  B/hr  ft2 
b  liner  section  area  at  section  n,  ft^ 

CP  b  coolant  heat  capacity,  B/lbm  * F 
Plk  b  coolant  flow  rate,  Ibm/hr. 


The  velocity  at  each  section  is  calculated  from 


V*  -  - -  (  9  ) 


ACpccm  HHO* 


where  V„  a 

m 

velocity,  ft/sec 

*nc„  = 

density,  lbm/ft* 

M  a 

coolant  flow  rate, 

Ibm/sec 

a 

coolant  flow  area, 

ft2. 

The  static  pressure  is  calculated  from 


>KSSS  8  "PBess^  -  OPS  -  DPSn-  "Dpv  . 

*  r»-i  n-i  n  ^ 


(  10  ) 


where  PR*ssx  » 

B 

b 

OWa  » 


OPy 


static  pressure  at  section  n,  psia 
static  pressure  at  section  n-l,  psia 
frictional  pressure  loss  in  exit  half  of 
section  (n-l) ,  psia 

frictional  pressure  loss  In  entrance  half 
of  section  n,  psia 

pressure  loss  due  to  acceleration  from  (n-l) 
to  n,  psia. 


The  pressure  loss  terms  In  equation  (  10  )  are  calculated  as  follows: 


For  DPS 


9./SV/©  Or„ 


,  N  -o.zs 


(  II  ) 
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where  ^PSn  „ 
DXn  . 
DM*  » 

«n  - 


frictional  pressure  drop  in  half  of  section  n,  psia 
section  length,  inches 
coolant  annulus  hydraulic  diameter,  ft 
Reynolds  number  at  section  n. 


For  DPV 


Dpv  /  gao»  v,*~ 

(n-0-*n  ^ 


(  12  ) 


whe  re  DPV  = 
?HOn  = 
^  = 
3«l  = 


acceleration  pressure  drop  from  (n-l)  to  n,  psia 
coolant  density  at  section  n,  lbm/ft' 
coolant  velocity  at  section  n,  ft/sec 
32.2  lbm  ft 

Ibf  sec2 


FORCED  CONVECTION  AND  SURFACE  BOILING  HEAT  TRANSFER 

The  forced  convection  heat  transfer  coefficient  is  calculated  from 

-0.2 

BARNn  .  S/34*Vn  *£„  P*n  (  ,3  ) 


where  PRn  =  Prandtl's  number  at  section  n. 

For  evaluating  REn  the  viscosity  is  calculated  from  the  following 
empirical  equation 


FMU„*  -Qigl£  _  ^4-S*  i6R/)TCn6-(m»")  TC^~  (  14  ) 


where  Fr*iUn  = 


viscosity  at  section  n,  Ibm/ft  sec 
film  temperature  at  section  n,  ®F 


For  calculating  PRn  the  following  empirical  equation  is  used 


i£o  _  0.002-m^ 


x  10 


)TC 


3 

n 


(  15  ) 


Equations  (  14  )  and  (  15  )  are  plotted  with  experimental  values  for  com 
pari  son  in  Figures  I  and  2  respectively. 
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Temperature,  F 
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Fig.  45  Prondtl  Number  i or  Water 
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The  liner  coolant  side  surface  temperature  in  forced  convection  is 
calculated  from 


baith* 


(  16  ) 


where  T$n  - 
T«*  - 
Q  n  * 


surface  temperature,  *F 
bulk  coolant  temperature,  *F 
heat  load  at  section  n,  B/hr  ft2 
heat  transfer  coefficient. 


Since  REn  and  PRn  are  functions  of  the  film  temperature  TCn,  and  TCn 
varies  with  TSn,  equations  (  13  ),  (  14  ) ,  (  15  ),  and  (  16  )  are  used  to 
solve  for  TSn  by  successive  approximation. 


At  each  section  the  saturation  temperature  is  first  calculated  based 
on  the  static  pressure.  For  pressures  to  400  psia  the  saturation  tempera¬ 
ture  is  calculated  from 


_  .3 _  *• 

TSAT.  x  240  +TfcCSS- «’?««»»  PRESS  - 

n  R  “ 


(mUO 


Z7S+-PKKS1, 


(  17  ) 


For  pressures  from  400  psia  to  1000  psia  the  following  equation  is  used 
TSAt,  *  390  +  0.16,  "PREtt*  (  18  ) 


For  pressures  above  1000  psia  the  saturation  temperature  is  calculated 
from 


"T5atk  >  4la  ■+  o,o8z.TitMs„ 


(  19  ) 


Equations  (  17  ),  (  18  ),  and  (  19  )  are  plotted  with  experimental  values 
of  the  saturation  temperature  in  Figure  3* 

If  the  required  surface  temperature  exceeds  the  saturation  tempera¬ 
ture  at  any  section  by  more  than  !0*F  a  boiling  heat  transfer  calculation 
is  made.  The  convection  heat  load  is  first  calculated  by  assuming  a  value 


13K 
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f°r  *n^  using 

*’*'*•  (TS‘>-'rm-)  (20) 

The  boiling  heat  transfer  load  Is  then  calculated  from 

1.S4 

»  0.2[TSn-(TiATn>,o)]  (  2)  ) 


(Qn) convection  +  (On) boiling  4  Q  another  value  of  TSn  Is  chosen 
and  the  process  repeated  until  the  proper  value  of  TSn  is  reached. 

The  burn-out  heat  flux  is  evaluated  using  each  of  the  following  three 
equations 


PBow.(4xiof^i» .Ota  qtso&.Jh- &oftaijaL»y  w (eaag») 


(  22  ) 


Q®2-n  *  7ooo(crr4u%„)  V* 

-w  -'S'V* 

»  (4*10  ^  I  i-.OZ  \/„ 


(  23  ) 

(  24  ) 


I 3TSOBn  t  “I" SAT^--|-I5n 

At  sections  where  there  is  no  boiling,  the  burn-out  heat  flux  is  set  equal 
to  0.0. 


PROGRAM  SWITCH  SETTINGS 

If  the  program  is  executed  with  Sense  Switch  I  off,  the  pressure  dis¬ 
tribution  is  first  calculated,  then  the  inlet  pressure  is  automatically  set 
to  give  a  static  pressure  of  10  psia  at  the  last  section.  This  automatic 
inlet  pressure  adjustment  does  not  occur  if  the  program  is  run  with  Sense 
Switch  I  on. 
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If  the  input  data  It  read  in  with  Sense  Switch  3  off  the  coolant 
annulus  widths  input  as  data  are  used  for  ail  calculations.  If  Sense 
Switch  3  is  on  when  the  data  is  entered,  the  coolant  annulus  widths  are 
automatical ly  calculated  to  give  a  flow  area  equal  to  DC00I  at  each  section 


INPUT  DATA  CARDS 

First  data  card  -  FORMAT  (50H  -  50  spaces  -  ) 

This  card  provides  for  input  and  output  of  up  to  50  alphabetic  charac¬ 
ters  to  provide  a  title  for  each  case. 

Second  data  card  -  FORMAT  (15,  FI0.5,  FI0.5,  FI0.3,  F10.5,  FI0.5) 

This  card  contains  data  in  order  as  follows: 


N  -  Number  of  liner  sections. 

FLR  -  Coolant  flow  rate,  Ibm/sec. 

TC00L-  Inlet  coolant  bulk  temperature,  *F. 

PRIN  -  Inlet  static  pressure  (plenum),  psia. 

AIN  -  Inlet  area  at  the  bend,  in. 

DH IN  -  Inlet  hydraulic  diameter  at  the  bend,  inches. 

Third  data  card  -  FORMAT  (15,  5X,  F10.5,  F10.5.  FI0.5) 

This  card  contains  data  in  order  as  follows: 


NRS  -  Number  of  sections  from  entrance  having  ribs. 

RN  -  Number  of  longitudinal  ribs  around  liner. 

RT  -  Rib  circumferential  thickness,  inches. 

OCOOL-  Coolant  flow  area  for  automatic  annulus  width 
calculations,  in2. 

Fourth  card  through  (3  +  n)th  card  -  FORMAT  (F10.5,  F10.5,  F10.5,  FI0.5, 

F10.0) 

These  n  cards  contain  the  following  data  in  order  for  the  n  liner 
sections: 

Fin  *  Liner  outside  radius  at  n1^1  section,  inches. 

D2n  *  Coolant  annulus  width  at  nth  section,  inches. 

THn  =  Liner  angle  with  axis  at  nth  section,  radians. 

DXn  »  Length  of  nth  section,  inches. 

Q.n  *  Liner  heat  load  at  nth  section,  B/hr  ft2. 


All  input  is  punched  out  as  a  check  on  the  input  data. 
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OUTPUT  NOTATION 

Output  values  ara  (abated  according  to  the  following  notation: 

N  -  Same  as  input. 

FL.RATE  -Coolant  flow  rata,  Ibm/sec. 


•F. 


ft2  *F. 

B/hr  ft]?. 

B/hr  ft2. 

B/hr  ft2. 

Note:  If  QB01  *  QB02  a  QB03  a  0.0  at  a  section,  it  shows  there  Is  no 
boiling  at  the  section. 

The  program  follows  on  pages  141,  142,  and  143,  and  sample  data  and 
output  on  pages  144  through  155. 


TC00L  -  Same  as  Input. 

PRIN  -  Same  as  Input  (may  be  different  value). 

AIN  -  Same  as  input. 

OH  IN  -  Same  as  Input. 

SEC.  -  Section  number 

T-BULK  -  Bulk  coolant  temperature,  *F. 

T-SURF  -  Liner  coolant  side  surface  temperature, 
T-SAT  -  Coolant  saturation  temperature,  *F. 

VEL.  -  Coolant  mean  velocity,  ft/sec. 

PRESS  -  Coolant  static  pressure,  psia. 

Q  -  Liner  heat  load,  B/hr  ft2. 

H  -  Surface  heat  transfer  coefficient,  B/hr 

QBOI  -  Burnout  heat  flux  from  equation  (  22  ) , 

Q.B02  -  Burnout  heat  flux  from  equation  (  23  ) , 

QB03  -  Burnout  heat  flux  from  equation  (  24  ) , 
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C  C  H  C  ROLAND  **  BACKSIDE  COOLING  PROGRAM 
C  SENSE  SWITCH  1  ON  BYPASSES  AuTO  PRESS  SETTING 

C  SENSE  SWITCH  3  ON  CAUSES  tUuALlZATION  OF  ALL  FLOW  AREAS 

DIMENSION  R1  <50)»02<50)»TH(50>,(j(50>  .COS ( 50  > • ACOOL (50  >  «  OH (50  > » 

1 AHOT (50 ) .G  <50 ) . V (50  > . DTB (50 ) <  TB ( 50  I »  TF ( 50  > . FMU < 50 ) .PR  < 50  )  • 

2BARH ( 50 ) .DTSI50  t.TS (5C ) .RE (50 ) .DPS (50  > .DPV(50 ) . VH (50) .RHO (50  )  . 
3PRESS(50) .TSAT (50) . DX (50 ,TSB<50 ) . TC (50 > . QBO < 50 > . QB2 ( 50 ) .UB3(50  > 

100  FORMAT  <  50H  ’ 

101  FORMAT  ( 1  OH  1 NPUT  DATA) 

102  FORMAT  <4X« 1 HN . 3X . 7HFL • R ATE , 3X . 5HT COOL , 6X » 4HPR I N . 
1BX.3HAIN.6X.4HDH1N) 

103  FORMAT  (  I  5 . F 1 0 • 5 .F 1 0 « 5 . F 1 0 « 3 . F 1 0 .5 • F 1 0 #5 ) 

104  FORMAT  (2X.3HSEC.6X.2HR1 .BX.2HD2.faX.2HTH.BX.2HDX.6X. 1HQ) 

105  FORMAT  <F 1 0.5. F IC.5.F 10.5.F1 0 .5.F 10.0 ) 

106  FORMAT  < I5.F10.S.F1C.5.F10.5.F10.5.F10.0) 

107  FORMAT  (6H0UTPUT) 

108  FORMAT  (  3HSEC  «  3X  .6HT-BULX  «4X  .6HT-SURF  «5X.  5HT-SAT  .6X.4HV=.L»  . 

1 5X . 5HPRESS ) 

109  FORMAT  <  I  2 . F 1 0 . 1 «F 1 0 . 1  . F 1 0 •  1  .  F 1 O  •  1  • F  1  0 • 1  > 

110  FORMAT  (  I  2 . F 1 0 • 0 .  F10.1.  F12.0.  F12.0.F12.0) 

111  FORMAT  (3HSEC.4X. 1HU.9X. 1 HH . 9X« 4H0B0 1 . BX . 4H0B02 . fax « 4H0803  ) 

112  FORMAT  (8HRI8.SECS.3X.BHN0*  R 1 BS • 2X . 1 OHR I B  THCKNS . 4X . 5HDC00L ) 

113  FORMAT  ( 1 5.5X.F 1 0.5.F I  0 »5.F 1 0 .3 > 

80  DO  40  L= 1 .20 

READ  100 
PUNCH  100 
PUNCH  101 
PUNCH  102 

READ  103.  N.  FLR.  TCOOL .  PR  I N .  AIN.  DH I N 
READ  113.  NRS.  RN.  RT.DCCOL 
PUNCH  1 03.N.FLR.TCOOL.PRIN. AIN.DhIN 
PUNCH  1 12 

PUNCH  113.  NRS.  RN.  RT.DCOOL 
PUNCH  1 04 

LV  =  0  \ 

DO  20  I  X* 1  «N 

READ  105.  Rl(IX).  D2(IX).  TH(IX).  DX (  I X )  .  Q(IX) 

C0S<  I  X  )  -COSF  (  TH  (  I,X)  ) 

42  IF  (SENSE  SWITCH  3)  45.20 
45  A=3. 1 4*C0S (IX) 

B  =  6  »  28*R 1  (  IX)-RN*RT 
B  1  -6  »  2B*R 1  (IX) 

C= -DCOOL 

IF  ( I X-NRS ) 180, 180. IB! 

180  D2 (  I  X )  =  (  -a  +  SQRTF <<B**2>-4.*A*C))/<2.*A) 

GO  TO  20 

181  D2 (  IX  )  =  (-Ol+SURTF <<dl**2)-4.*A*C)>/<2.*A> 

20  PUNCH  106,  IX,  R 1  (IX),  D2 (IX).  Tn(IX).  DX(1X),  C(  I X ) 

8  VIN= 144 .*FLR/ ( A IN*62.4 ) 

VHI  N  =  62.4*  (  VIN**c;  )  /  64  *  4 

V0= I 44.*FLR/ (62.4*  < DCOGL+RN*RT*D2 ( 1  )  >  > 

VH0*62.4* ( VO* *2 )/64 .4 
66  DO  7  I  =  1  «N 

IF  (I-NRS)  47.47.150 

4  7  ACOOL (  I  ) *D2 (  I  )*<6.28*R1  ( I  )*3.14*C0S<  I  >*D2( I  ) — RN*RT ) /  1  44  » 

DH  ( I  )  s24.* ACOOL ( I  >/ (6»28*R1  < I  >  +  3. 1 4* COS  < I  ) *D2 ( I  ) +RN* ( D2 ( I  )-RT  >  ) 
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GO  TO  1  ‘31 

150  ACOOL  < I  >»D2< I l*<6.2B*Rl  < 1  1 +3. 1 4*C0S < I  >*D2( 1  >1/144. 

151  DH ( 1  )*48.*AC00L< 1  > / C 12.56*R1 ( 1  1+6.28*D2< I 1*C0S< I  >  > 

AHOT ( I )«6»283*DX< I >*Rl ( 1 l/COS! 1 1 

G< I  l*FLR/ACOOL< 1  ) 

0TB ( I  1*0  < 1  )*AHOT ( I 1/ <7200 .*FLR* 144. ) 

IF  (1-1)  1,1,2 

1  TB  <  I  >*TCOOL+DTB  <  1  ) 

GO  TO  3 

2  TB< I  >*TB< I-l  >+OTB< 1-1  >+DTB< 1  ) 

3  RHO  < I  ) =62.6—  < 1 .5E-03 ) *  TB < I  1 - < 5 .3E-05 1  * <  TB(I>**2> 

V  < I  >«G< I  )  /RHO  < I  ) 

VH<  1  )=RHO< I  >  *  < V  <  1  1**2  1/64.4 

FMU  < I  1 *.045/TB< I  )- (4.5E-08  >*TB  < 1  1+ (2.E-1 1  I* (TB  <  I  1**2 ) 

50  RE  1 N*62 .4*V  1 N*DH 1 N/ < 12.*FMU< 1  > ) 

DPI NS=. 1 *VHI N/<RE IN** .25 1 
DP1NV'-VH0/144. 

RE  <  1  )*RHO  < I  ) *V  < 1  1  *DH <  I  )/FMU<  I  1 

DPS  <11  =  <  9. 15E-05  >*VH <  I  )*DX  < 1  )  / < DH ( I  )*COS< 1  )*<RE ( 1  )**.25)  1 
IF  < I-l  )  5,5,6 

5  DPV < I  1  =  < VH<  I  J-VH01/144. 

'PRESS!  I  ) *PR 1 N-DPS ( 1  1 -DPV  < I  ) -DP  1 NS-DP I NV 
GO  TO  7 

6  DPV < I  1= < VH(  I  )-VH <  l-l  )  )/144. 

PRESS ( 1  IMPRESS! l-l l-DPS! 1  >-DPV< I  >-DPS< I-l  1 

7  CON V 1 NUE 
LV*LV+1 

IF  (LV-l  )  81.81  .65 

ei  if  (sense  switch  n  65.82 

82  PRIN=PR IN-PRESS (N 1+10. 

GO  TO  8 

65  PUNCH  107 

PUNCH  1 02 

PUNCH  103.  N.  FLR,  TCOOL,  PRIN,  AIN,  DH I N 
DO  25  K= 1 ,N 

IF  (PRESS (K 1-400. 1  70,70.71 

70  TSAT  (K  1  =240. -.PRESS  (X  )  -  .  00  1  250*  ( PRESS  ( K  I  **2  1  -660  •/ <  2  •  75+PRESS  ( K  1  1 
GO  TO  79 

71  IF  <PRESS<K 1-1000. 1  72,72.73 

72  TSATIX 1=390. +. 16*PRE3S(K 1 
GO  TO  79 

73  TSATIK  1=470.  +  .082*PRESS < K  1 

79  CONTINUE 

LOP*  1 
LOT  *2 
LZ  =  0 

TSB (K 1 =TSAT ( K 1+ 1 0. 

93  TC (K  )* < TSB(K  1  +  TB (K 1  1/2. 

FMU (K 1  =  .045/TC (K 1 - ( 4 . 5L-08 ) * TC  <  K )  +  <2.t-| 1  1* (TC  <K 1**2 1 
PR <K )=4B0./TC<K >- <2 .E-03 >*TC <K )♦ (7.E-09 >* < TC(K 1**3 > 

BARH (K 1“5134.*V  <K 1* (  < FMU < K 1 / ( DH ( K 1 *RHO ( K 1  * V ( K I  1)**.2)/<PR(<> 

1 ** .667 1 

QC*BARH (K 1* ( TSB (K l-TB (X 1  1 
LZ=LZ+1 

IF  (OC-Q(Kl)  30, 95,4b 
48  IF  <  LOT— LOP  1  49,95,49 
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49  TSB  <K  )*TSB(O-20. 

GO  TO  93 

30  IF  (LZ-1 )51 ,51 .27 
27  TSB(K)=TSB(K>+5. 

LOP*LOT 
GO  TO  93 

95  BARH (K >=Q (K>/ < TSB (K  )-TB <K> > 

TS(K  )  =  TSB(0 
QBO ( K  )  *0. 

QB2 (K  )  =0. 

QB3 (0=0. 

GO  TO  25 

Ol  T  SB  ( O  =  <  (5.*Q(K) )**.259  l+TSAT <K >  +  l 0. 

120  TF  <K  )  =  <  TSB  <0+TB(K  >  )/2. 

FMU(O  =  .045/TF  (O  -  (4.5E-0S  )  *TF ( K  >  +  { 2 .E- 1  1  )  *  (  TF  (0**2  ) 

PR  (O  =  480  ./TF  (K)-(2.E-03)*TF (K>+  ( 7.E-09  )  *  (TF  (0**3  ) 

BARH (K ) *5 1 34  .*V (K  >*  <  (FMU(K ) / { DH IK) *62 .*V (K )  )  )**.2 )/ (PW (K  >**.667  > 
QTEST  =  BARH  <K  >*( TSB  (K)-TB(K)  )  +  .2*(  (  TSB  (K  l-TSAT  (O  - 1  0  .  >  **3 .86  ) 

IF  (QTEST-O(O)  29.29.91 

91  IF  (LOT-LOP)  92.99.92 

92  TSB  (0=TSB<  0-20. 

GO  TO  120 

29  TSB  <0=TSB<0+5. 

LOP=LOT 
GO  TO  120 

99  BARH(0=0(K)/(  TS6(0-TB(K>  > 

OTSUB=TSAT (K )-TB(K 1 

QBO(K)  =  (4.E  +  05)*( 1  .  +  .021 7*DTSUB  )  * ( 1 .+30.* 

1  ( PRE SS  ( K  1/3200.  1-31  .0*  (  (PRESS  (<  >/3200.  >**1.33)1*  (V(0**.333  > 

QB2 (K 1=7000.* (SQRTF (V (O ) >*DTSUB 

QB3(O  =  (4.E  +  05>*<  l  •  +  •  02*DTSU0  )*(V(K)**»333) 

TS(K  >=TSB(0 
25  CONTINUE 
PUNCH  108 
DO  10  I  L“  1 . N 

10  PUNCH  109.  1L.  To ( I L 1 «  VS(lL).  TSAT(IL).  VlIL>«  PRESS<1L> 

PUNCH  107 
PUNCH  1 1 1 
DO  35  IK« 1 «N 

35  PUNCH  110.  IK.  0 (  IK).  BARH (IK).  UBO(IK).  QB2  (IK). QB3 (IK) 

40  CONTINUE 
END 


143 


AEDC.TDR. 63-58 


H  C  ROLAND  *  FILM  COEFF  -  .30DIA  -  60  8/IN2SEC 
INPUT  DATA 

N  FL.RATE  TCOOL  PRIN  AIN  DHIN 

20  20,00000  S5.0C000  2000,000  ,57500  ,25000 

RIB, SECS  NO,  RIBS  RIB  THCKNS  DCOOL 


6  20,00000  .02000  ,20000 


SEC 

R 1 

D2 

TH 

DX 

Q 

1 

,70000 

•04B74 

.80000 

.10000 

30000000. 

2 

,60000 

.05725 

•8COOO 

•  10000 

30000000. 

3 

.50000 

.06917 

.80000 

•  1 0000 

30000000. 

4 

.40000 

.08687 

.80000 

.10000 

30000000. 

5 

.30000 

. 1 1 520 

.80000 

.10000 

30000000. 

6 

.20000 

.16449 

.80000 

•  1 0000 

30000000. 

7 

.15000 

.14812 

.50000 

.05000 

30000000. 

8 

.12500 

.15685 

•  1  0000 

.05000 

30000000. 

<9 

•  1  2500 

.15685 

.  1  COOO 

.05000 

30000000. 

10 

.16250 

.13797 

.  1  4000 

.50000 

25000000. 

1  1 

.22500 

. 1 1329 

.14000 

.50000 

20000000. 

12 

.28750 

.09517 

.14000 

•50000 

18000000. 

13 

.35000 

.08157 

.14000 

.50000 

1 6000000 • 

14 

.41250 

.07113 

.  1  4n00 

.50000 

1 4000000 . 

15 

.47500 

.06292 

•  1  4000 

.50000 

1 2000000. 

16 

.53750 

.05632 

.14000 

.50000 

1 0000000. 

1  7 

.60000 

.05093 

.14000 

.50000 

8000000. 

18 

•66250 

.04645 

•  1  4000 

.50000 

7000000. 

19 

.72500 

.04268 

•  1 4000 

.50000 

6000000. 

20 

.78750 

.03946 

.14000 

.50000 

6000000. 
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OUTPUT 


N 

fl.rate 

TCOOL 

PR  IN 

AIN 

DH1N 

20 

20.00000 

85.00000 

632.665 

.57500 

.25000 

SEC 

T-BULK 

T-SURF 

t-sat 

VEL. 

PRESS 

1 

85.9 

485.8 

404.7 

231  .9 

238.6 

2 

87.6 

487.4 

401.3 

232.0 

230.6 

3 

89.0 

489.4 

398.3 

232.0 

224.0 

4 

90.2 

492.0 

395.9 

232.1 

218.8 

5 

91 .1 

495.2 

394.1 

232.1 

214.6 

6 

91  .7 

493.8 

392.7 

232.1 

211.9 

7 

92.1 

496.0 

391.9 

232.2 

210.3 

8 

92.2 

497.7 

391  .6 

232.2 

209.6 

9 

92.3 

497.4 

391.3 

232.2 

209.0 

1  0 

93.0 

479.3 

389.3 

232.2 

205.0 

1  1 

94.3 

448.0 

385.0 

232.3 

196.7 

12 

95.8 

424,2 

379.4 

232.3 

186.3 

1  3 

97.4 

393.5 

372.1 

232.4 

1  73.5 

14 

99.2 

362.9 

362.9 

232.5 

158.3 

15 

100.9 

331.3 

351.3 

232.6 

140.6 

16 

102.6 

301.7 

336.7 

232.6 

120.2 

1  7 

104. 2 

268.6 

318.6 

232.7 

97.0 

18 

105.6 

250.7 

295.7 

232.8 

70.9 

19 

107.0 

235.0 

265.0 

232.8 

41  .9 

20 

1G8.4 

229.5 

198.0 

232.9 

9.9 

I 
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r 

4 

s 

1 

OUTPUT 


SEC 

Q 

H 

0801 

QB02 

QB03 

1 

30000000. 

75015.8 

43819626. 

33989235. 

18095875. 

2 

30000000. 

75035.4 

42606443. 

33449912. 

17847155. 

3 

30000000. 

74914.7 

41597287. 

32969585. 

17634880. 

4 

30000000. 

74649.6 

40785714. 

3261 1449. 

17460517. 

5 

30000C00 . 

74242.7 

40161026. 

32315738. 

17324169. 

6 

30000000. 

74620.9 

39706804. 

32099183. 

17224327. 

7 

30000000. 

73905.3 

39461055. 

31982089. 

17170344. 

8 

30000000. 

73991 .5 

39354934 • 

31932026. 

17147269. 

9 

30000000 . 

74068.2 

39260061 • 

31887560. 

17126774. 

10 

25000000. 

6471 1 .3 

38658761 • 

31606412. 

16997208. 

1  1 

2000C000. 

56541.2 

37403280. 

31014405. 

16724433. 

12 

18000000. 

54803.4 

35832543. 

30261063. 

16377399. 

13 

160C0000. 

54037.8 

33912359. 

2931521 1 • 

15941776. 

14 

14000000. 

53089.2 

0  . 

0. 

0. 

15 

1  2000000 . 

52100.7 

0. 

0. 

0. 

16 

10000000. 

50227.9 

0. 

0. 

0. 

17 

8000000. 

48665.2 

0. 

0. 

0. 

18 

70C0000  . 

48267.0 

0  • 

0. 

0. 

19 

6000000. 

46888.9 

0. 

0. 

0. 

20 

6000000 . 

49560.9 

7809965. 

9578363. 

6861783. 
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H  C  ROLAND  *  FILM  COEFF  -  .30DIA  -  60  8/IN2SEC 
INPUT  DATA 

N  FL.RATE  TCOOL  PRIN  AIN  DHIN 

20  20.00000  85.00000  2000.000  .57500  .25000 

RIB. SECS  NO.  RIBS  RIB  THCKNS  DCOOL 

6  20.00000  .02000  .20000 


EC 

R 1 

D2 

TH 

DX 

0 

1 

•70000 

.04874 

.80000 

. 1 0000 

30000000. 

2 

. 60000 

.05725 

.80000 

.10000 

30000000. 

3 

.50000 

.06917 

.80000 

.10000 

3000000C • 

4 

•  4  0  C  C  0 

.08687 

.80000 

•  1 0000 

30000000. 

5 

.30000 

•  1 1 52C 

.80000 

•  1 0000 

30000000. 

6 

.20000 

.16449 

.80000 

. 1 oooo 

30000000. 

7 

.15000 

.14812 

.50000 

.05000 

30000000. 

8 

•  1  2500 

.  1  5685 

.10000 

.05000 

30000000 • 

9 

.12500 

•  1 5685 

•  1 0000 

.05000 

30000000. 

10 

.16250 

. 1 3797 

•  1 4000 

.50000 

25000000. 

1  1 

.22500 

. 1 1 329 

•  14000 

.50000 

20000000. 

12 

.28750 

.09517 

.14000 

.50000 

18000000. 

13 

.35000 

.08157 

. 1 4000 

•50CCC 

16000000. 

14 

.41250 

.071  13 

.  1  400G 

.50000 

14000000. 

15 

.47500 

.06292 

.14000 

.50000 

1  2000000. 

16 

.53750 

.05632 

. 1 4000 

.50000 

10000000. 

17 

.60000 

.05093 

.  1  4000 

.50000 

8000000. 

18 

.66250 

.04645 

.14000 

.50000 

7000000. 

19 

.72500 

.04268 

.14000 

.50000 

6000000. 

20 

.78750 

.03946 

. 14000 

.50000 

6000000. 
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I 

ir 


OUTPUT 


N 

fl.rate 

T  COOL 

20 

20*00000 

85.00000 

SEC 

T-BULK 

T-SURF 

1 

85.9 

485.8 

2 

87.6 

487.4 

3 

89.0 

489.4 

4 

90.2 

492.0 

5 

91 .1 

495.2 

6 

91  .7 

493.8 

7 

92.1 

498.0 

8 

92.2 

497.7 

9 

92.3 

497.4 

1  0 

93.0 

479.3 

1  1 

94.3 

448.0 

12 

95.8 

424.2 

1  3 

97.4 

393.5 

14 

99.2 

362.9 

15 

IOC. 9 

331  .3 

16 

102.6 

301.7 

17 

104.2 

268.6 

18 

1  05.6 

250.7 

19 

107.0 

235.0 

20 

1C8.4 

229,5 

PR  IN 

AIN 

DHIN 

632.665 

.57500 

.25000 

t-sat 

VEL* 

PRESS 

404.7 

231  .9 

238.6 

401.3 

232.0 

230.6 

398.3 

232.0 

224.0 

395.9 

232.1 

218.8 

394.1 

232.1 

214.8 

392.7 

232.1 

211.9 

391  .9 

232.2 

210.3 

391  .6 

232.2 

209.6 

391.3 

232.2 

209.0 

389.3 

232.2 

205.0 

385.0 

232.3 

196.7 

379.4 

232.3 

186.3 

372.1 

232.4 

173.5 

362.9 

232.5 

158.3 

351.3 

232.6 

140.6 

336.7 

232.6 

120.2 

318.6 

232.7 

97.0 

295.7 

232.8 

70.9 

265.0 

232.8 

41  .9 

198.0 

232.9 

9.9 

148 


AEDC-TDR. 63-58 


OUTPUT 


SEC  Q 

H 

0801 

QB02 

QB03 

1 

30000000 • 

75015.8 

43819626. 

33989235. 

18095875. 

2 

30000000. 

75035.4 

42606443. 

33449912. 

17847155. 

3 

30000000. 

74914.7 

41597287. 

32989585. 

17634880. 

4 

30000000. 

74649.6 

40785714. 

3261 1449. 

17460517. 

5 

30000000. 

74242.7 

40161026. 

32315738. 

17324169. 

6 

300C0000 . 

74620.9 

39706804. 

32099183. 

17224327. 

7 

300C0000 • 

73905.3 

39461055. 

31982089. 

17170344. 

8 

30000000. 

73991.5 

39354934. 

31932026. 

17147269. 

9 

30000C00 . 

74068.2 

39260061 • 

31887560. 

17126774. 

10 

25000000. 

6471  1  .3 

38658761 • 

31606412. 

16997208. 

1  1 

20000000. 

56541.2 

37403280. 

31014405. 

16724433. 

12 

18000000. 

54803.4 

3583254 3. 

30261063. 

16377399. 

13 

16000000. 

54037.8 

33912359. 

2931521 1  • 

15941776. 

14 

1 40C0000 . 

53089.2 

0. 

0. 

0. 

1  5 

12000000. 

52100.7 

0. 

0. 

0. 

16 

10000000. 

50227.9 

0. 

0. 

0. 

1  7 

8000000. 

48665.2 

0. 

0. 

0. 

18 

7000000. 

48267.0 

0. 

0. 

0. 

19 

6000000. 

46888.9 

0. 

0. 

0. 

20 

60C0000. 

49560.9 

7809965. 

9578363. 

6861783. 
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! 


Si 


$ 

! 


f 

) 

I 

? 


H  C  ROLAND  *  FILM 

COEFF  - 

. 30D I  A  -  60 

B/IN2SEC 

INPUT 

DATA 

DH1N 

N 

fl.rate 

TCOOL 

PR  IN 

AIN 

20 

20.00000 

85.00000 

310.000 

.57500 

.25000 

RIB. SECS  NO.  RIBS  RIB  THCKNS  DCOOL 

6 

20.00000  .02000  .35000 

SEC 

R 1 

D2 

TH 

DX 

Q 

1 

.70000 

.08374 

.80000 

•  1  0000 

30000000. 

2 

.60000 

.09771 

.80000 

•  10000 

30000000. 

3 

.50000 

•  1 1683 

.80000 

•  10000 

30000000. 

4 

.40000 

•  14418 

.80000 

•  10000 

30000000. 

5 

.30000 

.18525 

.80000 

•  1 0000 

30000000. 

6 

.20000 

•24962 

.80000 

.  1  coco 

30000000  * 

7 

.15000 

.22433 

.50000 

.05000 

juoOOOOC . 

B 

.12500 

.23187 

.10000 

.05000 

30000000. 

9 

.12500 

.23187 

.  1 0000 

•05000 

30000000. 

10 

.16250 

.20938 

.  1 4000 

.50000 

25000000. 

1  1 

.22500 

.17798 

•  14000 

.50000 

20000000. 

12 

.28750 

.15335 

.14000 

.50000 

18000000. 

1  3 

.35000 

.13388 

. 1 4000 

.50000 

1 6000000. 

14 

.41250 

.  1  1830 

.14000 

.50000 

1 4000000. 

15 

.47500 

.10568 

.14000 

.50000 

1 20C0000 • 

1  6 

.53750 

.09531 

.  1 4000 

.50000 

! 0000000. 

1  7 

.60000 

.08668 

.  1 4000 

•50CC0 

8000000. 

1  B 

.66250 

.07941 

.14000 

•50000 

7000000. 

19 

.72500 

.07321 

. 1 40C0 

.50000 

6000000 • 

20 

.78750 

•C6787 

.14000 

. 50000 

6000000 • 
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OUTPUT 


N 

FL.RATE 

TCOOL 

20 

20.00000 

85.0C000 

SEC 

T-BULK 

T-SURF 

1 

85.9 

486.3 

2 

87.6 

485.3 

3 

89.0 

484.5 

4 

90.2 

488.9 

5 

91 .1 

488.4 

6 

91  .7 

488.0 

7 

92.1 

487.8 

8 

92.2 

487.7 

<5 

92.3 

487.6 

1  0 

93.0 

481.0 

1  1 

94.3 

467.8 

12 

95.8 

458.2 

1  3 

97.4 

448.0 

14 

99.2 

437.2 

1  5 

100.9 

420.6 

16 

102.6 

393.1 

17 

104.2 

351.3 

18 

1  05.6 

327.8 

19 

107.0 

298.7 

2C 

1C8.4 

299.2 

PR  IN 

AIN 

DH1N 

31C.000 

.57500 

.25000 

T-SAT 

VEL* 

PRESS 

365.2 

132.5 

162.0 

364.2 

132.5 

160.5 

363.4 

132.6 

159.2 

362.8 

132.6 

158.2 

362.3 

132.6 

157.4 

361  .9 

132.6 

156.7 

361.7 

132.6 

156.4 

361  .6 

132.6 

156.2 

361  .5 

132.6 

156.1 

360.9 

132.7 

155.2 

359.7 

132.7 

153.4 

358.3 

132.7 

151.2 

356.6 

132.8 

148.6 

354.6 

132.8 

145.5 

352.2 

132.9 

142.0 

349.5 

132.9 

138.0 

346.3 

133.0 

133.5 

342.8 

133.0 

128.4 

338.7 

133.0 

122.9 

334.2 

133.1 

1  16.7 
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OUTPUT 


SEC 

Q 

H 

QBO  1 

QB02 

QB03 

1 

30000000. 

74921  .4 

27787656. 

22510394. 

13410820. 

2 

30000000. 

75422.0 

274626J2. 

22299496. 

13303899. 

3 

3000000C. 

75848.3 

27191553. 

22 1 22063. 

13213947. 

4 

30000000. 

75244.0 

26971956. 

21977381 . 

13140600. 

5 

30000000. 

75513.3 

26bC0979. 

2186449 d. 

13083375. 

6 

30000G0C . 

7571 1 .1 

26675194. 

21782134. 

13041625. 

7 

30Q00CC0. 

75817.6 

266C664 1 • 

21737909. 

1301921 1 • 

8 

30000000. 

75862.6 

26576934. 

21719244. 

13009752. 

9 

3000000C . 

75902.1 

26550506. 

21702876. 

13001460. 

1  0 

25000000. 

64441 .0 

26389280 . 

21603645. 

12951 186. 

1  1 

200000CC • 

53554.4 

26066900. 

21406474. 

12B51 301 . 

12 

1 80000C0 . 

49672. 8 

25686454. 

2 1 1 75236. 

12734180. 

1  3 

16000000. 

45640,9 

25244280. 

20907643. 

12598662. 

14 

1 4000000. 

41418.7 

24743665. 

20607438. 

12446666. 

15 

12000000. 

37537.7 

24188255. 

20278402. 

122801 15. 

16 

10000000. 

34425. b 

23581916. 

19924250. 

12100904. 

1  7 

8000000. 

32368.6 

0. 

0. 

0. 

18 

7000000 . 

31512.0 

0. 

0. 

0. 

19 

6000 OCC . 

31295.6 

0* 

0. 

0. 

20 

60000C0 . 

31450.6 

0. 

0. 

0. 
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H  C  ROLAND  *  FILM 

COEFF  - 

•30DIA  - 

60  B/IN2SEC 

INPUT 

DATA 

N 

fl.rate 

TCOOL 

PR  IN 

AIN 

DH1N 

20 

25.00000 

85.00000 

3000.000 

.57500 

.25000 

RIB. SECS  NO.  RIBS  RIB  THCKNS 

DCOOL 

6 

20.00000  .02000  • 

16000 

SEC 

R1 

02 

TH 

DX 

Q 

1 

.70000 

.03919 

.80000 

.10000 

30000000. 

2 

.60000 

.04612 

.80000 

.10000 

30000000. 

3 

.50000 

.05589 

.80000 

• icooo 

30000000. 

4 

.40000 

.07059 

.80000 

•  1 cooo 

30000000. 

5 

.30000 

.09461 

.80000 

.10000 

30000000. 

6 

.20000 

•  13814 

.80000 

.10000 

30000000. 

7 

. 15000 

.  1  2450 

.50000 

.05000 

300000C0. 

8 

•  1 2500 

•  1  3320 

.10000 

.05000 

30000000. 

9 

.  12500 

.13320 

. 1 0000 

.05000 

30000000. 

10 

. 16250 

•  1  1  587 

. 1 4000 

.50000 

25000000. 

1  1 

.22500 

.09385 

.14000 

.50000 

20000000. 

12 

.26750 

.0781  1 

•  1  4000 

.50000 

1 8000000. 

13 

.35000 

.06653 

.  1  4000 

.50000 

16000000. 

14 

.41250 

.05775 

. 1 4000 

.50000 

14000000. 

15 

.47500 

.05093 

. 14000 

.50000 

12000000. 

16 

.53750 

.04549 

.14000 

.50000 

10000000. 

17 

.60000 

.04107 

.14000 

.50000 

8000000. 

18 

.66250 

.03741 

.14000 

.50000 

7000000. 

19 

.72500 

.03433 

.14000 

.50000 

6000000. 

20 

.78750 

.03172 

.14000 

.50000 

6000000. 

AEDC.TDR. 63.58 


OUTPUT 


N 

fl.rate 

TCOOL 

PR  IN 

AIN 

DH 

20 

25.00000 

85.00000 

3000.000 

.57500 

.251 

SEC 

T— BULK 

T-SURF 

t-sat 

VEL* 

PRESS 

1  ■ 

85.7 

414.3 

639.3 

362.4 

2064.6 

2 

87.0 

422.4 

637.4 

362 . 5 

2041 .7 

3 

88.2 

430.8 

635.8 

362.5 

2022.9 

4 

89.1 

439.6 

634.6 

362.6 

2008.1 

5 

89.9 

448.7 

633.7 

362.6 

1996.8 

6 

90.4 

463.0 

633.0 

362.7 

1988.7 

7 

90.6 

472.7 

632.7 

362.7 

1984.4 

8 

90  .8 

477.5 

632.5 

362.7 

1982.5 

9 

90.8 

477.4 

632.4 

362.7 

1980.8 

10 

91  .4 

426.5 

631  .5 

362.7 

1970.1 

1 1 

92.4 

369.6 

629.6 

362.8 

1947.2 

12 

93.6 

347.2 

627.2 

362.9 

1918.1 

13 

94.9 

319.3 

624.3 

363.0 

1882.3 

14 

96.3 

295.8 

620.8 

363.1 

1839.2 

1  5 

97.7 

271  .6 

616.6 

363.2 

1 788.5 

16 

99.1 

246.8 

61  1  .8 

363.3 

1730.0 

17 

1  00.3 

221.3 

606.3 

363.4 

1 663.3 

18 

101.5 

205.2 

600.2 

363.4 

1 588.2 

19 

1  02.6 

193.3 

593.3 

363.5 

1504.5 

20 

103.7 

195.7 

585.7 

363.6 

1412.0 

154 


AEOC.TDR.63-58 


OUTPUT 


SEC 

Q 

H 

QB01 

QB02 

QB03 

1 

30000000. 

91304.7 

0. 

0. 

0. 

2 

30000000. 

89462.9 

0. 

0. 

0. 

3 

30000000. 

87553.3 

0. 

0. 

0. 

4 

30000000. 

85594.3 

0. 

0. 

0. 

5 

30000000. 

83603.8 

0. 

0. 

0. 

6 

30000000. 

80504.4 

0. 

0. 

0. 

7 

30C00C00 . 

78529.1 

0. 

0. 

0. 

8 

30000000. 

77566.6 

0. 

0. 

0. 

9 

30000000. 

77612.6 

0. 

0. 

o. 

10 

25000000. 

74602.8 

0. 

0. 

0. 

1  1 

20000000. 

72153.7 

0. 

0. 

0. 

12 

18000000. 

70973.0 

0. 

0. 

0. 

1  3 

16000000. 

71314.5 

0. 

0. 

0. 

14 

14000000. 

70198.9 

0. 

0. 

0. 

1  5 

1  2000C00 . 

69015.4 

0. 

0. 

0. 

16 

10000000. 

67696.4 

0. 

0. 

0. 

1  7 

8000000. 

6611 3.3 

0. 

0. 

0. 

18 

7000000. 

67507.4 

0. 

0. 

0. 

19 

6000000. 

66131.4 

0. 

0. 

0. 

20 

600000C  . 

65196.0 

0. 

0. 

o. 
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